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COMPUTER PROGRAM FOR PRELIMINARY DESIGN AND 


ANALYSIS OF V/STOL TIP-TURBINE FANS 
by Henry C. Haller, Seymour Lieblein, and Bruce M. Auer 
Lewis Research Center 

SUMMARY 

An analysis of a tip-turbine driven lift fan was developed for use in preliminary de- 
sign or parametric investigation of tip- turbine lift fan systems applicable to V/STOL 
aircraft. An electronic digital computer program was developed to determine flow con- 
ditions and flow areas for a specific total thrust requirement over a range of fan stage 
and of turbine input variables. 

The program developed has provisions for treating single-stage fan rotors with 
either one- or two-stage tip turbines mounted on the fan periphery. Both turbine types 
have uncooled impulse rotors with the single-stage designs employing impulse, diffusing, 
or no exit stators. The program also treats tip-turbine velocity diagram determination, 
blade row loss analysis, partial admission (for the single-stage turbine), and considera- 
tion of turbine-to-fan leakage. The fan-stage inputs required are pressure ratio, effi- 
ciency, and inlet axial Mach number. 

Because of the detailed nature of the turbine loss analysis and provision for several 
stage types, the program is capable of handling turbine flows over a wide range of pres- 
sure and temperature levels. This allowed for solutions for a number of generator 
types or supply sources for the turbine flow. 

Sample calculations were made for a high-pressure turbine supply case to exhibit 
the potential of the program. 


INTRODUCTION 

The potential of the tip-turbine lift-fan concept for V/STOL applications has been 
investigated in numerous research efforts (refs. 1 to 7) and demonstrated by the XV- 5A 
fan-in-wii^ research aircraft (refs. 8 and 9). In the tip-turbine fan concept, turbine 
blades are mounted around the periphery of the rotor of a single-stage ducted fan to 



provide a compact drive arrangement. Lift fans employ relatively low fan-stage pres- 
sure ratios in order to obtain low discharge noise and low specific fuel consumption. 

Lift propulsion systems with tip-turbine lift fans can be used in the design of short- 
haul V/STOL transports (e.g. , refs. 10 and 11). An airplane of this type might employ 
multiple fans in the wings or in pods near the fuselage or out along the wings. The 
principal characteristic of the tip-turbine lift fan is its relatively short axial depth, 
which could constitute a favorable advantage in V/STOL aircraft design. 

Tip-turbine lift fans have been powered by the hot exhaust from a gas generator 
(turbojet engine) as in the case of the XV- 5A airplane. However, another drive possi- 
bility is a compressed-air generator powerplant that would supply relatively cold high- 
pressure air to the tip turbine. This concept could include burners ahead of the turbine. 
The compressed air drive approach might be considered if it is desired to minimize the 
size and weight of the fan turbine scroll or if large amounts of ducting are required be- 
tween the fans and the powerplants (e.g. , for engine-out safety). 

The heart of the tip-turbine lift fan is a single-stage ducted-fan rotor with turbine 
blades attached to the periphery. Other major components include the fan and turbine 
exit stators, the support frames, the inlet bellmouth, a scroll, which distributes the 
supply flow to the tip turbine nozzles, the exit louvers, which provide a means for vary- 
ing the fan exhaust flow direction, and a burner, if specified. 

This analysis is a cycle study that allows a determination of the gross characteris- 
tics of tip-turbine lift fans for use in preliminary design or for parametric analytical 
studies of lift-fan systems. In the analysis developed herein, overall fan-stage per- 
formance parameters and turbine supply flow state are program inputs. Fan component 
size and required turbine flow rate for a prescribed total fan thrust are the major out- 
puts. The tip turbine is treated in detail with an analysis that includes velocity diagram 
determination, blade row loss analysis, turbine-to-fan leakage, and scroll duct pres- 
sure losses. Four types of tip-turbine stage arrangements are investigated: a conven- 
tional single-stage turbine with impulse exit stators, a single-stage turbine with diffus- 
ing exit stators, a single-stage turbine without exit stators, and a two-stage turbine with 
impulse exit stators. In all four types of turbines, the rotors are designed for impulse 
operation. 

This report presents a description of the system analysis and a working computer 
program that includes input and output listings. A sample calculation is also included to 
point out program operating characteristics, and a comparison of the four turbine stage 
types is presented to illustrate the program output variations. 
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ANALYSIS 


Lift- Fan Model 

A schematic drawing of the basic elements of the lift-fan system showing the fan 
stages and the turbine-flow-path station des^nations is shown in figure 1. The config- 
uration illustrated is an outlet-stator fan stage with a single-stage turbine mounted on 
the periphery of the fan rotor. An inlet duct and a burner (not always required) are also 
shown in the figure. Compressed air or hot gas at total conditions p^ and T’j enters 
the duct and passes directly, or through a burner, to the entrance of the turbine scroll, 
which distributes the flow around the periphery to the turbine inlet stator. The flow is 
accelerated in the inlet stator and directed against the turbine rotor blades. Work is 
extracted from the rotor to drive the fan. The flow is then discharged either directly or 
through the exit stators, with the axial component of outlet velocity contributing to the 
overall thrust of the tip- turbine fan system. 

Ambient air at p^Q and TJ^q enters the fan rotor with velocity Vjj. The rotor 
produces an increase in static and total pressure from inlet to exit to accelerate the fan 
air flow to velocity Vj g, thus providing the majority of the system thrust. Louvers are 
considered at the discharge of the fan for thrust vectoring. 

A schematic cross sectional drawing of a tip-turbine lift-fan configuration as con- 
sidered in reference 6 is shown in figure 2. This figure also shows a planform view of 
the fan system. Symbols designating all quantities used in the description of the tip- 
turbine fan are given in appendix A. 


Approach 

The primary objective of the computer program development described herein was 
to afford a means of obtaining preliminary design characteristics of tip-turbine lift fans 
that could serve as thrusters for V/STOL applications. The program is capable of de- 
signing a tip-turbine fan to provide a prescribed total thrust level, fan plus turbine, with 
specified fan-stage inputs such as overall efficiency and pressure ratio, inlet axial Mach 
number, rotor tip speed, and rotor inlet hub-tip ratio. It was presumed that these inputs 
would be obtained from detailed fan- stage aerodynamic design procedures. 

Basic turbine inputs include the supply gas total temperature and total pressure, 
turbine stage arrangement, and specification of burning. Several types of turbine stage 
design were included to allow for good performance with respect to such factors as pres- 
sure level of turbine supply flow, low discharge noise, high stage efficiency, and reduced 
turbine-to-fan leakage. Turbine stage efficiency was determined in the calculation from 
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the resultant blade row velocity diagrams (mean radius values) in conjunction with pre- 
scribed blade row loss variations for stators (nozzle), rotors (buckets), and exit sta- 
tors. Loss relations for the flow in the scroll and for gas leakage from the turbine to 
the outer flow of the fan upstream of the fan rotor were also included (see fig. 3). Fur- 
ther provision for the effects of losses was made by the inclusion of loss inputs for the 
turbine supply flow duct, the fan inlet, the fan and turbine outlet ducts, and the exit 
louvers. Burner temperature and pressure drop were inputs. 

The principal outputs of the calculation program were the gas flow rate into the 
turbine required to produce the desired total thrust, the hub and tip diameters at the 
various fan and turbine stations, the diameter of the scroll inlet, and the maximum 
overall diameter of the unit. 


Turbine Stage Design 

Four turbine stage designs based on impulse rotor operation were included in the 
analysis. Impulse rotor design (constant static pressure across rotor) was specified 
in all cases for simplicity and for minimum inlet stator exit static pressure. Low inlet 
stator discharge pressure is desired in order to minimize gas leakage from the turbine 
flow path to the fan flow path in the axial clearance in front of and behind the fan-turbine 
rotor. 

The four-turbine- stage designs included are (1) a single-stage turbine with impulse 
exit stators, (2) a single-stage turbine with diffusing exit stators, (3) a single-stage 
turbine without exit stators, and (4) a two-stage turbine with impulse exit stators. A 
schematic diagram of the one- and two-stage tip turbine concepts is shown in figure 3. 
The figure also shows the potential leakage flow paths. 

A diffusing exit stator was considered for the single-stage turbine to provide for 
reduced discharge velocity and noise. The diffusing exit stator also results in reduced 
inlet stator exit pressure with corresponding increase in turbine specific work output 
and with reduced turbine-to-fan leakage. However, these gains are obtained at the ex- 
pense of reduced turbine discharge thrust and turbine efficiency. The stage with no exit 
stator was included to determine configurations with the minimum number of blade rows 
and minimum blade-row losses. Here, too. the objective is achieved at the expense of 
turbine discharge thrust. The two-stage configuration was introduced to maintain high 
overall efficiency levels and low discharge velocity (low noise) for cases with supply gas 
at high pressures. 

Partial admission was provided in the development for the single-stage cases as a 
means for increasing turbine rotor blade height, if desired, or for simulating effects of 
localized flow reduction for purposes of thrust control. 


4 



Development of Equations 


The development of the equations describii^ the design of the tip-turbine lift-fan 
system is presented in detail in appendixes B to G. The main program, which gives the 
overall fan design and calls the various specific turbine subroutines, is presented in 
appendix B. 

The main program calculation procedure yields an overall fan design, that is, fan 
component diameters and station thermodynamics. All analyses, resulting design 
equations, and program logic are given in detail in appendix B. The main program also 
outputs fan performance parameters such as thrust specific fuel consumption of the bur- 
ner (if used), bypass ratio (ratio of fan flow to turbine flow), and fan thrust augmenta- 
tion ratio (ratio of total thrust to the thrust that would be obtained from isentropic ex- 
pansion to ambient of the supply gas flow at p^ and Tp. 

Appendix C gives the procedures for determinii^ turbine losses. The losses are 
considered to comprise profile losses, secondary flow losses, and shock losses. Ap- 
pendix C details the approach to the determination of these losses by use of correlations 
obtained from the available literature as reported by numerous investigators (refs. 12 
to 14). 

The single-stage turbine design subroutines are given in appendix D. This appendix 
includes the detailed analysis of the three single-stage turbine types: impulse exit sta- 
tors, diffusing exit stators, and no exit stators. The loss analysis for diffusing exit sta- 
tors is also given in appendix D. 

The two-stage turbine analysis is presented in appendix E, and the leakage analysis 
required by all four turbine arrangements is given by appendix F. The various fluid 
property and computational subroutines are listed in detail in appendix G. A complete 
program listing is included and given in appendix H. 


COMPUTATIONAL PROCEDURE 

Calculations were performed using an iterative procedure programmed into an elec- 
tronic digital computer. The calculations were based on the relations and procedures 
described in the appendixes of this report. The calculation program comprises iteration 
loops within iteration loops as shown on the iteration flow chart (fig. 4) for the main 
computational routine. 

The program as written follows the following basic steps for obtaining the desired 
thrust for a tip-turbine lift fan; 

(1) Read in turbine flow properties, trial efficiency, initial value of flow rate, and 
exit stator exit Mach number. 
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(2) Read in fan-stage pressure ratio and efficiency, inlet axial Mach number, and 
rotor tip speed. 

(3) Calculate the shaft power and amount of turbine thrust. 

(4) Calculate fan thrust, fan size, and turbine blade height. (Note that the total 
thrust (fan plus turbine) does not equal the desired thrust at this point. ) 

(5) Calculate turbine losses, leakage, and revised efficiency. 

(6) Recalculate fan geometry and total thrust (note that total thrust does not equal 
desired thrust). 

(7) Repeat steps (1) to (6) using a second input value of turbine flow rate. 

(8) Use a linear interpolation of turbine weight flow and total thrust to obtain the re- 
quired turbine flow rate for the specified thrust. 

(9) Repeat steps (1) to (6) to obtain the correct total thrust, final fan turbine geom- 
etry, and cycle thermodynamics. 

The calculated thrust and desired thrust have been shown to agree within 1/2 percent 
based on this calculation procedure. 

The entire procedure is repeated for a range of input values of the turbine stator 
exit axial Mach number Mg. The program then has the option to print the results ob- 
tained for each input value of Mg or enter an optimization routine. 

The optimization routine first plots turbine output values of inlet stator exit angle 
0 !g, turbine tip diameter d^ rotor inlet relative Mach number M^ turbine rotor 
blade height H.^, turbine weight flow required w^^, and turbine louver exit axial velocity 
Vj^, all against exit stator exit Mach number Mg. Over the range of Mg used, the 
program then obtains a value of Mg that yields minimum turbine flow rate Wj. Before 
this value is accepted, a check is made on a g to determine whether the angle is less 
than the limiting value (taken as <75° or 1. 31 rad). Angles larger than 75° result 
in excessive inlet stator losses. The other outputs are visually inspected to insure that 
the turbine rotor blade height Hrp exceeds a minimum acceptable value. 

A sample computer output plot of the variation of these parameters for the single- 
stage turbine without exit stators is given in figure 5. The parameters are normalized 
to the values obtained for Mg = 0. 70 of table I. Similar curves can be obtained for the 
other three turbine arrangements investigated. 

The major inputs required for the main computational routine are listed in appen- 
dix B. Additional inputs required by the turbine subroutines are given in appendixes 
D and E. A complete listing of all program inputs is given in table I, which is a sample 
program input-output listing. The definition of the variables in table I are given in the 
nomenclature definitions of appendix I. 
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SAMPLE CALCULATIONS 


Sample calculations were run for the tip-turbine lift fan with the four turbine stage 
configurations over a range of turbine stator exit Mach number Mg. A total thrust 
(turbine plus fan) of 10 000 pounds was chosen for the representative cases. The tur- 
bines were specified as full admission and were to be within the limits of practical de- 
sign, that is, reasonable velocity diagrams, turbine rotor blade height, stator diffusion, 
and so forth. An available turbine supply flow pressure ratio Pj/Pq of 8. 0 was speci- 
fied along with a supply total temperature TJ^ of 1080° R. Burning was also specified 
to bring the turbine scroll inlet temperature T^^^ up to 1900° R. A list of the re- 
maining major input variables common to all four turbine types are 


Inlet supply duct pressure loss ratio, P^/Pj^ 0.95 

Burner efficiency, 77^^ 0. 98 

Burner total pressure loss fraction, A-p'/pg 0.06 

Fan-stage pressure ratio, P13/P11 1.20 

Fan inlet axial Mach number, 0. 55 

Fan rotor hub-tip ratio, Yp 0, 35 

Air flow supply duct inlet Mach number, Mj 0. 20 

Fan rotor tip speed, Uj, ft/sec 700 

Fan-stage efficiency, tJj, 0. 86 

Scroll Mach number, 0. 30 

The remaining inputs for the four turbine types are given along with a sample of the cor- 
responding output listings in table I which is for the specific case listed previously. 

The sample sets of data for the specific cases investigated were used to exhibit the 
comparative variations in geometry and performance for the four turbine stage configu- 
rations. Figure 6 shows a plot of the variation in required turbine flow rate with stator 
exit axial Mach number. It is seen from the figure that turbine flow rate does not vary 
greatly with stator exit axial Mach number. The figure also indicates the practical de- 
sign limitation imposed on inlet stator exit angle. For the inlet conditions investigated, 
the two- stage turbine yielded the lowest flow rate because of the inherently higher tur- 
bine total efficiency for these inlet conditions. The single-stage turbine with an impulse 
stator is shown as a single point because realistic turbine velocity diagrams (05 < 75°) 
are not obtainable for stator exit Mach numbers less than 0. 80. 

Another important aspect of turbine design is shown in figure 7 which is a plot of 
turbine exit flow velocity against stator exit axial Mach number. Because of noise con- 
siderations, it is desirable to reduce the turbine exit flow velocity to meet the require- 
ment for low discharge noise. Thus, designing a turbine to the maximum nozzle exit 
angle of around 75° yields minimum exit flow velocity in all cases. The two-stage tur- 
bine has the lowest velocity (fig. 7). 
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The velocities for the two-stage and the single-stage with impulse or diffusing exit 
stators all fall on essentially a single line because of the specification of axial dis- 
charge. The slight variation is due to varying efficiency and thus exit temperature. 

The exit velocity for the single-stage turbine without an exit stator is higher than the 
other three turbine designs because of the residual tangential velocity. Also, the h^her 
axial velocity across the rotor is inherent with this configuration. 

The single-stage turbine without exit stators and the single-stage turbine with im- 
pulse exit stators would probably be unacceptable designs for the supply gas pressure 
level because of the high discharge velocities. 

The calculated variation of turbine-tip diameter with exit stator exit Mach number 
Mg for the four turbine arrangements is shown in figure 8. The limit for a - 15^ is 
shown for each case. 


CONCLUDING REMARKS 

The program described herein is capable of handling turbine flows over a wide 
range of pressure and temperature levels because of the detailed nature of the turbine 
loss analysis and provision for several stage types. This capability allows for a wide 
range of generator types or sources for the turbine flow. The potential of the program 
was exhibited by a sample calculation for a high-pressure turbine supply case. Calcu- 
lations performed with the program compared favorably with the tip-turbine fan designs 
for use with an exhaust gas generator system presented in reference 6. Thrust, geom- 
etry and turbine performance agreed very closely with the values of the reference cases 
when specified for equal turbine flow rate. 

It is hoped that this working program will serve as a useful tool in tip-turbine fan 
preliminary design and for the comparison of various lift-fan propulsion systems for 
V/STOL aircraft applications, 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 24, 1970, 

721-03. 
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APPENDIX A 


SYMBOLS 


A 

flow area, in. 

S 

AR 

aspect ratio, H/c 


AUG 

thrust augmentation ratio 

H 

a 

arc length of turbine admis- 

HV 


sion, in. 


b 

exponent in seal leakage 



(eq. (F5)) 


B 

labyrinth tooth width 

h 

BPR 

bypass ratio 

Ah 

CC 

turbine rotor axial velocity 



multiplier 

IPR 

C 

constant 

J 

Cf 

flow coefficient 


cl 

seal clearance, in. 

K 

Cp 

specific heat at constant 

L 


pressure, Btu/(lbm)(°R) 

M 

C,p 

thrust coefficient 

n 

c 

blade chord, in. 

D 

diffusion factor 

P 

d 

diameter, in. 

Ap 

E 

energy loss, Btu/lb 

P 

e 

expansion function 

q 

FF 

flow factor, (w/A)(yRT'/p') 

R 

F 

thrust, Ibf 


f/a 

fuel-to-air ratio 

Re 

G 

fuel-to-air ratio function. 

r 


O' 

+ 

1 1 

SFC 


gravitational constant, 

32. 174 ft/sec^ 

height, in. 

lower heating value of fuel, 
Btu/lbm 

inlet enthalpy of the liquid fuel, 
Btu/lbm 

enthalpy, Btu/lbm 

specific work or enthalpy 
change, Btu/lbm 

ideal incompressible pressure 
recovery (eq. (B26)) 

mechanical equivalent of heat, 
778 ft-lb/Btu 

loss coefficient 

length, in. 

Mach number 

number of (blades, burners, 
etc.) 

power, Btu/sec 

pressure drop, psi 

pressure, psia 

velocity head, p'-p, psia 

gas constant, 1545/molecular 
weight, ft-lbf/(lbm)(°R) 

Reynolds number 

radius, in. 

specific fuel consumption, 
lbm/(hr)(lbf) 
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SPR 

seal pressure ratio 

V 

blade-jet speed ratio, U^/Vj 

s 

blade spacing, in. 

k 

energy loss coefficient, energy 

AT 

temperature drop, °R 


loss to ideal energy ratio 

T 

temperature, °R 

P 

density, Ibm/ft^ 

TRM 

two-stage turbine mean line 

a 

blade row solidity, c/s 


change factor 

w 

entropy function 

U 

blade velocity, ft/sec 


interpolation function 


ideal jet velocity. 

(U 

fan or louver loss coefficient 


y2gJAh.p ft/sec 

Subscripts: 

V 

absolute gas velocity, ft/sec 

A 

scroll arms 

W 

relative gas velocity, ft/sec 

a 

air 

w 

mass flow, Ibm/sec 

B 

turbine rotor blade row 

X 

two-stage turbine work split 

b 

burner 


factor 

c 

fan shroud 

Y 

hub-to-tip diameter ratio 

ch 

check 

Z 

ratio of leakage flow to main 

cl 

clearance 

a 

flow 

absolute gas flow angle mea- 
sured from axial direction. 

D 

d 

exit duct 
inlet duct 


deg 

eq 

equivalent or corrected 

/3 

relative gas flow angle mea- 

ex 

exhaust 


sured from axial direction. 

F 

fan 

y 

6 

deg 

ratio of specific heats 

ratio of inlet total pressure to 
sea level standard pressure. 

f 

m 

h 

fuel 

air or gas generator 
hub 


p'/14.7 

id 

ideal or isentropic 

V 

efficiency 

L 

louver 

e 

turbine admission arc, deg 

leak 

seal leakage 

X 

turbine speed-work parameter. 

M 

mixture 


Ut/av^ 

m 

mean 


gas viscosity, lbm/(ft)(sec) 

max 

maximum 
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N 

turbine stator blade row 

3 

burner exit 

P 

profile 

3sc 

scroll inlet 

pa 

partial admission 

4 

turbine inlet 

r 

relative 

5 

first turbine stator exit 

S 

exit stator exit 

5’ 

first turbine rotor inlet 

s 

sonic 

6 

first turbine rotor exit 

sc 

scroll 

7 

second turbine stator exit 

se 

secondary 

8 

second turbine rotor exit 

sh 

shock 

10 

fan inlet 

st 

stoichiometric 

11 

fan rotor inlet 

T 

turbine 

12 

fan rotor exit 

t 

tip 

13 

fan stator exit 

u 

circumferential component 

14 

fan louver inlet 

X 

axial component 

15 

fan louver exit 

0 

ambient air ' 

Superscripts: 

1 

air or gas generator supply 

t 

stagnation condition 

2 

inlet duct exit 

* 

reference 

2b 

burner inlet 

— 

overall 


APPENDIX B 


FAN-TURBINE MAIN PROGRAM 

The main program produces an overall fan design, that is, fan overall dimensions, 
total thrust, and thermodynamics. The fan- turbine design chosen for this analysis is 
shown in figure 1. This figure gives the station designations for both the tip turbine and 
the fan. For the two-stage tip turbine, figure 1 is modified by including station 7 for the 
interstage stator exit and station 8 for the second-stage rotor exit (fig. 9). 

The Fan-Turbine Design is initiated by assuming the following major inputs: 

(1) Inlet duct mass flow rate w^^ 

(2) Inlet duct inlet total pressure pj^ 

(3) Inlet duct inlet total temperature Tj 

(4) Turbine stage total-to-total efficiency rjrj, 

(5) Fan-stage pressure ratio PJ3/P11 

(6) Fan inlet axial Mach number M^j 

(7) Fan tip speed Up ^ 

(8) Fan-stage adiabatic efficiency 77^ 

The program nomenclature required for the main program and the turbine sub- 
routines are included in appendix I. 

This appendix details all the equations used to define the thermodynamics of the fan 
and turbine streams along with the descriptive equations for the tip turbine fan radial di- 
mensions. The calculational procedure is initiated by first analyzing the turbine stream. 
Then, on obtaining the shaft power output, the fan stream calculations are made. 


Turbine Calculations 

The inlet duct, as shown in figure 1 accepts flow from a gas or air generator source 
and directs it to the turbine at which point power is extracted. The turbine exit flow is 
then passed through an exit stator and louvers and ejected to supply a portion of the tip- 
turbine fan thrust. The turbine flow path is thus composed of a supply duct, burner, 
scroll, turbine, exit stator, and louvers. 

Inlet duct . - The inlet duct accepts flow from a gas or air source at an inlet total 
pressure and temperature, pj and TJ, respectively. A trial value of flow rate Wj is 
required to initiate the calculations to obtain the flow area Ag. 

The total pressure at the exit of the supply duct is obtained from 
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(Bl) 


P2 = Pi 



where p^/Pj is a duct loss factor and an input. With provision for heat loss in the 
supply duct, where AT^ is a program input, 

- (AT^) (B2) 

With no bleed flow from the duct, 

W 2 = Wj (B3) 


The flow factor at the duct exit is 


FF2 = 




1 +ZJl1m2 


(r2+l)/2(y2-l) 


w VrT' 
A p' 


where M 2 is a program input and y 2 is obtained from 


^2 


'p,2 


^P.2 

778.16; 


with the fluid property subroutines (appendix G) used to provide 


C„ Q = function 

P, ^ 



(B4) 


(B5) 


and 


R 2 = function 
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If high pressure air is supplied, 


f 


is set equal to zero. For both the air and hot-gas 


cases the burner efficiency 
sidered in these calculations. 


^'1 

is equal to 0, since the gas generator burner is not con- 
The inlet duct exit area is then obtained from 



(B6) 


Burner . - The diffuser section is between the inlet duct exit and the burner. The 
diffuser section flow is treated as an isentropic flow with the resulting equation for 
burner inlet area dependent on the number of burners n^ and the burner reference 
Mach number both of which are inputs. All losses in the diffuser section are 

lumped in the input burner pressure loss fraction. From continuity, 




(r2+l)/2(y2-i) 


(B7) 


and the diameter for a can type of burner is then 


(B8) 


The burner exit temperature (station 3) is set equal to the desired inlet stator inlet tem- 
perature, thus 


T’ =T^ 

The required burner fuel-to-air ratio is obtained from the expression 


(^3 


HVr^- 1 

(h^ - H°) 


(B9) 


(BIO) 
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where H° is the inlet enthalpy of the liquid fuel, and h^ and h^ are based on 
and Tg, respectively. This equation is given in the fuel-to-air ratio subroutine 
(appendix G) . 

The total mass flow of the gas mixture is then 


w 


3 



(Bll) 


where the burner exit flow state is given by the following expressions: 


^3 = 



where, from the fluid property subroutines. 


C o = function 
p, o 




and 



In equations (Bll), (B13), and (B14), the total fuel-to-air ratio is defined as 


(B12) 


(B13) 


(B14) 



The burner outlet total pressure is obtained from the expression 


P3 = P2 



(B15) 


where burner inlet pressure p^ is obtained from equation (Bl). The value of Ap'/pg 
is an input and is generally obtained from experimental results for specific applicable 
burner designs and is expected to be functions of reference Mach number and 
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temperature ratio Tg/Tg . The flow factor at the burner exit (constant area burner) is 


FF3 = 


^3 V'^3^3 
P3"^2b 


(B16) 


The burner exit Mach number Mg is obtained from FFg and the function FMM sub- 
routine, where Mg and FFg are related by an expression of the form of equation (B4). 

The burner exit also has a transition duct which accelerates the flow and connects 
to the scroll inlet. The friction pressure loss in this duct is also included in the overall 
burner pressure-drop fraction. Thus, for the calculation, Pgg^ P^* If a nonburning 
case is to be analyzed, the burner pressure-drop fraction is set equal to zero, and tem- 
perature T^ is set equal to Tg. This in effect bypasses the burner entirely. 

Scroll . - The flow, on leaving the burner, enters the turbine scroll. The scroll 
could have from one to four duct inlets where the number of duct inlets equals the num- 
ber of burners nj^ with one or two scroll arms per inlet n^, depending on the ducting 
configuration used. Some typical scroll configurations are shown on figure 10. Thus, 
the number of scroll arms n for a full admission turbine scroll is given by the ex- 

o 

press ion 


”sc = ”b(V 


(B17) 


Then, with the number of scroll arms n and the prescribed scroll inlet Mach num- 

s c 

ber Mggp known, the inlet area of a scroll arm is calculated using continuity given as 
(the scroll is assumed to be isothermal, T^ = "^Ssc “ 



(B 18 ) 


The inlet diameter of each scroll arm is then 


“^Ssc 



(B19) 


The loss in total pressure in the scroll is obtained as a function of the inlet dynamic 
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f 

( 

I 

I 

l! head q3gg- Thus, the total pressure at the scroll exit is 

= Pkc - (B20) 

where the loss coefficient K is an input and is approximated by the momentum loss 

s c 

associated with the number of turns the fluid is required to make. For a typical scroll 
that has three fluid turns, a K„„ value of 0. 6 was selected. The dynamic head is de- 
fined in terms of Mach number as 


^3sc " Pssc 


1 - 


yo ~ 1 2 

t o 3sc 


/gArg-l) 


(B21) 


Turbine exit stator . - A knowledge of the thermodynamic state of the gas flow enter- 
ing the turbine stators along with the fact that the turbine is not treated in detail in the 
main program, requires only turbine exit state conditions to determine the ideal specific 
power output of the turbine. Thus, it is necessary to supply the flow thermodynamics at 
the turbine exit stator exit. 

The exit stator exit axial Mach number Mg is a program input, and along with a 

trial value of = 1. 35 (A better value will be obtained when the stator exit total tern- 
o 

perature is determined. ) allows the calculation of the flow factor at the stator exit from 
the expression 


FFg = 




1 + 


y - 1 


(rg+i)/2(rs-l) 


M 


(B22) 


or for Vg = 1.35 


6. 591 Me 


FFg = 


1 + 0. 175 Me 


3. 3571 


(B23) 


Then it is desired to determine the total pressure at the stator exit p^, which then sets 
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the turbine pressure ratio and the total-to-total enthalpy drop across the turbine. 
It is specified first that 


Pd - Pl - Pq 


(B24) 


Next, it is required to obtain the exit Mach number in the duct Mj^. The duct is allowed 
to vary in flow area to provide flow deceleration if desired. Thus, the duct acts as a 
diffuser with the duct loss in total pressure given by the expression 


Pb 


= 1 - 


Ps 




(B25) 


where the loss term Ap^ is defined as the loss in total pressure. Using the data of 
reference 15 for an effective diffuser included angle of 10° (an angle that yields a 
10-percent loss in static pressure recovery) the loss term is taken as 


Apb 

— - = 0. lO(IPR) (B26) 


The ideal pressure recovery IPR is defined in terms of area ratio Aj^/Ag as 


IPR = 



(B27) 


The term Og/Pg equation (B25) is defined in terms of an input value of stator exit 
axial Mach number as 



Ps 




yg/(rg-i) 


(B28) 


Substitution of equations (B26) to (B28) into (B25) yields the expression for duct total 
pressure ratio 
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(B29) 



0. 1 


1 - 


^S/ JL 




Js 


The duct annular area ratio Aj^/Ag is defined for a 10° included angle and approxi- 
mated by a two-dimensional diffuser due to the very high hub-to-tip ratio as 

^ = 1 + 0. 1744/il^ (B30) 

Ag \kJj^ 

where the duct length to height ratio (L/H)jj is an input. The height H is defined at the 
stator exit. If no duct diffusion is specified, the value of (L/H)jj is set equal to zero. 

With Pj)/Pg from equation (B29) and Aj^/Ag from equation (B30), the duct flow 
factor is calculated from 


Ps ^s) 


(B31) 


Once the value of FFp is known, the duct exit Mach number Mj^ is obtained using the 
function FMM subroutine given in appendix G. After the Mach number at the duct exit 
is determined, the total pressure is obtained from the expression (for yg = 1.35); 

/ o \3. 86 

p^= PO^I +0.175 M^j (B32) 

Then, from equations (B29) and (B32), the total pressure at the stator exit is obtained as 


PS = 



(B33) 
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The values of turbine inlet and exit total pressures (p^ and p^) obtained from equa- 
tions (B20) and (B33), respectively, together with the total-to-total turbine efficiency r]rj, 
(an input for the initial calculation and an output of the turbine subroutine on subsequent 
cases) and the inlet total temperature enable the determination of the actual total 
enthalpy drop across the turbine Ah;p and the stator exit total temperature Tg. This 
calculation is performed in the subroutine called ADIABT given in appendix G. 

The gas constant at the stator exit Rg is obtained from the GASCON subroutine 
(appendix G) and at the stator exit is obtained from the GAMMA subroutine (appen- 
dix G) based on the determined value of stator exit total temperature Tg. From the new 
value of 7 g, a value of FFg is recalculated in equation (B22). The flow factor at the 
duct exit is then calculated from the total pressure and duct area ratio given previously 
according to equation (B31). Then, a new value of duct exit Mach number Mj^ is ob- 
tained from the FMM subroutine. The total pressure at the duct exit is then recalculated 
from the expression 


Pb- Po 




(B34) 


D 


where y-^ = 7 g. The stator exit total pressure p^ is recalculated in equation (B33) 
using the new value of pj^ and the previous value of Pj^/Pg from equation (B29). 

With the new value of pg and with p^, T^, and the calculation proceeds to the 
ADIABT subroutine to obtain new values of Ah and Tg. The entire procedure is re- 
peated until the stator exit temperature exhibits no further change. 

The static pressure at the stator exit is then recalculated as 


PS = 


Ps 


fl + 2 ^Jl 1 M^ 1 


(B35) 


where Mg is the input value. It should be noted that pg does not necessarily equal 
Pjj or Pq (it depends on A^/Ag). 

Turbine power output . - The turbine power output is obtained from 

= W 4 AlVp(l - Zg) (B36) 

where w^ is set equal to Wg and Ah^ is obtained from the ADIABT subroutine. In 
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the first pass through the main program the turbine leakage flow fraction Zg = 0. On 
subsequent passes, the value of Zg that is generated in the applicable subroutine is 
used. 

Turbine thrust calculation . - Thrust is calculated at the louver exit. It is assumed 
that the louver is in the flow passage of constant annular area. It is first required to 
determine the drop in total pressure across the exit louvers. The expression used is 




Pl = Pd 


yo "I p 
1 +_S m2 

I O •L' 


rs/(ys-i) 


(B37) 


where the total pressure loss coefficient across the louvers is a program input. 
This value of is also assumed to be applicable to the exit louver of the fan. In 
equation (B37) it is assumed that yj) = = yg- This is valid because y is dependent 

on the fuel-to-air ratio (which does not vary) and the total temperature (which is con- 
stant). Then from the louver total- to static-pressure ratio Pl/Pqj the louver exit 
Mach number is calculated from the expression 


M 


L 


r 


> 




, 2 

/^l\ _ 1 

\ 

rg- 1 

V^oj 


V 

- 

J 


1/2 


(B38) 


in which the positive root (absolute value) of equation (B38) is used. The exit velocity 
is then 

’^L = “lV’'sS''s'^L 

where Tj^ is obtained from 


T 


L ~ 



(B40) 
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The turbine thrust is 


F 


T " 


w^(l - Zs)VlCt 

g 


where C,p is the thrust coefficient, a program input. 

Turbine duct areas. - The annular area at the stator exit is 


(B41) 


Ag =W4(1 - Z5) 


Vt^ I 

^p'ff/s 


(B42) 


Then when Aj^/Ag from equation (B30) is known, the flow area at duct exit is obtained 
from 





(B43) 


Fan Calculations 

Once the turbine calculations are completed yielding the turbine station thermody- 
namics, turbine thrust, and the turbine power output, the fan is then analyzed. As pre- 
viously mentioned, only fan station thermodynamics will be considered because fan aero- 
dynamic details are not handled in this program. The fan calculations required the pre- 
scribed inputs of fan-stage pressure ratio Pi 3 /Pji> average inlet Mach number Mjj, 
fan tip speed Up and fan-stage adiabatic efficiency 77 ^. The fan discharge flow is 
considered axial in this analysis. 

Fan flow rate. - The fan upstream total temperature and total pressure are set as 


> (B44) 


The pressure drop in the inlet (effect of struts, bellmouth friction, etc. ) is read in as a 
total-pressure- loss coefficient 


Tio = Tq 


Plo “ Pq 
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= 


Aj^ 


ph-p 


11 


(B45) 


The total-pressure ratio of the rotor inlet is then 


Zii 

Pio 


1 - CO 


10 


^ 1 + 0.2 


so that the total pressure at the fan inlet is 


Pil 



(B46) 


(B47) 


With fan-stage pressure ratio an input, the total pressure at the fan-stage exit is 

Pb “ 

The fan actual enthalpy rise is now obtained using the ADIABT subroutine and the inputs 
T^l, p'jj, p^g, and T?p to obtain 


Ahp = - h^j (B49) 

The fan exit total temperature Tjg is also an output of the subroutine. Neglecting 
rotor bearing friction and windage losses, the fan power is determined directly from the 
previous final value of turbine power as 



The required fan air mass flow rate is then obtained from 


w 


11 " 



(B50) 


(B51) 
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Fan state conditions. - Rotor average inlet static pressure (for y = 1. 4) is 


Pll = 


Pll 


1 + 0. 2 M 


11 


3. 5 


The rotor inlet static temperature is 


11 


Ml 
1 + 0. 2 M 


11 


and the average axial inlet velocity is 

Vjl - 49.1 |/t^j (MjP 


(B52) 


(B53) 


(B54) 


The total pressures at each station are obtained next. The duct total pressure loss 
is an input to the program as the ratio Then, the duct exit total pressure is 


Pl4 = 



(B55) 


where pj^g is given in equation (B48). The pressure drop across the louvers is 


PH 


1 - CO, 


1 - 


2 ^ 

1 + O.2MJ4I 


(B56) 


where the total pressure loss coefficient across the louvers is the same value used 
for the turbine flow stream. 

The Mach number at the duct exit is obtained from the assumption that the 

static pressure at stations 15 and 14 is equal to ambient pressure (pj^^ = p^^ = Pq). 

Then, with the total- to static-pressure ratio at the duct exit P 14 /PJ 4 - is obtained 

from the expression. 
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Mi4 = 


i\Po/ 


(y-l)/y 


~i 1/2 


y - 1 


(B57) 


in which y=1.40. Now that P 15 /PJ 4 is known, the total- to static-pressure ratio at 
the louver exit can be obtained as 


Pi 5 W 4 / \Pl4/ 


(B58) 


Then, from equation (B57) the expression for louver exit Mach number (for y = 1. 40) is 


M 


15 


5.0 




0 . 286 


,1/2 


15 ^ 
.P0 7 


- 5.0 


(B59) 


The louver total temperature is obtained by assuming adiabatic flow through the duct and 
louvers; that is. 


»T» — T*» — 'T'» 

^15 ■ ^14 - ^13 


(B60) 


where TJ^g is obtained from the output of the ADIABT subroutine (eq. (B49)). The 
static temperature at the louver exit is then 


Ti5=7 — 

(1 + 0.2MJ5) 

and the discharge velocity is 

= (B62) 

For the fan stage, axial discharge is assumed, so that the louver exit velocity in equa- 
tion (B62) is the axial discharge velocity. 

The stator exit axial Mach number M^g is taken to be some function of the differ- 
ence between Mj 4 and Mj^g in the form 
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Mi3 = Mjj+Kj,(Mj5-Mj,) 


(B63) 


where Kp is an input. The stator exit static conditions are then 


Pl3 


Pi 3 


2 ^ 
1+0.2 Mjgj 


(B64) 


and 


Ti3 = 


T’ 

^13 


1 + 0. 2 M 


13 


(B65) 


The average exit axial velocity is 


Vi3=49.1Mj3^ 


(B66) 


The static temperature at the duct exit is 


T* 

^14 


14 


1 + 0. 2 M 


14 


(B67) 


And the duct exit axial velocity is calculated from the expression 

Vi4 = 49. 1 14 


(B68) 


Thrust . - The fan thrust is given by the expression (since Pjg = Pq) 




^ll^lS^T 


(B69) 


where the fan thrust coefficient is an input taken to be equal to the turbine thrust 
coefficient. 
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The tip-turbine fan total thrust is obtained as the sum of the fan and turbine thrusts 
(eqs. (B41) and (B69)) as 


F = + Fj, (B70) 

This value of total thrust, which is the result of the initial choice of turbine mass flow 
rate Wj, will generally not be equal to the value of thrust desired. 


Fan and Turbine Geometry 


After the calculation of the initial total thrust using the first value of turbine mass 
flow rate, fan and turbine areas and geometry are calculated. For the fan stage, at the 
rotor inlet (station 11), the flow area is 


A 


11 


wRT \ 

. pv 4i 


The stator exit area is (station 13) 

13 " 

Pi 3^1 3 

and the duct exit area is (station 14) 

■^14 

P14V14 


and the louver exit (station 1 5) 


A 


15 


^ 11^^15 

P15V15 


(B71) 


(B72) 


(B73) 


(B74) 


The fan rotor assembly is shown in figure 11. The fan rotor tip diameter is calculated 
from the expression 


dp 1.128 


(B75) 



where the fan hub-tip ratio at the rotor inlet (Yp = p) is a program input. The 

fan hub diameter is then 


^F, h " ^F, t^F 


and the fan rotor blade height 


^F “ ^ ^F, t^^ ■ ^F^ 


(B76) 


(B77) 


Between the fan rotor blade and the turbine blades is a shroud which acts as the 
carrier for the turbine blades. The carrier must be sized to provide an adequate struc- 
ture to withstand inherent thermal stresses and allow for fabrication and possible cool- 
ing requirements. In general, the ring lei^th increases somewhat with increasing fan- 
tip diameter. The value of rir^ height Ar^ will be an input. Then the turbine hub 
diameter is 


“^T, h “ ^F, t + ^ "^^c 


The turbine tip diameter for both full and partial admission is 


^T,t 


h2 ^ ^5 /360 


1/2 


(B78) 


(B79) 


where the turbine rotor annular area is obtained by first setting Ag = Ag. (On 
subsequent iterations, using the value calculated in the turbine subroutines, Ag equals 
the area at the turbine rotor inlet normal to the axial component of velocity. ) The arc 
of admission, 9 , in equation (B79) is in degrees and is an input. 

The turbine hub-to-tip diameter ratio is then defined as 


Y 


T 



(B80) 
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and the turbine rotor blade height is 


% t “ ^T, h^ 

At this point in the main program, the appropriate turbine subroutine (appendix D 
or E) is called to obtain an accurate value of turbine efficiency. The subroutine is 
called at this point because it requires the turbine tip diameter to obtain the turbine 
mean line rotational speed Urp; that is 


Uf t 

Ut = (cLp h + ^ (B82) 

The turbine subroutine goes into the detail of internal turbine analysis including blade 
loss and leakage flow considerations. After obtaining the new value of efficiency, the 
calculation returns to subroutine ADIABT to obtain a new value of Ah and total tem- 
perature Tg. The calculations then commence with equation (B36). This procedure is 
repeated until 



0.001 


This is done for each input choice of exit stator exit Mach number Mg and each value of 
inlet duct mass flow Wj used. 


Performance Parameters 
Corrected fan mass flow is obtained as 


11 


. 1/2 


^eq = ^ll 


,518.7. 


Pjl 

14.7 


(B83) 
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The corrected total thrust is 



F 


(B84) 


F 

6 Ph 

14.7 


The burner fuel flow is 


w. 


f,b 



and the burner thrust specific fuel consumption is 


SFC 



(B85) 


(B86) 


where total thrust F is from equation (B70). 

The fan bypass ratio is defined as the ratio of fan stage flow rate to turbine flow 
rate or 


W-. 

BPR = -ii. (B87) 

W4 

The fan thrust augmentation ratio is defined as the ratio of lift fan thrust to ideal thrust 
obtainable from the expansion of the supply gas to ambient pressure. This expression is 

AUG = — - — (B88) 

gg 


where 


F 


id,gg 


WjV 


id;g g 


g 


and the ideal expanded discharge velocity is 


(B89) 


(B90) 
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The ideal enthalpy drop 


Pi, 



Viy, and t] 


Ahj^ is obtained using subroutine ADIABT with the inputs 
where 77 is set equal to one. 


T' 

1 » 


Desired Total Thrust Determination 

On completion of all program calculations for an initial choice of turbine mass flow 
rate which yields a value of tip-turbine fan total thrust, the prc^ram repeats the calcu- 
lations using another choice of inlet duct mass flow Wj. This procedure yields two 
values of total thrust corresponding to the two choices of inlet duct mass flow. The pro- 
gram then uses a linear interpolation to obtain the correct value of Wj corresponding to 
the desired total thrust. The program is then rerun using the new value of Wj^ to obtain 
the turbine and fan designs. This procedure has been shown to yield less than 1/2 per- 
cent discrepancy in total thrust. 
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APPENDIX C 


TURBINE LOSS ANALYSIS 

The determination of turbine stage efficiency and turbine state conditions (appen- 
dixes D and E) in parametric studies that involve a relatively large variation in overall 
pressure ratio, inlet temperature, stage velocity diagram parameters, and turbine 
geometry requires a realistic estimation of the loss factors involved (rotors and sta- 
tors). In addition, the magnitude of the overall efficiency of the turbine will help deter- 
mine whether one or two stages is desirable for the application. 

A typical set of velocity diagrams for a single-stage impulse turbine are given in 
f^ure 12 (station designations correspond to those used in the program). Although the 
flow entering the inlet stator from the scroll may not be axial and may vary in approach 
direction around the periphery, for simplicity, the velocity at the inlet stator will be 
assumed to be axial. The error in inlet- stator- loss coefficient that this assumption 
causes is assumed to be negligible. 

If a two-stage turbine is used, the interstage stator will generally have a greater 
amount of turning than the inlet stator because the gas flow into the interstage stator is 
not axial. The velocity diagram of the inlet and exit absolute gas velocities for an inter- 
stage stator is given in figure 13. Thus, the turning angle of the interstage stator is 
AQ! = arj - Q!g. 

The loss analysis is divided into three parts: (1) the stator, which has a sizable 
static pressure drop (2) the impulse rotor blades, and (3) impulse exit stators. Diffus- 
ing exit stators are discussed in appendix D. 


Inlet and Interstage Stator Efficiency 


The efficiency of the turbine inlet stator is defined as 


^N 


Energy loss 
Ideal energy 


(Cl) 


where the ideal energy is obtained from the enthalpy change Ahj,^ in the isentropic 
expansion from inlet total pressure to outlet static pressure. The ratio of energy loss 
to ideal energy is referred to as the energy- loss coefficient or 

_ / Energy loss 

s ja ~ 1 

\Ideal energy 
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The loss coefficient of a row of stator blades is recognized, in general, to be the sum- 
mation of three principal loss components: (1) the profile loss (2) secondary- 

flow losses (?se)N> (3) shock losses (^gh)N* Thus, for a tip-turbine stator 

“ ^^p^N ^'^se^N ^^sh^N 

The profile loss (loss due to basic blade element) is depenent on a large number of fac- 
tors, such as, blade shape, turning angle, blade maximum and trailing edge thicknesses, 
solidity, incidence angle, and Reynolds number. For the subject stator, it will be as- 
sumed that the following conditions are met in the stator design: 

(1) Optimum solidity (e.g., Zweifel loading coefficient) 

(2) Optimum incidence angle (~0°) 

(3) Optimum blade profile shape 

For these conditions, the profile loss will depend primarily on the air turning angle and 
Reynolds number: 


(^p)N = (C4) 

The variation of profile- loss coefficient with air turning angle has been determined 
experimentally for typical stator blades, as shown in figure 14 (from refs. 12 and 13). 

A wide spread of loss values is indicated. For the range of calculations to be covered 
for varying pressure ratio across the turbine, the inlet stator outlet angle is not 
expected to vary by more than around 15*^ (e.g. , from around 60° to 75°). Equation (C4) 
then becomes 


In this case for simplicity, a single value of loss coefficient can be adopted to be repre- 
sentative of the range of o!g and thus eliminate the need for an iteration on The 

reference value of (?p)jsj g used in equation (C5) for the inlet stator is equal to 0. 030. 

For the two-stage turbine, the interstage stator will have a greater turning angle 
depending on the value of absolute rotor exit angle org (fig. 13). The range of turning 
angle for this case should also be relatively small (from around 75° to 85°). Once again 
a single representative (^ *)^ „ value of 0. 040 was chosen for the range of turning angle 
considered. Here again, the error associated with this assumption is assumed small. 
Thus, in both the inlet and interstage stators (?p)j^ is only a function of Reynolds 
number. 

Inasmuch as the calculation will include variations in turbine size, provision should 
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be included for reflecting an effect of Reynolds number on the profile loss. This will 
be done in the form 


(^p)N 

where the reference value of Reynolds number Re* used in the program is 1. 5x10 . 
For this analysis, if a relatively constant value of blade aspect ratio is considered as 
turbine size is varied, then Reynolds number can be defined in terms of the blade pas- 
sage height and expressed as in reference 14 as 


Re = 


w 


m 


(C7) 


where w is the turbine flow rate, r,p is the mean radius, and ji is viscosity. For 
simplicity, the viscosity can be expressed as a function of the total temperature at the 
stator inlet. Values of viscosity are obtained using the VISCON subroutine (appendix G). 

Experimental data for stators indicates that the loss coefficient for a given blade 
row stays essentially constant with outlet Mach number up to sonic velocity. However, 
as shock formation increases with supersonic Mach numbers, the loss will increase. 
Thus, a shock-loss coefficient can be established as a function of outlet Mach number, 
such that for 1.0 




(C8) 


The shock-loss coefficient given by equation (C8) for values of Mach number greater 
than 1. 0 can be established from experimental data. Such data are currently unavailable. 
Because of the difficulty in obtaining valid experimental data, a preliminary estimate 
was established from a plot of stator efficiency against stator exit Mach number shown 
in figure 15. The curve presented is based on limited analytical estimates available 
from unpublished Lewis calculations for stator designs with high exit Mach numbers. 

A representative shock-loss coefficient variation was obtained from this plot from 
the relation ^ = 1 - rj such that 


^ sh 

The corresponding 
was obtained as 


supersonic ~ ^^N^sonic ^^N^sonic ~ ^^N^supersonic 


plot of is shown in figure 16 for which an analytical variation 
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|^^. 0.0075(m2-m)„ 


(CIO) 


where is inlet stator outlet Mach number (M > 1). 

The losses due to secondary flows and casing friction are fairly complex and diffi- 
cult to express as a general functional relation. However, for simplicity and for the 
purposes of this analysis, the secondary- loss coefficient was accounted for by taking it 
to be proportional to the profile- loss coefficient. Thus, 

se^N ^ ^^se^N^^ p^N 

But, since secondary losses are not sensitive to blade Reynolds number effects, equa- 
tion (Cll) becomes 


se^N “ ^*^seM^ 

Values of (Kgg)jg> which is a function of turning, are indicated to vary from around 0. 5 
to somewhat over 1. 0 in the literature. The value selected for (Kgg)jvf was 0. 75. 

Then the total loss coefficient, from equations (C3), (C6), (CIO), and (C12), becomes 




' Re 
kRe* 


- 0.2 


<K^se)N 


+ 


sh^N 


(Cl 3) 


Rotor Blade Efficiency 

The rotor blade is designed for impulse flow and a shrouded tip. The efficiency of 
the turbine blades, which are treated as a cascade, is defined in terms of a kinetic en- 
ergy loss by the expression 


where the energy-loss coefficient for a row of rotor blades is defined as the ratio of 
the loss in kinetic energy to the ideal outlet kinetic energy. As in the case of the stator, 
the energy loss comprises the profile loss, secondary flow losses, and shock losses. 

The loss coefficient is then 


^ B ^ p^B ^ se^B ^ sh^B 


(Cl 5) 
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The variation of profile loss coefficient with air turning angle has been determined 
experimentally for rotor blades, as shown in figure 17 (from refs. 12 and 13). A wide 
range of values has been obtained. Only the data of Ainley are definitely known to be for 
impulse blades. Thus, the selected representative relation for )g against A)3 will 
be influenced primarily by these data. Figure 18 shows the selec^d curve used for 
turning angles up to 140°. In equation form, for Re*, 




0. 025 + 0. 035 



(C16) 


0 

where (5p)g corresponds to Re* = 1. 5x10 with Re defined as before in equation (C7). 
The turnii^ angle A/3 for simplicity and in order to avoid an iteration, is assumed to 
be equal to twice the inlet angle; that is. 


A/3 = 2/3 


inlet 


(C17) 


The change in the value of profile-loss coefficient due to this assumption is expected to 
be small compared with the spread in the data variation as indicated in figure 17. 

The influence of Reynolds number will be taken to be the same as for the inlet stator 
loss or 






2 


(Cl 8) 


The losses due to secondary flows are treated in a manner similar to the inlet stator 
procedure. The secondary- loss coefficient is taken to be proportional to the reference 
profile-loss coefficient. Thus, 


8se)B = (Kse>B(«?B (C19) 

Values of (Kgg)g that are a function, of turning, are indicated to vary from around 0. 30 
to 0. 50 (ref. 12). A (Kgg)g value of 0. 40 was used in the program. 

The shock- loss contribution can be established from an estimated variation of blade 
efficiency with rotor relative inlet Mach number. Since no data were available for this 
variation, an assumed shock- loss variation three times the magnitude of the stator shock 
loss was adopted as shown in figure 19. The resulting expression takes the form 
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where the constant is equal to 0.0225. Thus, the total loss coefficient, from equa- 

tions (Cl 5) and (C18) to (C20), is 


>B “ 


Re \ 
Re*j 


- 0.2 


+ (K J 


se'B 


+ (l 


sh^B 


(C21) 


Exit Stator Total- Pressure- Loss Coefficient 

The loss coefficient for the impulse exit stator is obtained from the approach used 
for the rotor blades. The profile loss is composed of mostly surface friction and dif- 
fusion. In terms of a kinetic energy loss, the profile-loss coefficient is 

(^p)g = f(A«,Re) (C22) 

The relation of (4 *)g 3,nd Aof is taken to be the same as for the rotor blade (eq. (C16)) 

2 5 

(f *)q = 0.025 + 0.035/— V (C23) 

\ 100 / 

Since the primary purpose of the exit stator is to straighten the flow to axial, the amount 
of turning is just equal to the exit stator inlet angle Og or a^; that is, 

AQ'=|ffgl (C24) 

is used with equation (C23) for the exit stator. Since Mach numbers in the stator are al- 
ways less than one, the shock loss term is not needed. Thus, with Reynolds number and 
secondary losses. 


where (K^g)s = (Kgg)g = 0. 40 and Re* = 1. 5x10®. 

The exit stator total-kinetic -energy- loss coefficient is converted to an exit 
stator total-pressure-loss coefficient Kg for ease of calculation. Conversion of the 
kinetic- energy- loss coefficient to a total-pressure-loss coefficient is developed in ref- 
erence 12 (section 2.7). This expression in general form is given as 


Re 

Re* 


- 0.2 


+ (^se)s 


(C25) 


37 




1 + yM 


(C26) 


where Mg is exit stator outlet Mach number and the total- pressure- loss coefficient Kg 
in the preceding expression is defined as Ap’/(p' - P)g- Then, the final expression for 
total-pressure-loss coefficient is obtained from equations (C25) and (C26) as 




Re 


a. 5X10 


- 0.2 




1 + yM*^ 


(C27) 


Since 


(P’ - P)g = (P’ - P)e or 8 


then 


Mg = M 


6 or 8 


for an impulse stator. Thus, the total-pressure loss factor was applied to the stator in- 
let, so that, for the exit stator. 


Ap^ - Kg(p' - p)g 3 


(C28) 
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APPENDIX D 


SINGLE-STAGE TURBINE SUBROUTINES 


This subroutine calculates the value of turbine total-to-total efficiency that is re- 
quired by the main program. The single-stage turbine analysis presented investigates 
three types of turbines: (1) with impulse exit stators, (2) without exit stators, and 
(3) with diffusing exit stators. These three types of turbines are described in sub- 
routines, TURB-1, TURB-2, and TURB-4, respectively. The turbine design- point effi- 
ciency includes the effects of blade-row losses, partial admission and its associated 
losses, and turbine seal leakage which is included as a reduction in turbine mass flow. 

It is assumed that the velocity diagrams are not affected by either partial admission or 
leakage flow losses. 

The turbine calculations are for the mean line since the turbine hub-to-tip radius 
ratio is of the order 0. 90 or greater. The turbine mean line velocity diagrams are de- 
termined by the choice of exit stator exit Mach number, wheel speed, inlet stator spout- 
ing velocity, and the designation of turbine type. The turbine station designations used 
in the analysis for the single-stage turbines are given in figure 1. 

The inputs required from the main program for all three cases are 


Turbine mean line radius 

Turbine mean line velocity 

Total temperature of the gas into the inlet stator 
Total pressure of the gas into the inlet stator . . 
Total mass flow of gas into the inlet stator. . . . 

Exit stator exit axial Mach number 

Exit stator exit static pressure 

Fan tip diameter 

Fuel-to-air ratio at the burner exit 



Mg 

Ps 



In addition to these inputs, several geometry and loss-coefficient inputs are required for 
the leakage and blade-row efficiency determinations, respectively. 


Subroutine TURB-1 - Impulse Exit Stators 

Assumptions . - The velocity diagram for this stage is shown in figure 10. The 
diagram is determined by the choice of inputs. The design assumptions for the impulse 
exit stator case are constant static pressure through the rotor and exit stator 
(Ps = Pg = Pg) constant annular area across the rotor. In addition, the rotor inlet 
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axial velocity is initially prescribed by the condition c = and a trial value of in- 

X. j 0 k) 

let stator efficiency ri-^ is prescribed. 

Rotor inlet conditions . - The inlet stator exit velocity is obtained from the actual 
total- to static -enthalpy drop across the inlet stator, Ahj^ = h^ - hg as illustrated in 
the enthalpy-entropy diagram of figure 20. The value of Ahj^ is determined in the 
ADIA.BT subroutine described in appendix G. The specific inputs required for the deter- 

f 

mination of Ahj^ are T^, p^, — , p^, and 77 ^. The ADIABT subroutine also de- 

a 4 

termines the value of the static temperature at the inlet stator exit Tg. 

The velocity at the inlet stator exit is then calculated from the expression 


V5-|/2gJ Ahj^ (Dl) 

Values of and at the inlet stator exit are obtained from the subroutines of ap- 
pendix G. Inputs to the gas property subroutines are temperature and fuel-to-air ratio. 
Thus the value for sonic velocity is 


The inlet stator exit Mach number is then calculated from 


(D2) 


M 


5 ~ 



(D3) 


Once Mg is calculated, the new value of inlet stator efficiency is obtained from 




(D4) 


where the loss coefficient is obtained from equation (C13) of appendix C. With the 
value of T 7 ^ known, the subroutine iterates on Mg until 




0.002 


(D5) 


The rotor inlet velocity diagram is then established in the following manner; The in- 
let stator exit angle cHg (measured from the axis) is calculated as 
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-1/M 

Q!p. = COS —1 = cos M 

^ \v 


(D6) 


where g is the initial prescribed value (=Mg). The rotor inlet axial velocity is then 


Vx, 5 = V 5 cos 


p7) 


and the tangential velocity component is 


Vu, 5 = ^5 sin 


(D8) 


Then, when the rotor mean line velocity is known, the tangential component of the 
gas may be determined from 


Wu,5 = V,,5-Ut 


(D9) 


The relative velocity at the rotor inlet is 




(DIO) 


The rotor relative inlet blade angle measured from the axis is 


= tan 


l/'^u 


(DU) 


The total pressure at the rotor inlet is obtained from pg and Mg as 


P5 = P5 


1 + - (r - i)M^ 

2 


75/(75-1) 

5 


(D12) 


and the density is 


P5 = 


144 p \ 
RT /g 


(D13) 


41 



Rotor blade efficiency and exit conditions . - The blade efficiency Tjg is obtained 
using the loss coefficient approach developed in appendix C. The blade efficiency is 

77b = 1-|b (D14) 

where the rotor blade loss coefficient is calculated from equation (C21) of appen- 
dix C. The relative gas velocity at the rotor exit is then, by definition 


We=W5Vr?B 


(D15) 


The energy loss in the rotor blade is calculated from the expression 


E 


B 


^5 6 

2gJ 


or in terms of rotor blade efficiency 


E 


B 



(1 V-q) 


(D16) 


(D17) 


From this value of the energy loss, the rotor exit fluid properties are calculated by as- 
suming that the gas specific heat Cp g is obtained using rotor inlet conditions as the in- 
puts to subroutine SPHTCP (appendix’ G). This assumption is valid because a relatively 
small variation in static temperature exists across the impulse rotor. Thus, 


E 


^6 = ^5 + 


B 


"p,6 


The static density is, since pg = pg. 


^6 “ 


144 p I 
TR /g 


(D18) 


(D19) 


For constant annular area across the rotor, the rotor exit axial velocity from con- 
tinuity is 
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(D20) 



The relative tangential component of the outlet flow is 




(D21) 


and the component g is 


Vu,6=W„,6 


Un 


Then the relative flow angle at the exit is 


^6 



and the absolute angle is 


“6 



The absolute flow exit velocity is then 



and the absolute exit Mach number Mg is obtained as 


Mg 



(D22) 


(D23) 


(D24) 


(D25) 


(D26) 


The relative exit Mach number is 




The total temperature at the rotor exit is then 


(D27) 


T6 = T6 


1 +-(r - i)M^ 


(D28) 


and the total pressure is 


Pe ^ ^6 


l +-(y - 1)M^ 
2 


yg/Crg-i) 

6 


(D29) 


Exit stator . - The exit stator is specified as impulse, that is, constant static pres- 
sure. The total pressure at the stator exit is determined from 

Ps = Pg - Kg(p’ - p)g (D30) 


where the exit stator total-pressure-loss coefficient is calculated as described in 
appendix C. 

The exit stator exit axial Mach number is calculated by noting that pg = pg so that 



(D31) 


The value of Mg calculated is then compared with the value of Mg specified in the 

main program. K the values are different, the starter value of M which is initially 

0 

obtained as M ^ = M^,, is adjusted according to 

X, D O 


M 


X, 5 


1 = M 

n+1 X, 5 


n 


^S(pr escribed) 

_ ^S(calculated)_ 


(D32) 
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The entire calculation is repeated until convergence is obtained according to 



0.001 


(D33) 


After the final stator exit Mach number is determined, the stator exit fluid proper- 
ties are calculated. The static temperature is, since Tg = Tg, 

T* 

Tg = ^ (D3^) 

1 +i(y - 1)M^ 

L 2 Jg 

and the stator exit velocity, which is prescribed as axial, is then calculated as 


Vg = (u-^gyRT^^ (D35) 

The total pressure at the stator exit is calculated knowing the value of pg and the input 
value of Mo from the expression 


Ps = Ps 


1 + i(y - 1)M^ 

2 


yg/Crs-i) 

s 


(D36) 


Turbine-stage efficiency . - The turbine- stage efficiency required for the main pro- 
gram is a total-to-total efficiency since the main program only deals with total condi- 
tions. This value of turbine efficiency will replace the assumed value used as an initial 
value by the main program as discussed in appendix B. The stage efficiency is calcula- 
ted initially for the case of no leakage (i. e. , Wg = w^). This calculation yields the maxi- 
mum attainable efficiency for the previously calculated velocity diagrams. By definition, 
the turbine- stage efficiency is 


77rp — 



(B37) 


where, for a constant turbine pitch diameter, the no- leakage actual power is 
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_ Ut(Wu , - Wu 6)W4 


(D38) 


P 


T 


gJ 


and the ideal power based on total-to-total conditions is 

^T, id " "^4^^^T^id 

As illustrated in the h-S diagram of figure 20 


(D39) 


(D40) 


The ideal total-to-total enthalpy drop from turbine inlet to exit stator exit is calculated 
using the ADIABT subroutine. Inputs required are p^, T^, and p^ along with 

= 1. 0. The subroutine also issues as output the isentropic total temperature at the 
exit stator exit (Tg)-^. The efficiency is then obtained using equations (D37) to (D40) as 


Ut<Wu, 5 


W„.6> 


gJ(Ah^)j^ 


(D41) 


This value of efficiency is then adjusted to account for the turbine flow leakage at the 
inlet stator exit to the fan passage. This leakage reduces shaft power output by reducing 
the mass flow across the rotor. The adjusted efficiency is then given as 

^T, leak "" ’^T^^ ' ^5^ 

The value of Z^, which is defined as the ratio of the turbine seal leakage flow rate to 
turbine inlet flow rate, is given as 


^leak 

Z (D43) 

W4 

Details for computing the leakage fraction are given in appendix F. 

The total-to-total efficiency as given by equation (D42) corresponds to the definition 
of efficiency used in the main program (appendix B). 

Partial admission effects on design-point efficienc y. - Additional losses in power 
will arise in partial admission due to (1) windage losses (centrifuging of gas in the non- 
active part of the wheel), (2) filling and emptying losses (expansion and loss of momen- 
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turn of the gas in the blade passages at the end of the sector of admission), and (3) mixing 
losses (mixing of the working gas with stagnant gas in the blade passages at the ends of 
the admission arc). The method developed by Stenning (ref. 16) is used to estimate these 
losses. 

The filling and emptying loss affects the momentum of the gas leaving the rotor 
blades so that, in effect, the relative velocity leaving the blade is reduced. Thus equa- 
tion (D15) becomes 


We 



where 



In equation (D45), a is the arc length of admission given by 


a = Trd 


T, m 



(D44) 


(D45) 


(D46) 


and s is the blade spacing at the mean line, which is obtained from the blade solidity 
o = c/s and aspect ratio AR = Hrp/c to give 



ct(AR) 


(D47) 


The values of solidity and aspect ratio are program inputs. 

The effect on the actual power due to pumping and mixing of the fluid in the rotor is 
given as a modification of equation (D38). 


U 

= — (W, 


w 


gj 


u,5 - '^u,6^'"4 


)w.- 


1.4 u;^c W4 
2gJ V^a cos Q!g 


Simplification results in 


U w 

p _ t"^4 

T 

gJ 


1.4U^c 

m _ w ^ 

u, 5 u, 6 2 V -a cos O'. 


(D48) 


(D49) 


47 


The revised turbine total-to-total efficiency expression replacing equation (D41) is, then, 


— 


Un 


(Ah^)idgJ 



1.4 u;^c 
2Vga cos Q!g 


CD50) 


Turbine annulus flow areas . -- After completing the velocity diagram and efficiency 
calculations, turbine annulus areas were obtained. The inlet stator exit annulus area for 
full admission is giverrby the expression 


A 


5 


144 w^ 
^x,5^5 


(D51) 


Rotor inlet annulus area is for constant axial velocity assuming no change in static con- 
ditions 


^ 5 ' " ~ ^ 5 ^ 


and, for prescribed constant annulus area. 


^6 



The exit stator exit annulus area is 


RgTsW4(l - Zg) 


(D52) 


(D53) 


(D54) 


Turbine design parameters . - The speed-work parameter (as defined in ref. 14) is 
given as 


X = 




AV 


u 


(D55) 


But, for constant turbine pitch diameter 


AV^ = AW^ 
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thus, 


A = 


Un 


Wu,5-W„,6 


(D56) 


Another parameter often used in describing turbine efficiency is the blade-jet speed 
ratio, defined as 


V = 



(D57) 


where the velocity term Vj is the ideal jet velocity corresponding to the ideal expansion 
from inlet total to exit static conditions across the turbine 



(D58) 


Subroutine TURB-2 - No Exit Stators 

In this design the whirl component of the exit velocity is lost because it is not avail- 
able to produce turbine thrust. This subroutine differs from the TURB-1 subroutine in 
that Mg = g, Pg = Pg, Tg = Tg, and Vg = g. Specifying a value of Mg in the 
main program sets both the rotor inlet and exit axial velocity for this case. Definitions 
of stage work and turbine total efficiency are not altered. The calculation procedure and 
sequences for this case are the same as those of the impulse stator case of TURB-1. 


Subroutine TURB-4 - Diffusing Exit Stator 

In order to get the maximum power from impulse turbines along with obtaining low 
exit velocities, the use of diffuser exit stators is required. This type of vane will be 
designed to, not only turn the flow, but also to permit a reduction in velocity such that 
Vs < Vg. The velocity diagrams for this case are the same as those shown in figure 12. 
For this type of stage and for an ambient exit pressure (pg = Pg), the rotor inlet and exit 
static pressures will be subambient. 

The amount of diffusion in a blade row will depend on the loss that can be tolerated. 
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The losses in the blade were expressed in terms of the total- pressure- loss coefficient: 


1-p^ 




P6 


^ -Ps 
P6 


1 - 


Pe 


P6 / , 1 

1 - 1 


(D59) 


The value of the total-pressure-loss coefficient was estimated from the relation 

.3. 7D 




2agK| e^ 


cos ]3, 


(D60) 


where D is the diffusion factor, is the stator solidity (taken to be 1.2), and KJ is 

o 

equal to (Kg cos ^0/2cTg)£j_Q (taken to be 0. 005). The diffusion factor for turning back to 
axial is defined as 


D = 



(D61) 


The value of exit velocity Vg can be calculated from equation (D61) knowing rotor exit 
flow conditions (Vg), an input value of exit stator solidity, and an input value of allow- 
able diffusion factor. In order to insure that a solution is satisfactory, the program iter- 
ates on diffusion factor D until the calculated value of exit stator exit Mach number 
agrees with the prescribed value to within 0. 001. This procedure then yields a specific 
value of diffusion factor (program output) for each choice of exit stator exit Mach num- 
ber. The value of diffusion factor calculated for a choice of Mg may be beyond the 
practical limit thus limiting the range of exit stator exit Mach number. 
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APPENDIX E 


TWO-STAGE TURBINE SUBROUTINE 

For turbine designs with large pressure ratios, it may be desirable to design a two- 
stage turbine in order to reduce the pressure ratio per stage and increase stage and 
overall efficiencies, thus reducing required turbine weight flow. This section provides a 
calculation of turbine total-to- total efficiency for use in the main program. The two- 
stage design is based on full admission, impulse rotors, and near-equal stage work with 
each individual stage made up of a stator and rotor. In addition, there is an impulse exit 
stator to turn the gas flow to the axial direction to increase turbine thrust. 

The turbine design point efficiency includes the effects of blade row losses. In addi- 
tion, it is adjusted to include the reduction in turbine mass flow due to turbine-to-fan 
seal leakage. It is assumed that the velocity diagrams are not affected by this reduction 
in flow rate. Only the leakage path upstream of the first-stage rotor is considered in the 
analysis. 

The two- stage turbine is mounted at the tip of the fan rotor blade as shown in fig- 
ure 3. The turbine station designations for the two-stage turbine design are given in fig- 
ure 9, and the typical velocity digrams are shown in figure 13. Provision is made in 
the analysis for a variation of mean diameter between first- and second-stage rotors. 
Most inputs required from the main program are identical to those required for the 
single-stage designs given in appendix D. The rotor inlet axial Mach number for each 
stage is at first set equal to the exit stator exit Mach number, that is, M^^ g = M^ rj =Mg, 


First Stage 

The static pressure at the first-stage rotor inlet was determined by the expression 

P 5 = X Vp^Pq (El) 

where X is an input that enabled the program to obtain variable work split between the 
stages. For example, a value of X equal to 0.90 results in a 55-percent - 45-percent 

work split for the first and second stages, respectively. With the calculated value of p^ 

f ^ 

and the inlet conditions T^, p^, — , and an initial trial input value of 77 ^^^ g, the 

a 4 

subroutine ADIABT is used to obtain inlet stator outlet condition Tg and the enthalpy 
drop 
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5 ' ■ ^5 


(E2) 


for the expansion from and to Pg. The inlet stator exit velocity is then 


Vj = V^SJ 5 

The values of the gas constants Rg and are obtained from the subroutine GASCON 
and GAMMA, respectively (appendix G). The Mach number at the inlet stator exit is then 


M 


5 " 



(E4) 


A new value of inlet stator efficiency is then obtained from the relation 


5 " ^ “ ^N, 5 

and the inlet stator loss relations of appendix C. The value of g is obtained from 
equation (Cl3) of appendix C. Once a value of 5 has been obtained, the program 
iterates on Mg until 


M 


5ln 


- M 


5ln+l 


< 0. 002 


(E 6 ) 


The velocity diagram calculations are conducted by initially setting 

“x. 5 = “s 

Then 


1 \ 


Og = COS"^( ) 

(E8) 

^ \M / 


5 


Vx, 5 = ^5 cos og 

(E9) 

= "5 

(ElO) 
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The relative tangential velocity component is 


Wu,5 = V5-Ut,6 


and 



Then 


W< 






and 


M 


r, 5 


W 

^rgRT^ 


The total pressure at the rotor inlet is 


P5 = PS 


1 + i (y - 1)M^ 
2 


_n 75/(75-1) 


j 5 


and the density is 



After completing the calculation for the rotor inlet, the exit conditions are 
by calculating the losses across the rotor blades. The relative exit velocity is 


W6=Wg 



where the rotor blade efficiency ?7g g is calculated from 


(Ell) 


(El 2) 


(El 3) 


(El 4) 


(El 5) 


(E16) 

initiated 

(E17) 
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(El 8) 


The loss coefficient is obtained from equation (C21) of appendix C. The values of g 
and Wg permit the calculation of rotor energy loss: 


E 


B,6 


W? - 
5 6 

2gJ 


(El 9) 


Next. C is obtained from the SPHTCP subroutine of appendix G using rotor inlet con- 
p, b 

ditions, Pg and Tg. Then 


For an impulse rotor, 


^E 


^6 = ^5 + 


B 


P/6 


(E20) 


P6-P5 


(E21) 


and the density is 


p WliiP' 

® V RT y 


(E22) 


For the condition that Ag = Ag and for continuity 


^x,6 = Vx,5:^ 


(E23) 


The tangential component of the flow is (see velocity diagram fig. 13) 




W.. . -l/W^ - vl 6 


(E24) 


and the component g is 
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Vu, 6 = <"u - »T>6 

Then the relative gas angle at the rotor exit is 


^6 



and the absolute angle is 


"6 



Thus, 


Vg = (V^ cos o)g 

The rotor exit absolute Mach number is calculated from the expression 


(E25) 


(E26) 


(E27) 


(E 28 ) 


Mg = 



The relative Mach number is then 



The total temperature at the rotor exit is 






(E29) 


(E30) 


(E31) 


and the total pressure is 
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P6 = P6 


1 + 1 (r - i)M^ 


76 /( 76 - 1 ) 


(E32) 


These fluid properties now constitute the inputs to the second-stage stator. 


Second Stage 

The velocity diagram for the second-stage stator and rotor are shown in figure 13. 
For the stage, 


P7 - Pg - Pg 


(E33) 


The enthalpy drop across the second-stage stator was first calculated with a trial 
starter value of stator efficiency rj and inputs to the ADIABT subroutine (appendix G) 
T^, P^, - , 77 ^, and p,^. This yield’ed T,^ and 
^ 4 


7 = *^6 ■ N 

Then the stator exit velocity is 

Vrj = y^^TZh^ (E35) 


The Mach number is then calculated as 



(E36) 


A new value of stator efficiency ^ obtained from the relation 

%, 7"1"^N, 7 

The energy-loss coefficient rj is calculated in the same manner as the loss coeffi- 
cient of equation (E5). 
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The relation used for adjustir^ the mean radius of the second- stage rotor is 


where (TRM) is an input which is generally less than 1. 0 because of the fan duct contrac- 
tion. The pitch line velocity is then 

g = (TRM)U^^ g (E39) 


Once the second-stage stator efficiency g has been calculated, the program iterates 
on stator exit Mach number until 


M 


7|n 


- M 


7 n+T 


< 0,002 


(E40) 


At the second-stage rotor inlet, the axial flow Mach number M ,, is taken to be 

X, / 

equal to the axial Mach number into the first-stage rotor M„ c- The second-stator 

Xj D 

absolute exit angle a,j is then calculated as 


arj = cos 


V M 


(E41) 


And 


Vx 7 = (V cos a)ri 
7 = (V sin a)rj 


(E42) 

(E43) 


P^ = P7 


1 +l(y - 1)M^ 


77/(77-1) 


(E44) 


J7 


In the preceding equation p,^ is obtained from equation (E33), The flow density is 


_ /144 p 
\ RT 


(E45) 
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The remaining velocities and angles of the velocity diagram at rotor inlet are 


Wu,7 = V7-%,8 


^7 







For the rotor outlet, the relative velocity is obtained from 


(E46) 


(E47) 


(E48) 


(E49) 


Wg 



(E50) 


In this equation, the second-stage rotor blade efficiency is obtained from the expression 


8 = ^ ■ ^B, 8 


(E51) 


for which the loss coefficient is obtained in a manner similar to the first- stage rotor- 
blade loss coefficient. The value of second-stage rotor blade efficiency from equa- 
tion (E51) enables the determination of energy loss through the rotor based on relative 
velocities 


E 


B, 8 


^8 

2gJ 


(E52) 


then C o is obtained from subroutine SPHTCP of appendix G using T«. Then, 
p, 8 l 

^E, 

Tg = T^ +1 



(E53) 
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By specification, the static pressure is 


Ps - P7 


and the flow density at the rotor exit is 

I 144 p\ 

From continuity and the condition that Ag = A,^ the exit axial velocity is 



the tangential component is then 




and the component g is 


Vu, 8 = 


U^) 


T^8 


The rotor exit angles are 


^8 



O! 


8 



and 



(E54) 


(E55) 


(E56) 


(E57) 


(E58) 


(E59) 


(E60) 


(E61) 
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Then 


Mg = 



and 




The rotor exit fluid properties are 


and 



P8 = P8 


y 


— M 


\y8/(yQ-i) 


(E62) 


(E63) 


(E64) 


(E65) 


Exit Stators 

The velocity diagram for the exit stators is given in figure 13. The exit stator is 
specified as impulse, that is, Pg = pg. The exit stator total-pressure-loss coefficient 
Kg, is calculated in a manner similar to that of the single-stage turbine. Equation (C27) 
of appendix C gives the expression for Kg. The total pressure at the stator exit is 

Ps = Pg - Kg(p» - p)g (E66) 


Then knowir^ Pg and pg, the exit axial Mach number is 



(E67) 
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The calculated value of Mg is then compared with the value of Mg specified in the 
main program. K the values are different, an iteration is conducted based on an adjust- 
ment of the starter value of M c and the second-stage axial Mach number M „ 

Xj D Xy I 

according to 


M, 


X, 5|n+l 


(prescribed) 

■'^S(calculated) 


and 



n+1 



5 


n+1 


This procedure is repeated until 


(^s|n+l 



< 0.001 


(E68) 


(E69) 


(E70) 


From the final value of exit stator exit axial Mach number, the exit fluid properties 

are 


T^ = T^ (adiabatic stator) 


(E71) 


and 



1 + i (y - 1 )m2 

2 


The exit stator exit fluid density is 




/ 144 p \ 

V RT /g 


Then the exit stator exit axial velocity is calculated as 


(E72) 


(E73) 


(E74) 


Vg = (mV 7R gT)g 


61 



Overall Turbine Efficiency 


A calculated value of total-to-total efficiency based on blade-row losses is then de- 
veloped to replace the assumed value initially used in the main program. The determi- 
nation is based on the enthalpy- entropy diagram for the two-stage turbine shown in fig- 
ure 21. The total-to-total efficiency is defined as 





(E75) 


The actual power Prp is made up of the first- and second-stage power inputs, which for 
the case of no leakage and full admission, is 


Pt - 


4-6 


+ P, 


6-8 


where 


P 4-6 - %, 6 ^u, 5 


Wu,6) 


w. 


gJ 


and 


P 


6-8 


Ut, 8^Wu, 7 - - ^6) 

gJ 


(E76) 


(E77) 


(E78) 


where Zg is the interstage stator leakage factor defined as the ratio of the flow that by- 
passes the interstage stator (see fig, 3) to the through flow. The value of Zg is a pro- 
gram input. It is assumed that the interstage leakage flow does not affect the stage 
velocity diagrams or flow areas. From equation (E39), equation (E76) can be rewritten 
as 


Pt g 


W 


u, 


c - W c + 

5 u, 6 


(TRM)(W^_,-W„_8)(1 - Zg) 


W 


_4 

gJ 


(E79) 


The ideal power for the two- stage turbine based on total-to-total conditions is (from 
21 ) 
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(E80) 


^T, id “ " ^4^4 “ ^S, id^ 


The isentropic total-to-total enthalpy drop is obtained using the ADIABT subroutine of 

f 

appendix G with the inputs, p^, T 4 , Pg, — 
total-to-total efficiency becomes 


77 ^, and ?7 = 1.0. Thus, the overall turbine 


U, 




’T ,6 


5 - 6 + (TRM)(W„, 7 - s)(l - Zg) 


P'i - •'s, id>SJ 


(E81) 


The total-to-total efficiency as given by equation (E81) corresponds to the definition of 
efficiency used in the main program. This value of efficiency is then adjusted to account 
for the turbine-to-fan leakage which reduces the actual power output by reducing the 
mass flow. 

The two- stage total-to-total turbine efficiency adjusted for leakage is then 


'^T, leak “ ^T^^ " ^5^ 


(E82) 


The value of in equation (E82) is defined as 


w 


^5 = 


leak 


w> 


(E83) 


This expression is obtained using the method of appendix F. 


Turbine Annulus Flow Areas 


The first-stage stator exit annulus area for full admission is given by 


144 w^ 
^x, 5^5 


(E84) 


The rotor inlet annulus area which includes the turbine leakage term (to give the same 
velocity across the station) is given by 
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A., = A.(l - ZA 


(E85) 


and as previously stated, rotor outlet annulus area is 


^6 ~ 


(E86) 


The second- stage stator exit annulus area is 


^7 ~ 


- ^5) 
^7^x, 7 


(E87) 


The second-stage rotor inlet area is the same as the inlet stator exit annulus area, 
and it was also prescribed that rotor inlet area equals rotor exit area, so that 


As - A,^ 


(E88) 


Exit stator exit annulus area is obtained from 






(E89) 


Turbine Design Parameters 

For a multistage turbine, the overall speed-work parameter as defined in refer- 
ence 14 is 


A = 


U. 


A V 


u 


(E90) 


Then for a two-stage turbine. 


AV=W .-W c+W „ 
u u, 5 u, 6 u, 7 


W. 


u, 8 


and equation (E90) can be rewritten as 
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(E91) 


Un 


X = 


(W 1--W c+W „-W o) 
' u, 5 u, 6 u, 7 u, 8' 


The blade- jet speed ratio is (where TRM = 1, for simplicity) 


Un 


V -■ 




(E92) 


where Vj is the ideal jet velocity corresponding to the ideal expansion from inlet to exit 
static conditions across the turbine. Equation (D58) of appendix D is the expression used 
for obtaining Vj . 
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APPENDIX F 


LEAKAGE ANALYSIS 

The principal problem with leakage seals for tip-turbine fan-rotor wheels (see 
fig. 3) is the difficulty in maintaining close seal clearances due to poor dimensional con- 
trol inherent in this configuration. Seals applicable for such use could probably be of the 
abradable labyrinth or honeycomb types which allow for local rubbing contact without 
seizing or destruction of the seal. 


Leakage Flow Factor 

Leakage of the primary flow from the turbine to the tip region in front of and behind 
the fan rotor not only reduces turbine power output and thrust but also deteriorates fan 
performance. Thus, it is important to be able to predict the magnitude of the leakage 
flow for the geometry in question. This analysis will consider only the leakage for the 
first-stage stator because it has a direct flow path to the inlet tip region of the fan rotor 
as shown in figure 3. 

The basic relation for the leakage flow parameter was obtained from the relations 
given in references 17 and 18. The ideal leakage flow rate is given by 

^leak = ^<^f®V^ 

where A is the seal clearance geometry area (equal to irdp, ^cl), p is the density of the 
gas at the seal inlet, is the orifice flow coefficient, p' is the seal inlet total pres- 
sure, and e is the expansion function. Rewriting equation (FI) to form the leakage flow 
factor results in 


Ap’ 



(F2) 


In this equation, e was determined from experimental data on a round sharp-edge ori- 
fice. Values of e are tabulated in table I of reference 17 and plotted against p„.^/p' in 
figure 7 of the same reference. These values of e were determined based on a constant 
value of flow coefficient. Also, reference books give values for sharp-edge orifices 
from 0. 60 to 0. 624; an average value, Cj = 0. 615, was used in conjunction with the ex- 
pansion function e for reference 17 (single labyrinth). 
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Comparison of resultant plots of the product Cj and e of references 17 and 18 pre- 
sented in figure 22 showed very good agreement. Thus, for calculation purposes, the 
flow parameter based on equation (F2) was used in which the effective orifice coefficient 
Cj is set equal to 0, 615, 

For computer program use, it is desirable to have an analsdical variation for the ex- 
pansion function. It was found that the variation of e in reference 17 could be expressed 
reasonably well by the relation (for B/cl = 1,0) 


e = 




(F3) 


where e^„^ varied with the number of chambers in the seal. An expression for e„„,^ 

m 3-x 

as a function of the number of labyrinth restrictions was deduced from the data 

of reference 17 and given as 


e 

max 


0,925^ 


0. 50 


n 


leak , 


(F4) 


Variations in due to other choices of B/cl are handled by varying the pressure 
ratio exponent in equation (F3) in the following manner 


e = 



(F5) 


The variation of the exponent b with element width-to-clearance ratio B/cl was ob- 
tained from reference 17 and is shown plotted in figure 23, 

The final expression for the seal leakage flow rate function is obtained using equa- 
tions (F2) to (F4) as 



(F6) 
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Seal Pressure Ratio 


Since the approach velocity upstream of the seal orifice is expected to be quite low 
(relatively small flow rate and large area), the total pressure is essentially the local 
static pressure. Thus, for the inlet seal, in terms of the program nomenclature (see 
fig. 3) 


P’ = P5 (F7) 

and the temperature T' in equation (F6) is the total temperature, that is, T' = Tg. 

The seal exit pressure p^^ is defined as the local static pressure at the inlet bell- 
mouth at the point of the seal. This pressure is given by the expression 

Pex = Pq “ ^leak^Pio " ^10^ 

where Pjq ~ ^10 average value of the inlet velocity head at the mean fan radius 

and is a constant obtained from potential flow solutions of typical inlet bellmouth 

sections for lift fans (e.g., ref. 19). Calculations indicated a representative 
value of 1.2, since the axial velocity at the tip of the rotor inlet is generally larger than 
at the midpoint. The average static pressure in equation (F8) is obtained from (for 
X = 1. 4) 


Pio 


1 


^1 + 0. 2 M 



3. 5 


where p^Q and Mjq are inputs from the main program. 


(F9) 


Leakage Rate 

Once C^e and expansion pressure ratio Pgj^/P’ are obtained, the leakage flow fac- 
tor equation (F6) can be rewritten in the following form by inserting A = n^d_ . cl and 

X* j t 

noting that p’ and T' are replaced by pg and Tg, respectively, 
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(FIO) 


Since the actual value of runnir^ clearance will not generally be known, the value of cl 
was prescribed by the empirical expression 


cl.Ceidp_, (FlI) 

_4 

The factor in practice has values of around 2. 5x10 for the rub-in type of seal. 

The leakage ratio based on turbine flow rate is then 


2 (F12) 

»4 

The turbine-to-fan leakage for the two-stage turbine design (fig. 3) was treated in 
tie same manner as the single-stage design. Inherent in the two-stage turbine is a 
i^her inlet stator exit static pressure and possible mechanical design problems (larger 
learance gap) that would increase the leakage flow rate. These conditions may warrant 
somewhat larger value of for use in equation (Fll). 

The leakage flow not only affects the performance of the turbine by reducing the 
rimary weight flow, but also affects the performance of the fan. The leakage flow dis- 
arbs the fan flow at the tip region which causes a reduction in fan efficiency r/^,. The 
‘ffect on 77p, is not included specifically in the program; however, the value of 7]y 
which is an input) can be chained to account for the leakage flow if the variation is 
snown. 
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APPENDIX G 


FLUID PROPERTY SUBROUTINES 


In addition to the major subroutines required by the main program, there are 12 
smaller routines that describe the fluid properties of both primary (turbine) flow and sec- 
ondary (fan) flow. The subroutines are written in terms of a basic routine called POLY. 
This subroutine is a working tool for the other routines because it describes enthalpy, 
entropy function, and specific heat as a function of temperature along with the reverse 
situation of temperature as a function of enthalpy or entropy function. The subroutines 
presented in this section are general in that they are valid for both fuel and air mixtures 
and air, A list of the fluid property subroutines described in this appendix are 

(1) Function POLY (7) Function SPHTCP 

(2) Function FULAER (8) Subroutine ADIABT 

(3) Function TFROMS (9) Function GAMMA 

(4) Function GASCON (10) Function ENTROP 

(5) Function TFROMH (11) Function VISCON 

(6) Function ENTHAL (12) Function FMM 


Function POLY 


The thermodynamic properties (specific heat, enthalpy, and entropy function) for air 
and for stoichiometric fuel and air mixtures for several fuels have been tabulated by Hall 
and Weber in reference 20. These properties for air and for a stoichiometric mixture of 
air with hydrocarbon fuel have been prc^rammed as polynomials in absolute temperature 
(*^R) in function POLY. 

Formulas are also derived in reference 20 that permit the calculation of the ideal 
fuel-to-air ratio for any given inlet air temperature and final combustion temperature. 
Ideal fuel-to-air ratio corresponds to complete combustion assuming no dissociation. 
Procedures are also indicated for calculating the thermodynamic properties of fuel-air 
mixtures with fuel-to-air ratios less than stoichiometric. 

POLY yields gas properties based on relations for the specific heat of a mixture, 
entropy function of a mixture, and enthalpy of a mixture. These relations are 


C„ = C 


p,a 


+ G»//(Cp) 


(Gl) 


cp = (p^ + G\l/{cp) 


(G2) 


h = h^ + G;//(h) 


(G3) 
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where the term G is a function (modified for incomplete combustion) of the efficiency of 
burning and the fuel-to-air ratio. The expression for G is 



The value of ri^ is the same as that used in the burner. The functions Cp (p^, 
hg^, 3^re all tabulated in reference 20 and are expressed as poly- 

nomials in the absolute temperature (°R) in function POLY. The POLY functions are: 

For air 


POLY (T, 1) = h 
POLY (T, 2) = cp 
POLY (T, 3) = C 
POLY (h, 7) = T 
POLY (^, 8) = T 

For JP fuel in air 


POLY (T, 4) = i//(h) 
POLY (T, 5) = x}^(<p) 
POLY (T, 6) = i//(Cp) 


Function FULAIR 

This function evaluates fuel-to-air ratio as a function of burner inlet and outlet tem- 
perature and burner inlet fuel-to-air ratio. The function FULAIR gives only the ratio of 
fuel added in the burner to the air (not the fuel-to-air mixture) entering the burner. 

Then the total burner outlet fuel-to-air ratio will simply be the sum of the input and the 
computed fuel-to-air ratio. An additional input required is the burner efficiency 
The function FULAIR is given as 



Function FULAIR jTj, Tg, - 




where the calculation procedure is 


(l)fi- =0.02 
'2T 



(3) hg = ENTHAL 




(4) hj = ENTHAL 



(5) 



(hg - hj) 


1 +1 


77^^HV(h2 - H°) 


(6) If 




1. 0x10"® 


0 




and return to step (2). 


If the difference is <0 then do the following 


(7) FULAIR = 



(8) Return to main program 
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Function TFROMS 


This function evaluates temperature from entropy function using an iterative proce- 
dure. The function TFROMS is given as 


Function TFROMS 



f 

a 




The calculation procedure for this function follows: 


(1) G = 

(2) T^ = POLY(<p^, 8) 

(3) <p^ = POLY(T^, 2) 

(4) T^ = POLY(<^^, 8) 

At this point in the program, solve for T by the false position method which obtains 

oi 

approximations for T. Thus, for the mixture, TFROMS = T . 

A 3 , 




POLY(T^, 2) + (G) POLY(T^, 5)' 

<x a. 



Function GASCON 

This function evaluates the gas constant as a function of fuel-to-air ratio. 

Function GASCON 
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and the GASCON = R. 


R 


53.349 + 55.088- 
a 


1 + 


Function TFROMH 

This function evaluates temperature from enthalpy and the fuel-to-air ratio of the 
mixture. This function also requires the burner efficiency and the function POLY. An 
iterative process is required to evaluate the mixture temperature. The function 
TFROMH is given as 


Function TFROMH 



f 

a 



where the calculation procedure is 



(2) = h^ 

(3) T^ = POLY(h^, 7) 

(4) i//^(h) = POLY(T^, 4) 



(6) T^ = POLY(h^, 7) 

Now, solve for T^^ using the false position method, which obtains approximations for 

T . Then, for the mixture, TFROMH = T . 
a ’ ’ a 
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Function ENTHAL 


This function gives the enthalpy as a function of temperature, burner efficiency, and 
fuel-to-air ratio. Function POLY is used in the procedure: 


Function ENTHAL |t, 77 ^, — 


(1) h^ = POLY(T, 1) 

(2) i//j^(h) = POLY(T, 4) 



(4) + 

(5) ENTHAL = h^^. 


Function SPHTCP 


This function gives specific heat at constant pressure as a function of temperature, 
fuel-to-air ratio, and burner efficiency. 


Function SPHTCP 



f 

a 



(1) Cp^ ^ = POLY(T, 3) 



(3) i//(Cp) = POLY(T,6) 

(4) Cp,M = Cp,a-^G'^(Cp) 

(5) SPHTCP = Cp j^. 
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Subroutine ADIABT 


This subroutine evaluates an expansion or compression of a gas taking into account 
initial temperature and pressure, fuel-to-air ratio, adiabatic efficiency, and final pres- 
sure. Outputs of the subroutine are final temperature and enthalpy change. The enthalpy 
change for an expansion process will be a negative value. 

(1) <p^ = ENTROP(^Tj, i, 

(2) R = GASCON^ij 

(4) T2_y=TFROMs(^2, i, 

(5) hj = ENTHAl(^Tj, i, 

(6) h2_id = ENTHAL(T2^ia. L, 


For a compression process (p 2 > Pi) 

(7) Ah = 


^^2, id " 


For an expansion process (p 2 < Pi) 


(8) Ah = ~ ^1^ 

The final enthalpy is 


(9) h2 = hj + Ah 
and the exit temperature 


(10) T 2 = TFROMH^h2, 
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Function GAMMA 


This function evaluates the ratio of specific heats according to the perfect gas 
definitions . 


Function GAMMA 



f 

a 



(1) Cp = SPHTCP^T, K 


(2) GAMMA = 

(C - [GASCON ]/778. 16) 


Function ENTROP 

This function yields entropy function as a function of temperature, burner efficiency, 
and fuel-to-air ratio. 


Function ENTROP 



(1) <p^ = POLY(T,2) 



(3) xpicp) = POLY(T,5) 

( 4 ) <Pyi" 

(5) ENTROP = (p^ 


Function VISCON 

This function determines the viscosity of air from a table of viscosity against tem- 
perature. The function VISCON is an interpolation procedure (see appendix H for pro- 
gram listing). 
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Function FMM 


This function solves for Mach number for an input value of flow factor. The equa- 
tion is solved by an iterative procedure. This expression is 

FF 

Inputs are FF, y, g and the output is M. 
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APPENDIX H 


program listing 


SIBFTC VTOLl 
C 

COMMON/DATl/ ADAS, OMEGAS, 

1 DTLOSS, ETABR,ETAF,ETAN,ETAT,FKV, 

2 FM11,FM2 , FM2B,FM35. FNB , FNS , FPR , H2P , HHV, P I , 

3 P1413, PAMB, POPS, PP2B,R,RHRT, RING, TAMB, 

4 T1,T4,UFT, WIOBRtWFWAI ,WLBAR,CT, 

5 CM,FKSB,FKSN,XIALS,CCM,AEXP, 

6 THARC.COQS, AR,F14 

COMMON/DAT2/ ITURB, lETAT, IWR I TE , P I E , FK, TRM, X I AL 6 , , 

1 FLLC,WIS2,P20P1,W1S1,FN3 

C 

DIMENSION FN3J(50),P2nPlJ(50),WlJ(50),FMSJ( 100) 

DIMENSION ALPH5JI 100 ) ,DTTJ( 100) , FMR5J I 100) , TBHJ ( 100 ) , 

1 VLJJJ(IOO), 

1 WICMP(IOO) ,W1TST(5),FMSTST(5),ALPTST(5) 


IF lOUT IS READ IN AS 1, COMP WILL BE REAPEADETLY 
CALLED AND OUTPUT OBTAINED 
IF lOUT IS READ IN AS 0, 

MINIMUMS WILL BE OBTAINED, AND MOST OF THE 
OUTPUT WILL NOT BE PRINTED - IE. , ONLY 
MINIMUM CASES WILL BE PRINTED OUT. IN THIS CASE 
lOUT WILL BE USED TO SUPPRESS OR NOT SUPPRESS 
PRINTING 

IN THE CALLING VECTOR, IF lOUT = 0, DO NOT WRITE 
OUTPUT 

IN THE CALLING VECTOR, IF lOUT = 1, WRITE OUTPUT 
JP7 IS THE INDEX IN THE FMSJ PROFILE 
CORRESPONDING TO FMS = .7 
THE VALUES IN THE FMSJ INPUT VECTOR MUST BE IN 
INCREASING ORDER 

SUB. COMP INTERATES ON ETAT, AND INITIAL GUESSES 
FOR ETAT AND ETAN ARE NEEDED, FOR ANY TWO 
RUNS USING IDENTICAL INPUT, EXCEPT FOR 
DIFFERENT ETAT AND ETAN INPUTS, THE RESULTS 
SHOULD BE IDENTICAL, BUT INFACT, COMPARE TO 
ONLY FOUR OR FIVE DECIMAL PLACES 
IF IPLDTT = 0, AND IDUT = 0, THEN PLOTS WILL BE 
DRAWN 

IF IPLOTT = 1, PLOTS WILL NOT BE DRAWN 


1 READ(5,51) ITURB, IETAT, IWRITE , lOUT, JP7, I WGHT, IPLOTT 
C 

IWRTSV = IWRITE 
C 

IF(IETAT) 160,160,57 
160 GO T0( 2, 161, 3,2) , ITURB 
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C 

3 REAIH5.50) FK, TRM, XI AL6, FLLC 
WRITF(6,5) 

WRIl E( 6,5 55 I ITURB,IETAT, IWRITE, TOUT, JP7, IWGHI, IPLOTT 
WRITE(6,60) FK,TRM,XI AL6,FLLC 
C 

GO ro 88 
57 WRI 1E( 6,58) 

WRITE! 6,555) I TUR8 , I E T AT , I WR I T6 , lOUT , J P 7, I WGHT 
GO TO 38 

161 WRIIE(6,162) 

WRITE(6,555) ITURB,IETAT,IWRITE, lOUT, JP7, I WGHT 
GO TO 88 
2 WRITE(6,54) 

WRITE(6,555) I TURB , I ET AT , I WR ITE , I OUT , JP 7 , I WGHT 
C 

88 IF(IWRITE) 56,56,87 
87 REA1H5,50) W1S2 
WRITE(5,89) W1S2 


56 REAU(5,50) TLH, OMEGAS, 

1 DTLOSS, ETABR,ETAF,ETAM,FTAT,FKV, 

2 FM11,FM2, FM2B,FM3S, FNB , FMS , FPR , H2P , HHV , P 1 , 

3 P1413, PAMB. PP2B,R,RHRT,RING, TAMB, 

4 T1,T4,UFT, W10BR,WFWA1 ,WLPAR,CT, 

5 CM.FKS8, FKSN,XIALS,CCM,AEXP, 

6 TH4RC,C00S,AR,F14,ZIL 

NOTE - IF IWRITE DOES NOT EQUAL 2, THEN FN3 IS 
NOT used. BUT ONE VALUE OF FN3 MUST BE READ 
INTO THE PROG. TO SATISY THE DO LOOP 


READ(5,102) JFN3,(FN3J(J),J=1,JFN3) 

REAi.'(5 , 102) JP20P1,( P20P1J ( J ) , J= 1 , J P20P 1 ) 

REA0(5,102) Jwl, ( WIJ ( J ) , J=l, JWl ) 

REAU(5,102) JFMS,(FMSJ(J),J=1,JFMS) 

107 WRITE{6,101) (FN3JU ) , J=l, JFN3) 

WR1TE(6,103) {P2GP1J(J), J=1, JP20P1 ) 

WRITE! 6,104) ! W IJ ! J ) , J= 1 , J W 1 ) 

WRITE!6,105) !FMSJ!J ) , J=1, JFMS) 

C 

WRITE! 6,66) 

WRI TE! 6,66) 

C 

WRITE!5,55) TLH, OMEGAS, 

1 DTLOSS, ETABR,ETAF,ETAN,ETAT,FKV, 

1 FM11,FM2 , FM2B,FM3S, FNB, FNS , FPR , H2P , HHV , P 1 , 

2 P1413, PAMB, PP2B,R,RHRT,RING, TAMB, 

3 T1,T4,UFT, W10BR,WFWA1,WLBAR, CT, 

4 CM,FKSB,FKSN,XI ALS,CCM,AEXP, 

5 THARC,COOS,AR,F14 
C 

ADAS = 1. + .1744*TLH 
C 

DO 92 192 =1,JP20P1 
P20P1 = P20P1J!I92) 
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1> JWi 


C 


00 93 193 = 

WlSl = W1 J ( 193) 

DO 91 191 = 1, JFN)3 
FN3 = FN3J(I91) 

W1 = WlSl 
IWRITE = IWRTSV 

PIE = 3.1415926 


IFdOOT) 70,70,71 
71 DO 94 194=1, JFMS 

FMS = FMSJ( 194) 

W1 = WlSl 
IWRITE = IWRTSV 

94 CALL C0MP( lOUT. IWGHT,FMS,ALPH5,DTT,FMR5,TBH,W1,FN, VL) 
GO TO 91 


70 DEG75 = 75.*PIE/180. 

TO OBTAIN VECTORS TO PLOT 

DO 25 J=1,JFMS 

IFU-JP7) 21,22,21 
22 lOUT = 1 
21 FMS = FMSJ(J) 

W1 = WlSl 
IWRITE = IWRTSV 

CALL COMP ( I OUT , IWGHT , FMS , ALPH5 , DTT , FMR5 , TBH , W1 , FN , VL ) 
ALPH5JU) = ALPH5 
DTTJ(J) = DTT 
FMR5JIJ) = FMR5 
TBHJU) = TBH 
WICMP(J) = W1 
VLJJJ(J) = VL 
25 lOUT = 0 


TO FIND MINIMUM Wl 


JMIN = 1 

DO 27 J=2,JFMS 

IFIWICMPI JMIN) - WICMP(J)) 27,27,26 

26 JMIN = J 

27 CONTINUE 


TO SEE IF WICMP(JMIN) IS AN END POINT 
OR AN INTERMEDIATE POINT 
TO FIND WI MINIMUM, TWO POINTS 

ARE TAKEN ON EITHER SIDE OF WICMPIJMIN) 

IFIJMIN - 2) 30,31,31 
30 FMS = FMSJIl) 
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WITSTW = WlCMP(l) I 

ALPTST(l) = ALPH5J(1) 

JATST = 1 
GO TO 3A 

31 IFCJMI'g - JFMS) 32,33,33 

33 FMS = FMSJ(JFMS) 

WlTSTk^ = WICMP(JFMS) 

ALPTST(JFMS) = ALPH5J(JFMS) 

JATST = JFMS 

GO TO 34 

32 DEL13 = (FMSJ(JMIN) - FMSJtJMIN - l))/3. 

DEL23 = (FMSJ(JMIM * 1) - FMSJ ( JMI N ) ) /3 . 

JJ. = 0 

SIGN = 1. 

DEL = DEL13 

DO 35 1=1,2 
SIGN = - SIGN 
lOUT = 0 

' D336J=1,2 

JJ = JJ +1 

FMS = FMSJ(JMIN) + S I GN*DEL *FLOAT I J ) 

W1 = WlSl. 

IWRITE = IWRTSV 

CALL COMPt TOUT, I WGHT , FMS , ALPH5, DTT , FMR5 , TBH, W1 , FN , VL ) 
. ALPTST(JJ1= ALPH5 
FMSTSTUJ) = FMS 

36 WITST(JJ) = W1 ‘ 

35 DEL = DEL23 

FMSTST(5)= FMSJ(JMn') 

W1TST(5) = WICMPIJMIN) 

ALPTST(5) = ALPH5J(JMIN) 

TO FIND A SECOND AND MORE REFINED W1 

JMIN5 = 1 ■ 

DO 37 J=2,5 

IF(W1TST( JMIN5) - WlTSTUn 37,3.7,36 
38 JMIN5 = J t 

37 CONTINUE 
JATST = JMIN5 

FMS = FMSTSTt JMIN5) 

WITSTW = W1T5T1JMIN5) 

34 IF ( ALPTST ( JATST) - DEG75) 28,28,29 
28 lOUT = 1 

WRITE(6,9) FMS 
W1 = WlSl 
IWRITE = iwRTSV 

CALL :0MP( IQUT, IWGHT, FMS,ALPH5,DTT,FMR5, TBH, W1 ,FN, VL) 
lOUT = 0 
GO TO 39 


ALPTSTI JATST) IS GREATER THAN 75 DEG. 

TO FIND THE INTERVAL IN WHICH ALPH5J IS 

75. DEG., USING THE READ IN, COARSE 




INTERVALS 
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c 


29 IFtJMIN 
73 JMINl = 


- 1) 74,74,73 
JMIN - 1 

OD 399 J=1,JMIS1 
JJ = JMIN - J 

IF{ ALP^5J{ JJ) - DEG75) 40,40,399 

40 J751= JJ’ 

J752 = JJ + 1 
GO TO 41 

399 CONTINJE 

INTERVAL NOT FOUND TO TRY J LARGER THAN JMIN 

74 IFIJMIV - JFMS) 75,76,76 

75 JMIN2 = JFMS - JMIN 
JJ = JMIN 

DO 42 J =1,JMIN2 
JJ = JJ + 1 

IF( ALPH5J ( JJ ) - DEG75) 43,43,42 
43 J751 = JJ - 1 
J752 = JJ 
GO TO 41 
42 CONTINUE 

76 WRITF(S,44) FMS J ( JMI N ) , W 1C MP ( JM I N ) , ALPH5 J ( JM I N ) 

GO TO 39 

TO COMPUTE FMS WHICH CORRESPONDS TO ALPH5 = 75 DEG. 

41 FMSll = (DEG75 - ALPH5 J I J751 ) ) / ( AL PH5 J ( J 752 > - ALPH5J ( J751 ) > 

FMS = FMS11*(FMSJ( J752) - FMSJU75D) + FMSJ(J751) 

lOUT = 1 
WRITE(5,45) FMS 
W1 = WlSl 
IWRtTE = IWRTSV 

CALL C3MP{ lUUT, IWGHT , FMS , ALPH5 , DTT , FMR5 , TBH , W1 , FN , VL ) 

TO MAKE CALCOMP PLOT 

TO NON-DIMENSIONALIZE THE POINTS TO bE PLOTTED 
BY DIVIDING BY THE VALUES AT FMS = .7 
39 DO 46 J=1,JFMS 

ALPH5JIJ) = ALPH5JU )/ALPH5J( JP7) 

OTTJ(J) = DTTJ(J)/DTTJ( JP7I 
FMR5JIJ) = FMR5JI J)7FMR5J( JP7) 

TBHJ(J) = TBHJl J)/T3HJUP7) 

WICMP(J) = wlCMPI J)/W1CMP( JP7) 

46 VLJJJ(J) = VLJJJI J)/VLJJJ( JP7) 

IF(IPLOTT) 200,200,91 

200 CALL HALPLTI JFMS,FMSJ, ALPH5J,0TTJ,FMR5J,TBHJ,WICMP, 

1 VLJJJ.FPR, Ti,T4,UFT,Pl, ITURB) 

91 CONTINUE 
93 CONTINUE 

92 CONTINUE 
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IF(IWRITE) 56,1,1 


FM2 =E12.5,5X, 


FNS =E12.5,5X, 


R =E12.5,5X, 


T1 =E12.5,5X, 


50 F0RmaT(9E8.5) 

51 F0RMAT(18I4) 

555 FORMAT! 1H0,8H ITLIRR = I 4, 5X , 7H I FT AT = I 4, 5X , 8H I WR I TE =14, 5X, 

1 5HI0UT =I4,5X,6H JP7 =I4,5X,8H IWGHT =14) 

5 FORMAT! 1H1.27H TWO STAGE TURBIfTE WITH OGV) 

54 FORMAT! IHl, 30H SINGLE STAGE TRUBINE WITH OGV) 

58 FORMAT! IHl, 02H INPJT FOLLOWS - TURBINE ROUTINE IS SKIPPED SINCE 
1 THERE IS NO IMTERATIGN ON ETAT) 

7H TLH =E12.5,5X,7HOMEGAS=E12.5,5X,7HDTLOSS=E12.5, 

/IH , 

7HETABR =E12.5,5X,7H ETAF =E12.5,5X, 
7H ETAN =E12.5,5X,7H ETAT =E12.5/1H , 

7H FKV =E12.5,5X,7H FMll =E1?.5,5X,7H 

7H FM2U =E12.5,5X,7HFM3S =E12.5/1H , 

7H FNB =E12.5,5X,7H 
7H FPR =E12.5,5X,7H H2P =E12.5/1H , 

7H HHV =E12.5,5X,7H PI =E 1 2 . 5 , 5X , 7HP 1 41 3 =E12.5,5X, 

7H PAMB =E12.5/1H , 

lyx, 5X,7H PP2B =E12.5,5X,7H 

7H RHRT =E12.5,5X,7H RING =E12.5/1H , 

7H TAM6 =E12.5,5X,7H 
7H T4 =E12.5,5X,7H UFT =E12.5/1H , 

7HW10BR =E12. 5, 5X, 7HWFWA1 =E12.5,5X, 
7HWLBAR =E12.5, 5X,7H CT =E12.5/1H , 

CM =E12.5,5X,7H FKSH =E12»5,5X,7H FKSN =E12.5,5X, 

7HXIALS =C12.5,5X,7H CCM =612.5/1H , 7H AEXP =E12.5, 
5X.7HTHARC =E12.5,5X,7H COOS =E12.5,5X,7H AR =E12.5, 

5X,7H F14 =E12.5) 

60 F0RMAT!1H , 7H FK =E12.5,5X, 7H TRM =E 1 2 . 5 , 5X , 7HX I AL6 =E12.5, 

1 5X,7H FLLC =E12.5) 

89 FORMAT ( IHO, 7H W1S2 =E12.5) 

102 FORMAT! 14, ! 9E 8.5 n 

101 FORMAT! 1H0,20H FN3 PROFILE FOLLOWS/ ( 8E 15 .5 ) ) 

103 format ( 1H0,22H P20P1 PROFILE FOLLOWS/ ! 8E 1 5 . 5 ) ) 

104 FOR’'<AT ( IHO, 19H W1 PROFILE FOLLOWS/ ! 6 E 1 5 . 5 ) ) 

105 FORMATIIHO, 20H FMS PROFILE FOLLOW S/ ! 8E 1 5 . 5 ) ) 

162 FORMAT! IHl, 36H SINGLE STAGE TURBINE - STATORLESS) 

66 FORMAT (IHK) 

44 FORMAT ( IHO, 104H FOR THE GIVEN FMS RANGE, ALPH5 IS NEVER LESS THAN 
175 DEGREES - FOR THIS RANGE, THE MINIMUM W1 FOLLOWS/ 

2 7H FMS =E12.5,5X,7H Wi =E 1 2 . 5 , 5X , 7HALPH5 =E12.5) 

45 FORMATIIHO, 50H THE CASE CORRESPONDING TO ALPH5 = 75 DEG. FOLLOWS/ 

1 IH ,7H FMS =E12.5) 

9 FORMATIIHO, 39H OUTPUT FOR THE MINIMUM Wl CASE FOLLOWS/IH , 

1 6H FMS =E12.5) 

END 


55 FORMATIIH 
1 
2 

3 

4 

5 

6 

7 

8 
9 
1 
2 

3 

4 

5 

6 

7 

8 
1 
2 


7H 
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$IBFTC TUi<Bl 


C 


2 

3 

4 

5 

6 

7 

8 
9 
1 
2 
3 
1 


MAY 1968 

SUBi^OJTINE TURB1(T4,WA,P4,PSS, OFT.DT, FMS.ETAN, 

F3, FMIO ,P10,PAMB, 

TS5,V5,FM5,ALPH5, VZ5, VU5,P5,U,BETA5, 

R5, FMR5,R6»ETAB.TS6, VZ6,FMZ6,FMR6,BETA6, 
RU6,QMEGA, ALPH6f P6,PS,T6tFMSCAL tWLEAK, 

PHI ,PS10,SPR,ZLEAK,ETATL1, ETAT, 

FLAW,F,'-IUQR,£TABR, 

CM , FKSB,FKSN,XIALS,CCM, 

AZ5tFM6,FMZ5,PSFAN,RE,RU5,TS,TSS, V0,VU6,XIM 
XIN, XI P,XIS,Z£TAAS,ZETAB,ZErAIM, ZETAP.ZETASS 
ZETAS, TSP.DHTI I , C P6 , DH , AEXP , DTR , 

TBH.THARC, COOS tAR, F14) 


WRITTEN FOR HANK HALLER MAY 1968 

SINGLE STAGE TURBINE WITH OGV 


MAIN REFERS TO THE CALLING PROGRAM 
WFWA3P IN MAIN IS F3 IN TURBDP 
FMll IN main is FMIO IN TURBUP 
Pli IN main is PIO IN TURBDP 
PSCHK IN MAIN IS PS IN TURBOP 
PSll IN MAIN IS PSIO IN TURBDP 


G = 32.174 
FJ = 778. 

PIE = 3.1415926 
PIE180 = 180. /PIE 

RE = 2 4.*W4/0T/V I SCI T4 ) 
REXP2 = (RE/1.5E6)**. 2 
XIP = XIALS/REXP2 
XIS = FKSN*XIALS 
PS5 = PSS 
PS6 = PSS 

GASKON = 6ASCUNIF3) 

R = GASKON 
FM5C = 2.*FJ/R 
DLH8C = 2.*G*FJ 
RH06C = 144.*PS6/R 


THE INITIAL M5 VALUE FOLLOWS 


FM5 = 1. 

KCT = 0 
C 

11 IFIFM5 - 1.16,7,7 
6 XIM = 0. 

GO TO 8 
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7 XIM = D.0075*{FM5*FM5-FM5) 

8 XIN = XIP + XIS + XIM 
ETAN = 1.- XIN 


CALL ADIABT(T4,P4,ETAN, F3. TS5 , PS5 , DH, ETABR ) 
GAM = 3AMMA(TS5, F3, ETABR ) 

FM5 = SORT! FM5C*ABS( 0H)/TS5/GAM) 


IF(KCT-1 )9, 10,10 
9 KCT = 1 

FM5SV = FM5 
GO TO 11 

10 IF{ ABS ( FM5-FM5SV)-.002 ) 14, 14, 12 
12 FM5SV = FM5 

KCT = KCT+1 
IF(KCT-100)11,11,15 

15 WRITE(S,16) FM5,FM5SV,ETAN 

16 FORMAT! IH , 36H IN SUB TURBDP, FM5 CAN NOT BE FOUND/IH , 

1 6H FM5 =E12.5,5X,8H FM5SV =E12.5,5X,7H ETAN =E12.5) 

BEGIN STATOR EXIT MACH NUMBER AXIAL COMPONENT FLOW 
FMZ5 CALCULATION 

14 CP6= SPHTCPI TS5,F3, ETABR ) 

FMZ6C = GAM*G*R 
GAMl = GAM-1. 

GAMll = GAMI/GAM 
GAM2 = GAM/GAMI 
GAM3 = 2. /GAMl 
GAM4 = (1. -GAMI/GAM 
GAM5 = -GAM4 

V5 = S3RT(2.*G»FJ*ABS(DH)) 

P5 = PS5* ( 1 .+0. 5*( GAM-1. )*FM5*FM5) **GAM2 
RHOo = 144.*PS5/TS5/R 
FMR5C = SQRT(GAM*G*R*TS5 I 

BEGIN FMZ5 COMPUTATION 


KMS = 0 

FMZ5 = FMS*CCM 
34 ALPH5C = FMZ5/FM5 
ALPH5= ARCOS (ALPH5C) 

VZ5= \/5*ALPH5C 
VU5= y/5*SIN(ALPH5) 

AZ5 = 144.*W4/VZ5/RH05 
RU5= VU5-U 

betas = ATAN2IRU5.VZ5 ) 

R5 ^ SQRT(VZ5*VZ5+RJ5*RU5) 

FMRS= R5/FMR5C 

TO COMPUTE THE BUCKET EFFICIENCY, ETAB 


C 


OELB =2.*BETA5 

ZETABC = PIE18O*0ELB/ 100. 
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ZETABS = 0,02b *■ 0.03A*ZETABC*ZETABC*SORT(ZETABC ) 
ZETAP = ZETABS/REXP2 
ZETAS = FKSB*ZETA8S 
IF(FMR5-1.) 42,43,43 

42 ZETAM = 0.0 
GO TO 44 

43 ZETAM = CM*{ FMR5*FMR5 - FMR5) 

44 ZETAB = ZETAP + ZETAS + ZETAM 
ETAB = 1.- ZETAB 

IF(rHARC - (360. - .01)) 60,61,61 

60 SPACE B TBH/COOS/AR 
ARCL = PI E*0T*THARC/360. 

FKPA = 1. - SPACE/3. /ARCL 
R6 = R5*S0RT(ETAB)*FKPA 
GO TO 62 

61 R6 = R5*S0RT(ETAB) 

62 OLHB = (R5*R5-R6*R6) /DLHBC 
TS6 = TS5+DLHB/CP6 

RHOb = RH06C/TS6 

VZ6 = 144.*W4/RH06/AZ5 

FMZ6S = SQRr(FMZ6C*TS6) 

FMZ6 = VZ6/FMZ6S 

FMR6 = R6/FMZ6S 

RU6 = SQRT(R6*R6-VZ6*'VZ6) 

BETA6= ATAN2 (RU6.VZ6 ) 

VU6 = RU6-U 

ALPH6 = ATAN2( VU6,VZ6 ) 

OELALP= PIE180*ABS(ALPH6 ) 

DLALPC = OELALP/100. 

ZETAAS = 0.025 + 0 . 0 34*DLALPC*DLAL PC*SQRT ( DLALPC ) 
ZETASS = ZETAAS* ( 1./REXP2 + FKS8) 


51 OMEGA = ZETASS*(1. + 0 . 5*GAM*FMS*FMS ) 

52 FM6 = FMZ6/C0S(ALPH6) 

,P6 = PS6* (l.+0.5*(GAM-l. )*FM6*FM6)**GAM2 
PS = P6-0MEGA*( P6-PS6 ) 

T6 = TS6* ( 1 .+0. 5*GAM1*FM6*FM6) 

FMSCAL = SQRT(GAM3*( ( PS/PSS)**GAM11- 1.)) 

IF( ABS(FMSCAL-FMS) -.001) 31,31,30 
30 KMS = KMS+1 4 

IF(KMS-25) 32,32,33 

32 FMZ5 = FMZ5*FMS/FMSCAL 
GO TO 34 

33 WRITE(6,35) FMS , FMSCAL , FMZ 5 

35 FORMAT! IH , 90H IN SUB. TURB68,FMZ5 CAN NOT BE FOUND - 
ICONTINUES USING THE LAST VALUE OBTAINED/IH , 

2 6H FMS =E12.5,5X,7HFMSCAL=E12.5,5X,7H FMZ5 

FMZ5 IS FOUND 


31 TS = T6 

TSS = TS/ ( l.+.5*GAMl*FMS*FMS ) 


PROGRAM 

=E12.5) 
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vs = FMS*SORT (GAM*R*G*TSS) 

PS10= PIO/ { 1 .+0.2*F»U0»FM10 ) »*3.5 
DELPS = 1.20*PAMB*(P10-PS10)/P10 
PSFAN = PAMG-DELPS 
SPR = PSFAN/PS5 
C 

PHI = 3.92S*SURT( l.-SPR**AEXP) 

DELLIT = 0.00U25*UFT 

WLEAK = PI E*DTR»DELL I T*S«RT( G/R)*PS5*0.615»PH1 /SORT! TS5) 

C 

A1 ZLEAK = WLEAK/W4 

CALL ADIAflT(T4, PA, 1.0, F3.TSP.PS, OHTTI ,ETABR) 

C 

C CHANGES MADE 3-7-69 

C 

IF(rHARC - (360. - .01)) 63,64,64 

63 ETATL2 = RU5 + RU6 - .5ttF14*U*U*C00S*SPACE/V5/ARCL/ALPH5C 
ETATLl=U*ETATL2/AbS{DHTTI)/G/FJ 

GO TO 55 

64 ETAFLl = U*(RU5 + RU6 ) / A ES ( DHTT I ) /G/ F J 
C 

65 ETAT = ETATL1*(1. -ZLEAK) 

FLA'^’ = U/(RU5+RU6) 

VO = V5/S0RT( E7 AN) 

40 FNUoR = U/Vu 
C 

c 

RETURN 

END 


1.IBFTC TURB2 


C 


C 

C 

C 

C 

C 

c 

c 

c 

c 


2 

3 

4 

5 

6 

7 

8 
9 
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FEBRUARY 1969 

SUBIOJTINE TURB2( T4,W4,P4,PSS, DFT,DT, FMS,ETAN, 

F3, FMIO ,P10,PAMB, 

TS5, V5,FM5,ALPH5, VZ5, VU5,P5,U,BETA5, 

R5 , FMR5,R6,ETAB,TS6, VZ6.FMZ6, FMR6,BETA6, 
RU6,0MFGA, ALPH6,P6,PS,T6,FMSCAL ,WLEAK, 

PHI ,PS10, SPR, ZLEAK, ETATLl, ETAT, 

FLAMjFNUBR.ETABR, 

CM,FKSH,FKSN,XIALS, 

AZ5 ,FM6,FMZ5,PSFAN,RE,RU5, TS, TSS,V0, VU6.XIM 
XIN,XIP,XIS, ZETAB,ZETAM,ZETAP, 

ZETAS,TSP,OHTTI, CP6, DH, AEXP , DTR , 

TBH, THARC, COOS ,AR,F14) 


, 


STATORLESS TURBINE SINGLE STAGE 
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MAIN REFERS TO THE CALLING PROGRAM 
WFWA3P IN MAIN IS F3 IN TURBDP 
FMll IN MAIN IS FMIO IN TURBDP 



onoo n on on o noon o o oooooo 


Pll IN MAIN IS PIO IN TUKBDP 
PSCHK IN MAIN IS PS IN TURBOP 
PSll IN MAIN IS PSIO IN TURBOP 


G = 32.174 
FJ = 778. 

PIE = 3.141B926 
PIE18D = 180. /PIE 

RE = 24.*W4/DT/VISC( T4) 
REXP2 = (RE/1.5E6)**.2 
XIP = XIALS/REXP2 
XIS = FKSN*XIALS 
PS5 = PSS 
PS6 = PSS 

GASKDN = GASC0N(F3) 

R = GASKON 
FM5C = 2.*FJ/R 
DLHriC = 2.*G*FJ 
RHOfaC = 144.*PS6/R 


THE INITIAL MS VALUE FOLLOWS 


FM5 = 1. 

KCT = 0 

11 IF(FM5 - 1. )6,7,7 

6 XIM = 0. 

GO TO 8 

7 XIM = 0.0075*«FM5*FM5-FM5) 

8 XIN = XIP + XIS + XIM 
ETAN = 1.- XIN 


CALL ADIABT(T4,P4,ETAN, F 3 , T SS , P S 5 , DH , E T A BR ) 
GAM = GAMMA(TSS,F3,ETABR) 

FM5 = SORT(FMdC*ABS( DH)/TS5/GAM) 


IF(KCT-l)9,i0,10 
9 KCT = 1 

FM5SV = FM5 
GO TO 11 

10 IF{ ABS(FM5-FMbSV)-.00 2)14, 14,12 
12 FM5SV = FM5 

KCT = KCT+1 
IFIKCT-lOOIll, 11,15 

15 WRITE(5,16) FM5,FM5SV,ETAN 

16 FORMAT! IH , 36H IN SUB TURBOP, FM5 CAN NOT BE FOUNO/IH , 

1 6H FM5 =E12.5,5X,8H FM5SV =E12.5,5X,7H ETAN =E12.5) 

BEGIN STATOR EXIT MACH NUMBER AXIAL COMPONENT FLOW 
FMZ5 CALCULATION 
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14 CP6= SPHTCP(TS5fF3,ETABR) 

C 

FMZ6C = GAM*G*R 
GAMl = GAM-1. 

GAMll = GAMI/GAM 
GAM2 = GAM/GAMI 
GAM3 = 2. /GAMl 
GAM4 = ( 1 .-GAM)/GAM 

GAMS = -GAM4 

VS = S9RT ( 2.*G*FJ*ABS ( OH ) ) 

PS = PS5*( l.+0.5*{ GAM-1. )*FM5*FM5)**GAM2 
RHOS = 144.*PS5/TS5/R 
FMRSC = SQRT(GAM*G*R*TSS ) 

BEGIN! FMZS COMPUTATION 

KMS = 0 

FMZS = FMS 

34 ALPHSC = FMZ5/FM5 
ALPH5= ARCOS (ALPHSC) 

VZ5= V5*ALPH5C 
VU5= V5*S IN( ALPHS ) 

AZ5 = 144.*n!4/VZS/RH05 
RU5= VU5-U 

BETAS = ATAI\I2(RUS,VZ5) 

RS = SORT ( VZS*VZS+RU5*RU5) 

FMRS= R5/FMK5C 

TO COMPUTE THE BUCKET EFFICIENCY, ETAB 
DELH =2.*BEIA5 
ZETABC = PIE J80*DELB/ 100. 

ZETABS = 0.025 + 0.034*ZETABC*ZETABC*SQRT(ZETABC) 
ZETAP = ZETABS/REXP2 
ZETAS = FKSe*ZETABS 
IF(FMR5-1.) 42,43,43 

42 ZETAM = 0.0 
GO TO 44 

43 ZETAM = CM* ( FMR5*FMR5 - FMR5) 

44 ZETAB = ZETAP + ZETAS + ZETAM 
ETAB = 1.- ZETAB 


CHANGES MADE 3-7-69 

IFITHARC - (360. - .01)) 60,61,61 

60 SPACE = TBH/CQOS/AR 
ARCL = PI E*DT*THARC/360. 

FKPA = 1. - SPACE/3. /ARCL 
R6 = RS*SQRT(ETA8)*FKPA 
GO TO 62 

61 R6 = R5*S0RT ( ETAB) 

62 DLHB = (R5*R5-R6*R6) /OLHBC 
TS6 = TS5+0LHB/CP6 

RH06 = RHQ6C/TS6 

VZ6 = 144.*W4/RH06/AZS 

FMZ6S = SQRT(FMZ6C*TS6) 

FMZ6 = VZ6/FMZ6S 
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FMR6 = R6/FMZ6S 

RU6 = SQRT(R6*R6-VZ6*VZ6 ) 

BETA6= ATAN2(RU6,VZ6) 

VU6 = RU6-U 
C 

ALPH6 = ATAN2CVU6.VZ6 ) 

C 

c 

c 

52 FM6 = FMZ6/C0S(ALPH6) 

P6 = PS6*(1.+0.5*(GAM-1. )*FM6*FM6)**GAM2 
PS = P6 

T6 = TS6*( l.+0.5*GAMl*FM6*FM6J 
FMSCAL = FMZ6 
C 

PS=PSS*(1 .+.5*(GAM-1 . )*FMSCAL**2)**GAM2 
IF(ABS(FMSCAL-FMS) -.001) 31,31,30 

30 KHS = <MS+1 
IF(KMS-IOO) 32,32,33 

32 FMZ5 = FMZ5*FMS/FMSCAL 
GO TO 34 

33 WRITE(6,35) FMS , FMSC A L , FMZ5 

35 FORMArUH , 90H IN SUB. TURB68,FMZ5 CAN NOT BE FOUND - PROGRAM 
ICONTINJES USING THE LAST VALUE 08TAINED/1H , 

2 6H FMS =E12.5,5X,7HFMSCAL=E12.5,5X,7H FMZ5 =E12.5) 

C 

C FMZ5 IS FOUND 

C 

31 TS = TS6*(1. + 0.5*(GAM - 1 . ) *FMZ6*FMZ6 ) 

TSS = TS6 

VS = FMS*SQRT(GAM*R*G*TSS) 

PS 10= PIO/I l. + 0.2*FM10*FM10) **3.5 
DELPS = 1.20*PAMB*(P10-PS10)/P10 
PSFAN = PAMB-DELPS 
SPR = PSFAN/PS5 
C 

PHI = 0.925*SORT( l.-SPR**AEXP) 

DELLIT = 0.00025*0FT 

WLEAK = PIE*DTR*DELLIT*S0RT(G/R»*PS5*0.615*PHI/S0RT(TS5) 

C 

41 ZLEAK = WLEAK/W4 

call A0IABT(T4,P4, 1.0,F3,TSP,PS,OHTTI ,ETABR) 

C 

C CHANGES MADE 3-7-69 

C 

IFITHARC - (360. - .01)) 63,64,64 

63 ETATL2 = RU5 + RU6 - .5*F14*U*U*COOS*SPACE/V5/ARCL/ALPH5C 
ETATL1=U*ETATL2/ABS( DHTTD/G/FJ 

GO TO 65 

64 ETATLl = U*(RU5 + RU6 ) /ABS ( DHTTI )/G/FJ 
C 

65 ETAT = ETATL1*( 1. -ZLEAK) 

FLAM = U/(RU5+RU6) 

VO = V5/SQRT(ETAN) 

40 FNUBR = U/VO 
C 
C 

RETURN 

END 

91 





SIBFTC TURB3 


SUBROUTINE FURB3 ( , WA , PA, P SS , DFT, FMS,ETAN, 

F3, FMIC, PIO, PAMB.FKS, SWIO.FNSML, 

ETABR,FK, 

ALPH5,ALPH6,ALPH7,ALPH8,AZ5,AZ7,AZ8,AZS. 

BET A5, BETA6, BET A7 , 8ETA8, CP6, CP8 , DELBTl , 
DELBT2,DHN1,DHN2,0LHR1,DLHB2,DT,ETAB1, 
ETAH2, ETAN1,ETAN2,ETAT,ETATL1, 
F<M,FKK1,FKM2,FM5,FM6,FM7,FMH,FMR5, 

FMR6, FNR7, FMR 8 , FM Z 5 , PM Z 5 S V , F M Z6 , FMZ 7 , 
FMZ7SV,FMZ8,GAM5,GAM6,GAM8,P5,P6, 

1 P7,P8,PHI,PS,PS5,PS6,PS7,PS6,PS1C5,R5,R6,R7, 

2 R8 , RH0‘5,RHn6,RH07,RHQ8,RU5,KU6,RU7, 

3 RUa.SHP.SPR, 

A T6,T8,TS,TSS,TSP, TS5 , T S6 , TS 7 , TS8 , U , 

5 V5,V7,VS, VU5, VU6,VU7,VU8,VZ8, VZ5, 

6 VZ6,VZ7, VZS,WIDEAL,WLEAK, XIBETl, 

7 XI 8ET2 ,ZLE AK, 

8 0MEfiA,CM,FKSB,FKSN,XIALS,TRM,XIAL6,CCM,FLLC,ZIL) 

C 

C TWO STAGE TURBINE WITH OGV 

C 

G= 32.17A 
FJ = 778. 

PIE = 3.1A15926 
PIElflO = 180. /PIE 
FKSl = FKS 
FKS2 = FKSl 
R = GASC0N(F3) 

RGG = 3*R 

ETA JCl = (1. +FKSS)* XIALP 
PS5 = FLLC*SORT( PA*P AMB ) 

RHOSC = lAA. *PS5/R 
DLH ’,C = 2.*G*FJ 
WAIA = 1AA.*WA 
AZSC = R*WA/PSS 
FM3 = FMS 
C 

C STARTING VALUES FOR FMZ5 AND FMZ7 

C 

C 

REA = 2A. *WA/UT/VISG { TA) 

REXP2 = IREA/1.5E6)**.2 
XIPA = XIALS/RCXP2 
XISA = FKSN*XIALS 
C 

KCT =0 
C 

C BEGIN FMS AND ETAN COMPUTATION 

C 

ETANl = ETAN 

15 call ADIABT(TA,PA,ETAN1,F3,TS5,PS5,DHN1,ETABR) 

DHNi = ABS(DHNl) 

V5 = S0RT(DLHBC*DHN1 ) 

GAMS = gamma { TS5, F3 , ETABR ) 

FM5SV = FMS 

FMS = V5/S0RT(GAM5*RGG*TS5) 

IF( ABS(FM5SV-FM5 1 -.002) 12.12,13 
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13 KCT = KCT+1 

IF(KCr - 100)14,14,16 

16 WRirE(5,17) FM5,FM5SV,ETAN,ETAN1 

17 F0RHAT( IH, 5UHIN TURB2 AT STATEMENT NO. 13, FM5 CAN NOT BE FOUND/ 

1 5H FM5=E12.5 ,5X,7H FM5SV=E 12 . 5 , 5X , 6H ETAN=E12.5, 

2 5X,8H ETA.Nl =E12.5) 

GO TO 12 

14 IF(FM‘5-1. )6,7,7 

6 XIM4 = 0. 

GO TO 8 

7 XIM4 = .0075*( FM5*FM5 - FM5) 

8 XI, N4 = XIP4 + XIS4+ XIM4 
ETANl = 1. - XIN4 

GO TO 15 

12 GAM51 = GAM5-1. 

GAM52 = GAM5/GAM51 

COMPUTE STARTING VALUES 

FOR FMZ5 AND FM27 COMPUTATIONS 

FMZ5 = FMS*CCM 
FMZ7 = FMZ5 
FMSCAL = 0 
KMS =0 

BEGINNING OF THE FMZ5 AND 
FMZ7 COMPUTATIONS 

42 FMZ505 = FMZ5/FM5 

ALPH5 = ARC0SIFMZ505) 

VZ5= V5 *FMZ5Q5 

VU5 = V5 * SORT! 1.-FMZ505*FMZ505) 

P5 = PS5 *(!.+. 5*GAM51*FM5*FM5)**GAM52 

RH05 = RH05C/TS5 

AZ5 = W414/VZ5/RH05 

RU5 = VU5 -U 

BETA5= ATAN2(RU5,VZ5 ) 

R5 = S3RT(VZ5«VZ5+RU5*RU5) 

FMR5 = R5/SQRT(GAM5*RGG*TS5) 

OEL8T1 = 2.*BETA5 

ZETABC = PIE180*DELST1/100, 

ZETABS = 0.025 + .034*ZETA8C*ZETABC*SQRTI ZETABC ) 

ZETAP = ZETABS/REXP2 
ZETAS = FKSB*ZETABS 
IF( FMR5 -1. ) 52,53,53 

52 ZETAM = 0.0 
GO TO 54 

53 ZETAM = CM*(FMR5*FMR5 - FMR5) 

54 ZETABl = ZETAP + ZETAS + ZETAM 
ETABl = 1. - ZETABl 

R6 = R5*SQRT ( ETABl ) 

DLHBl = (R5*R5-R6*R6) /DLHEC 
CP6 = SPHTCP(TS5,F3, ETABR) 

TS6 = TS5 + DLHB1/CP6 
PS6 =PS5 
AZ6 = AZ5 

RH06 = 144. *^PS6/R/TS6 
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VZ6 = W414/RH06/A26 

GAMd = GAMMA( TS6, F3, ETAOR) 

GAMol = GAM6-1. 

GAM62 = 0.5*GAM61 
GAM63= GAM6/GAM61 
FMZ6SQ=SQRT! GAM6*RGG*TS6 ) 

FMZf>= VZ6/FMZ6S0 
FMR6 = R6/FMZ6SQ 
RU6 = SQRT(R6*R6-VZ6*VZ6 ) 

BETA6 = ATAi'^2(RU6,VZ6) 

VU6 = RU6-U 

ALPH6 = ATAN2(VU6,VZ6) 

FM6 = FMZ6*SQRT ( VU6*VU6+VZ6*VZ6) /VZ6 
T6 = TS6*( 1.+GAM62*FR6*FM6) 

P6 = PS6*{ l.+GAM62»FM6*FM6)**GAM63 
PS7 =PSS 
FM7 = FMS 
ETAN2 = ETAN 

VISCT6 = VISC(T6) 

RE6 = 24.*W4/(VISCT6*DT*TRM) 

REXP6 = ( RE6/1.5E6)»» .2 
XIP6 = XIAL6/REXP6 
XIS6 = FKSN*XIAL6 

KCT =D 

26 CALL ADIA8T{T6,P6,ETAN2f F3tTS7,PS7,DHN2,ETABR) 

0HN2 = ABS(DHN2) 

GAM7 = GAMMA(TS7,F3, ETABR) 

V7 = S3RT ( DLHBC*DHN2 ) 

FM7SV = FM7 

FM7 = V7/S0RT(KGG<GAM7»TS7) 

IF( ABS(FM7-FM7SV) -.002) 21,21,22 

22 IF(KCr-lOO) 23,23,24 

23 1F(FM7-1.) 35,56,56 

55 XIM6 = 0.0 
GO TO 57 

56 XIM6 = .0075*(FM7*FM7 - FM7) 

57 XIN2 = XIP6+ XIS6 + XIM6 
ETAi-j2 = 1.- XIN2 

GO TO 26 

24 WRITE (6,25) FM7 , FM7 S V , E TAN2 , GAM7 

25 FORMATllH , 5HFM7 =E12.5,5X,8H FM7SV =E12.5,5X,8H ETAN2 =E12.5, 

1 5X,6H GAM7=E12.5) 

21 GAM71=3AM7-1. 

GAM?2=0.5*GAM7l 
FK = GA'17 

GAM73=3AM7/GAM71 
FMZ^C = FMZ7/FM7 
ALPH7 = ARC0S(FMZ7C) 

VZ7 = ^(7*FMZ7C 

VU7 = V7»S0RT( l.-FMZ7C*-»2) 

P7 = PS7* ( 1 . +GAM72*FM7*FM7 ) **GAM73 

RH07= 144. *PS7/R/TS7 

AZ7 = W414/VZ7/RH07 

02 = TRM*U 

RU7 = VU7-U2 

BETA7 = ATAN2(RU7, VZ7 ) 
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R7=SQRT (VZ7*VZ7+RU7*RU7) 

FMR7 = R7/SQRT (GAM7*RGG*TS7) 

DELBT2 = 2.*BETA7 
C 

ZETABC = PIE180*OELBT2/100. 

ZETABS = 0.025 + 0 .034»Z ET ABC*ZETABC*SQR T ( ZETABC ) 

ZETAP = ZETABS/REXP6 
ZETAS = FKSB*ZETABS 
IF(FMR7-1.) 58,59,59 

58 ZETAM = 0.0 
GO TO 50 

59 ZETAM = CM*(FMR7*FMR7 - FMR7) 

60 ZETAB2 = ZETAP + ZETAS + ZETAM 
ETAB2 = 1.- ZETAB2 

C 

R8 = R7*SQRT(ETAB2) 

DLHB2 = (R7*R7 -R8*R8 ) /DLHBC 
CP8 = SPHTCP(TS7,F3,ETABR) 

TS8 = TS7+DLHB2/CP8 
PS8 = PS7 

RH08 = 144. *PS8/R/TS8 
AZ8 = AZ7 

VZ8 = W414/RH08/ AZ8 
GAM8 = GAMMA(TS8,F3,ETABK) 

GAM81 = GAM8-1. 

GAM82 = 0.5*GAM81 
GAM83 = GAM8/GAM81 
GAM84 = GAM31/GAM8 
FMZ8C = SQRT(GAM8*RGG*TS8) 

FMZ8 = VZ8/FMZ8C 
FMR8 = R8/FMZ8C 
RU8= SQRT(R8*R8-VZ8*VZ8 ) 

BETA8 = ATAN2(RU8,VZ8) 

VU8=RU3-U2 

T8=TS8* ( l.+GAM82*FM8*FM8 ) 

ALPH8 = ATAN2(VU8,VZ8 ) 

FM8 = FMZ8»SQRT(VU8*VU8+VZ8*VZ8)/VZ8 
P8 = PS8*(1.+ GAM82*FM8*FM8)**GAM83 

DELALP = PIE180*ABS( ALPH8) 

DLALPC = DELALP/100. 

ZETAAS = 0.025 + 0.034*DLALPC*DLALPC*SQRT ( DLALPC ) 

ZETASS = ZETAAS*l 1./REXP6 + FKSB) 

OMEGA = ZETASS*( 1.0+0. 5*GAM8»FMS*FMS) 

PS = P3 - OMEGA*( P8-PS8) 

T8 = TS8*( l.+GAM82*FM8*FM8 ) 

TS = T8 

FMSCAL=S0RT(2./GAM81*( ( PS/PSS ) **GAM84 -1. ) ) 

DEBUG AZ8, VZ8,GAM8,GAM84,FMZ8C,FMZ8,RU8,BETA8, VU8,T8, ALPH8,FM8,P8, 
1 XIP,WSBAR,PS,T8,TS,FMSCAL 

IF(ABS(FMS -FMSCAL) -.001) 40,40,41 
41 KMS = XMS+1 

IF(KMS-IOO) 45,45,43 
45 FMZ7SV = FMZ7 
FMZ5S\/ = FMZ5 
FMZ5 = FMZ5*FMS/FMSCAL 
FMZ7 = FMZ5 
FMSCAL = FMS 
GO TO 42 
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43 WRITE (6,44) F MS , EMS C AL , FMZ 5 

44 FQRMAT(1H , ^8H I 'N TURH2 AT STATE. NO. 41, FMS CAN NOT BE FOUND/ 

1 IH ,6H FMS = Ei;.’.5, 5X.9H FMSCAL =t:12.3,5X, 

2 6H FMZ5=E12.5) 

END OF FMS COMPUTATION 

40 TSS^TS/ ( 1 .+bAMH2*FMS»FMS j 
VS= FMSftSORTICjAMBe-RGlT^fTSS) 

VZS=VS 

AZS= AZSC*TSS/VS 
DTR=XI BET2 
P10=P1 J 
FMl 1 = FM10 

PSlJ= P10/(l. + 0.2»FM10*FM10)*-*3.5 
PS10=PS10 

DEEPS = 1.20*PAMB»(P10-PS10)/P10 
PSFAN = PAMh-DFLPS 
SPR = PSFAN/PS5 
AEXP=XI 3ET1 

PHI = D.925»S(JRT( 1,-SPR»»AFXP) 

DECEIT = 0.00025*DFT 

WEEAK = P I E*DTR»DEEE I T*S0RT{G/R)*PS5*0.61S*PHI/S0RT( TS5) 

ZEEAK = WEEAK/W4 

CAEE ADIABT(T4,P4, 1.0,F3,TSP, PS , DH4S , E T AbR ) 

DH4S = ABS(DH4S) 

WIDtAE = W4*DH4S 

SHP = U* ( RU5 + RU6 + TRM*( 1 .-Z IE ) » (RU T+RU8 ) ) *W4/ ( G^^FJ ) 

ETATE1= SHP/WIDEAE 
ETAT = ETATE1*( i. -ZEEAK) 

46 RETURN 
END 


SIBFTC TURB4 


C 

C 

C 

C 

C 

C 

C 


2 

3 

4 

5 

6 

7 

8 
9 
1 
2 
3 
1 


SUBROUTINE TURH4(T4,K4,P4,PSS, DFT,DT, FMS.FTAN, 

F3, FM10,P10,PAMR, 

TS5 , V5,FM5.AEPH5, VZ5, VU5, P5,U,BETA5, 

R5 , FMR6,R6,ETAB,TS6, VZ6, FMZ6,FMR6,BETA6, 

RJ6 jOMEGA, ALPH6, P6, PS, T6,FMSCAE, WEEAK, 

PHI ,PS10,SPR,ZEEAK,ETATE1, ETAT, 

FE AM,FNUBR,ETA6R, 

CM, FKSB,FKSN,XIAES,CCM, 

AZ5 ,FM6,FMZ5,PSFAN,RE,RU5, TS,TSS, V0,VU6,XIM 
XI N,XI P,XIS,ZETAAS,ZETAB,ZETAM,ZETAP,ZETASS 
ZETAS, TSP.DHTTI , CP6 , DH , AEXP , DTR , 

TBH.THARC, CGOS ,AR, F14) 


MAIN REFERS TO THE CAEEING PROGRAM 
WFWA3P IN MAIN IS F3 IN TURBDd 
FM ll IN MAIN IS FMIO IN TURBOP 
Pll IN MAIN IS PIO IN TURBDP 
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PSCHK IN MAIN IS PS IN TUR80P 
PSIl IN MAIN IS PSIO IN TUR8DP 


G = 32.174 
FJ = 778. 

SIGS=1.2 

PIE = 3.1415926 

PIE18D = 180. /PIE 

RE = 24. ttW4/or/VISC( T4 ) 
REXP2 = (RE/1.5E6)*».2 
XIP = XIALS/REXP2 

XIS = fksn*xials 

PS5 = PSS 
PS6 = PS5 

GASKON = GASCON! F3) 

R = GASKON 
FM5C = 2.*FJ/R 
DLHWC = 2.*G»FJ 
RH06C = 144.»PS6/R 


IHE INITIAL M5 VALUE FOLLOWS 


FM5 = 1. 

KSD=1 

S0=. 

95 KMS=0 

101 FMZ5 = FMS*CCM 
KCT = 0 

11 IF(FM5 - 1. )6,7,7 

6 XIM = 0. 

GO TO 3 

7 XIM = D.0075*( FM5*FM5-FM5) 

8 XI N = XIP + XIS + XIM 
ETAN = 1.- XIN 


CALL ADIABT(T4,P4,ETAN, F 3, TS5 , PS5» OH, ETA8R ) 
GAM = 3AMMA(TS5,F3,ETABR) 

FM5 = S0RT(FM5C*A8S( DU)/TS5/GAM) 


IF(KCT-1)9, 10,10 
9 KCT = 1 

FM5SV = FM5 
GO TO 11 

10 IF ( ABS (FM5-FM5SV)-. 002)14, 14,12 
12 FM5SV = FM5 

KCT = KCT+1 
IF(KCT-IOO) 11,11,15 

15 WRITE(5,16) FM5,FM5SV,ETAN 

16 FORMAT! IH , 36H IN SUB TURBDP, FM5 CAN NOT BE FOUND/IH , 
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1 6H FM5 =E12.5,5X,8H FM5SV =E12.5,5X,7H ETAN =E12.5I 

BEGIN STATOR EXIT MACH NUMBER AXIAL COMPONENT FLOW 
FM25 CALCULATION 

lA CP6= SPHTCPI TS5,F3,ETABR ) 

FMZ5C = GAM»G*R 
GAMl = GAM-1. 

GAMll = GAMI/GAM 
GAM2 = GAM/GAMI 
GAM3 = 2. /GAMl 
GAMA = (l.-GAM)/GAM 
GAMS = -GAMA 

V5 = S3RT(?.*G*FJ*A3S(DH) ) 

P5 = PS5*(1.+0.5*(GAM-1. )*FM5*FM5)**GAM2 
RH05 = 1AA.»PS5/TS5/R 
FMRiC = SORT ( GAM*G*R»TS5 ) 

BEGIN FMZ5 COMPUTATION 

3A ALPH5C = FMZ5/FM5 
ALPH5= ARCGS ( ALPH5C) 

VZ5= V5*ALPH5C 
VU5= V5*S 1 Nl ALPH5) 

AZ5 = 1 AA. *WA/VZ5/RH05 
RU5= VU5-U 

BETAS = ATAN2(RU5,VZ5 ) 

R5 = SORT ( VZ5*VZ5+RU5*RU5) 

FMRS= R5/FMR5C 

TO COMPUTE THE BUCKET EFFICIENCY, ETAB 
DELB =2.*BETA5 
ZETABC = P IEl8n*nELB/ 100 . 

ZETABS = 0.025 + 0.03A*ZFTABC*ZETABC*SQRT { ZETABC ) 

ZETAP = ZETABS/REXP2 
ZETAS = FKSB*ZETABS 
IF(FMR5-1.) A2,A3,A3 
A2 ZETAM = 0.0 
GO TO AA 

A3 ZETAM = CM*(FMR5*FMR5 - FMR5) 

AA ZETAB = ZETAP + ZETAS + ZETAM 
ETAB = 1.- ZETAB 


IFITHARC - (360. - .01)) 60,61,61 

60 SPACE = TBH/COOS/AR 
ARCL = PI E*0T*THARC/360. 

FKPA = 1. - SPACE/3. /ARCL 
R6 = R5*S0RT( ETAB)*FKPA 
GO TO 52 

61 R6 = R5*SQRT ( ETAB) 

62 DLHB = (R5*R5-R6*R6) /DLHBC 
TS6 = TS5+DLHB/CP6 

RH06 = RH06C/TS6 

VZ6 = 1AA.»WA/RH06/AZ5 

FMZ6S = SORT( FMZ6C*TS6) 
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FMZt) = VZ6/FMZ6S 

FMR6 = R6/FMZ6S 

RU6 = iQRT(R6*R6-VZ6*VZ6) 

BET<\6= ATAN2(RU6,VZ6) 

VU6 = RU6-U 

ALPH6 = ATAN2C VU6.VZ6) 

51 OMEGA=.008^fSIGS*EXP( 3 . 7* SD ) /COS ( ALPH6 ) 

52 FM6 = FMZ6/C0S(ALPH6 ) 

P6 = PS6*{ l.+0.5*(GAM-l. )*FM6*FM6)**GAM2 
PS = P5-0MEGA*1P6-PS6) 

T6 = TS6*( l.+0.5*GAMl»FM6*FM6) 

V6=S0RT(VU6**2 + VZ6**2 ) 

VS=V6» ( l.-SD+SIN( ALPH6)/2./SIGS» 

IF(KMS.EQ.O) rSS=TS6 
FMSCAL=VS/SQRT (GAM*32 . 17*R*TSS) 

TS = T6 

TSS=TS/ ( l.+GAMl/2.*FMS**2) 

PSSN=PS/{ l.+(GAM-l.) /2.*FMS**2)**GAM2 

PS5=PS5*14.7/PSSN 

PS6=PS5 

IF(ABS{PSSN-14.7).LT..01) GO TO 31 
PSS=PSSN 

30 KMS = KMS+1 

IFlKMS-100) 32»32,33 

32 GO TO 101 

33 WRITE(6,35) FMS, FMSCAL , FMZ5 

35 F0R“AT(1H , 90H IN SUB. TURB68.FMZ5 CAN NOT BE FOUND - PROGRAM 
ICQNTINUES USING THE LAST VALUE OBTAINED/IH , 

2 6H FMS =E12.5»5X,7HFMSCAL=E12.5,5X,7H FMZ5 =E12.5) 

FMZ5 IS FOUND 

31 GO TO ( 80,81,90) ,KSD 

80 FMS1=FMSCAL 
SD1=SD 
SD=.6 

KS0 = 2 
GO TO 95 

81 FMS2=FMSCAL 
SD2=SD 

KSD = 3 

SD=( SD2-SD1 )/ ( FMS2-FMS1 ) *( FMS-FMSl ) +SD1 
GO TO 95 
90 CONTINUE 

DEBUG FMS, FMSCAL, SD,PSS,PS5, VS 
TS = T6 

PS 10= PIO/ ( 1.+0.2*FM10*FM10)**3.5 
DELPS = 1.20*PAMB*(P10-PS10) /PIO 
PSFAN = PAMB-DELPS 
SPR = PSFAN/PS5 
IFISPR.GT.l. ) SPR=1. 

PHI = D.925*S0RT( l.-SPR**AEXP) 

DELLIT = 0.00025*DFT 

WLEAK = PICttDTR*DELLIT*SQRTlG/R»*PS5*0.615*PHI/S0RT{TS5) 

41 ZLEAK = WLEAK/tJ4 

CALL ADIABT(T4,P4, 1.0,F3,TSP,PS,OHTTI,ETABR) 

IFITHARC - (360. - .01)) 63,64,64 
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63 FTAFL2 = RDo + RU6 - . 5 * F 1 4* U*U*COnS* S P ACE / V 5 / ARC L / AL PH5C 
ErA[Ll=U*CTATL2/ABS(DHTTI)/G/FJ 

GO ro 55 

64 ETATLl = U»(RU5 + RU6 ) / ABS ( DHTT I ) /G/ F J 

65 ETAr = ETAFLIM l.-ZLEAK) 

FLAV = u/ ( RU5+KU6 ) 

VO = vs/soRrcErAN) 

40 FMUiiR = G/VU 


RETUR'J 

END 


I8FTC HAPir 


SUBROjriNE HALPLTlJFMS.FMSJt AL PH5 J , DTT J , FMR 5 J , TBHJ,W1CMP, 

1 VLJJJ.FPR, n,T4,UFT,Pl, ITURB) 

COM'-inM/CLPLUT/XPFN, YPFN, NX ,.NY, I PFiV , XL ABEL (lO).YLABEL(lO) 

DATA BCDX/6HSTATUR,5H EXIT , 6HMACH N,6H0., MS/, 

1 5CDY/6HPARA. ,6HRATI0,,6H BASE ,6HPT. MS,6H = 0.5/, 

2 L A0EL/6HALPH5 ,6HDTT ,6HMR5 ,6HT8H ,6HW1 , 

3 5HVL / 

DATA MAM1/6HFPR = /,NAM2/6H T1 = /,NAM3/6H T4 = /,NAM4/6HUFT = /, 
1 'lAMS/BH PI = /, NAM7/6HI TURB = / 

DIMENSION ALPHbJI 100 ) , DT TJ ( 1 00 ) , FMR5 J ( 1 00 ) , TBH J ( 100 ) , 

1 VLJJJ(100),W1CMP(100),FMSJ(100) 

D I MENS ION Y J ( 100) ,3CDX( 4) ,BCDY{ 5 ) ,LABEL(6) 


CALL lOENT 

CALL PLOKO. ,.5,-3) 

FITURB = ITURB 
OX = .05 

XMI ) = FMSJ ( 1 ) - DX 

FLONGX = ABS(FMSJ( JFMS) - FMSJ(1))/DX + 1. 
YMIN = .8 

OY = . 1 
FLONGY = 8. 

YH I OH = . 1 


TO DRAW THE X AXIS 

CALL AXI S ( 0. ,0. ,BCDX , -24 , FLONGX , 0 . 0, XMIN,DX, lu. ) 
TO DRAW Y AXIS 

CALL AXIS (0. ,0. ,8CDY, 30, FL0NGY,90.,YMIN,DY, 10. ) 
TO DRAW GRID LINES 
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63 ETATL2 = RU5 + RU6 - . 5*F14*U*U*C0nS*SPACE/V5/ARCL/ALPH5C 
ETAtH=U*ETATL2/ABS( DHTTI) /G/FJ 
GO ro 65 

6A ETATH = U*(RU5 + R J6 ) /AHS ( DHTTI ) /G/ FJ 
C 

65 ETAT = ETATL1*( I.-ZLEAK) 
flan; = U/(RU5 + RU6) 

VO = V5/S0Rr(ETAN) 

40 FNUBR = U/VU 
C 
C 

RETURN 

END 


SI8FTC HAPLT 
C 

SUBR(3JT INE HALPL T ( JF M S , FMS J , ALPH5 J , DTT J , FMR 5 J , TBHJ ,W1CMP, 

1 VLJJ J.FPR, n,T4.UFT,Pl , I TURR) 

C 

COMMON/CLPLUT/XPEN, YPEN, NX.NY, IPEN, XL ABEL ( 10 ) , YL ABEL ( 10 ) 

C 

DATS BDOX/6HSTATOR.6H EXIT , 6HMACH N),6H0., MS/. 

1 3CDY/6HPARA. ,6HRATI0,,6H BASE ,6HPT. MS.6H = 0.5/, 

2 LA8EL/6HALPH5 ,6HDTT .6HMR5 ,6HTBH ,6HW1 , 

3 6HVL / 

C 

data NAM1/6HFPR = /.NAM2/6H T1 = /,NAM3/6H T4 = /,NAM4/6HUFT = /, 
1 NAM5/6H PI = /, NAM7/6HI TURB=/ 

C 

DIMENSION ALPH5J ( 103 ) .DTTJI 100) ,FMR5J ( 100) , TBHJ ( 100) , 

1 VLJ JJ ( 100 ) , WICMPI 100) , FMS J (100) 

DIMENSION YJ(100),BCDX{4) , BODY ( 5 ) , LABEL ( 6 ) 

C 

CALL IDENT 

CALL PL0T(0.,.5,-3) 

C 

FITURB = ITURB 

DX = .35 

XMIN = FMSJ ( 1 ) - DX 

FLONGX = ABS ( FMSJ ( JFMS ) - FMSJ(1))/DX + 1. 

YMIN = .8 

DY = . 1 

FLONGY = 8. 

YHIGH = .1 


TO DRAW THE X AXIS 

CALL AXIS(0.,0., BCDX,-24,FLONGX,0.0,XMIN,DX, 10. ) 
TO DRAW Y AXIS 

CALL AXIS (0. ,0. , BCOY, 30, FLONGY.90., YMIN.DY, 10. ) 
TO DRAW GRID LINES 
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noooo o nnoooo 


FLCjXX = FLONGX + .l*FLnNGX 
JLOrJGX = IFIX(FLGXX) 

C 

X = 0. 

DO 60 J=l,JLnNGX 
X = X + 1. 

CALL PL0T(X,0.,3) 

60 CALL PL0T(X,FLQNGY,2 ) 

FLGYY = FLONGY + .1*FL0NGY 
JLOMGY = IFIX(FLGYY) 

Y = 0. 

C 

DO 62 J=l,JLO'JGY 

Y = Y + 1. 

CALL PLOT (0. , Y , 3 ) 

62 CALL PLOT ( FLONGX, Y,2 ) 

GRID LINES ARE DRAWN 


TO DRAW THE FIVE LINE PLOTS 

DO 300 JJ=1,5 

GO TO ( 31 , 32 ,33. 34,35 ) , J J 

31 DO 36 J=1,JEMS 

36 YJ(J) = ALP^^5J(J) 

GO TO 30 

32 DO 37 J=1,JFMS 

37 YJU) = OTTJ(J) 

GO TO 30 

33 DO 38 J=1,JFMS 

38 YJU) = FMR5JU) 

GO TO 30 

34 DO 39 J=1,JFMS 

39 YJU) = TBHJU) 

GO TO 30 

35 DO 40 J=1,JFMS 

40 YJU) = W ICMP I J ) 

GO TO 30 

30 DUMY = DUMY 

300 CALL LINE(FMSJ,YJ,JFMS,1,0,1.XM1N,0X,YMIN,DY) 
LINE PLOTS ARE DRAWN 
TO LABEL CURVES 

X = (FMSJ(l) - XMIN)/DX/2. 

DO 45 JJ=1,5 

GO TO (46,47,48,49,50) ,JJ 

46 Y = ALPH5J(i) 

GO TO 51 

47 Y = DTTJ( 1 ) 

GO ro 51 

48 Y = FMR5J ( 1 ) 

GO TO 51 

49 Y = TBHJ( 1 ) 

GO ro 51 
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50 Y = WlOMPdJ 
GO TO 51 

12 'f = VLJJJ ( 1 ) 

51 DUMY = DUMY 

Y = (Y - YMrN)/DY + .1 

45 CALL SYMBOL(X,Y,YHIGH,LABEL( JJ»,0.0,6) 

C 

C TO DRAW LEGEND PLOT 

C TO MOVE TO NEW PLOT TO WRITE LEGEND 


C 

CALL PLOT ( 10. ,0. ,-3) 

CALL SYMBOL (0.0, 9. 0, . 1 5 , NAMl , 0. 0, 6 ) 

CALL NJMBER(1.05,9.0,.15,FPR,0.0,3> 

CALL SYMBOL (0.0, 8. 75, . 1 5 ,NAM2 , 0. 0 , 6 ) 
CALL NJMBER(1.05,8.75,.15,T1,0.0,3) 

CALL SYMBOL (0.0, 8. 50,. 15, NAM 3, 0.0, 6) 
CALL NJMBER(1.05,8.50,.15,T4,0.0,3) 

CALL SYMBOL(0. 0,8. 25, . 15 ,NAM4,0.0, 6) 
CALL N J MB ER( 1.05, 8. 25,. 15, OFT, 0.0, 3) 
CALL SYMBOL (0.0, 8.00, . 1 5 ,NAM5, 0 . 0 , 6 ) 
CALL NJMBER(1.05,8.00,.15,P1,0.0,6) 

CALL SYMB0L(0.0,7.75,.15,NAM7,0.0,6) 
CALL NJMBER( 1 . 05 , 7 . 7 5 , . 1 5, F I TURB , 0 . 0 , 6 ) 
C 

CALL PLOT( 10. ,0. ,-3) 

C 

RETURN 

END 


$I8FTC SFaLSE 

SUBROUTINE F AL SE ( X , Y , CO VTST, J , KCOUNT ) 

C 

C 

C 

C TO FIND X BY THE METHOD OF FALSE POSITION SUCH THAT 

C FIX) = Y =0 

C 

CCC CCC REFERS TO WHAT MUST BE DONE IN CALLING PROG. 

C 

CCC COVTST IS THE CONVERGENCE CRITERION AND IS 

COMPUTED IN THE CALLING PROGRAM 

THE CONDITION ABSIY) LESS THAN COVTST 
MUST BE MET 

Y IS ALWAYS COMPUTED IN THF CALLING PROG. 

TO BE USED IN FALSE 

J IS AN INDICATOR 

J MUST BE INITIALLY SET EQUAL TO -1 IN THE 
CALLING PROGRAM 

J = 0,1,2 ARE SET IN FALSE 
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C J = -1 MEANS X = XFIRST IS BEING USED IN FALSE 

C J = 0 MEANS X = XSECONU IS BEING USED IN FALSE 

C J = 1 MEANS X IS TO BE COMPUTED IN FALSE AND 

C RETURNED TO MAIN 

C J = 2 IS A FLAG SENT TO THE CALLING PROGRAM SAYING 

C WE HAVE CONVERGENCE 

C 


ccc 

IN 

THE 

CALLIN3 

PROGRAM, TWO READ IN X INITIAL GUESS ARE 

c 


NEEDED, XFIRST AND XSECOND. THESE TWO GUESSES 

c 

p 



ARE USED ONLY IN THE CALLING PROG. 

ccc 

IN 

THE 

CALLING 

PROG. 

ccc 


IF 

J = - 

1 , SET X = XFIRST 

ccc 


IF 

J = 0 

, SET X = XSECOND 

ccc 


IF 

J = 1 

, USE THE X FROM FALSE 

ccc 


IF 

J = 2 

, STOP. X IS FOUND 


c 

C KCOUNT IS ONLY USED IN FALSE. IT IS IN THE CALLING VECTOR 

C TO TELL THE CALLING PROG. HOW MANY INTERATIONS WERE 

C NEEDED 

C 
C 

c 

IF(J)10,15,20 

C 

10 XI = X 
Y1 - Y 
J = 0 

KCOUNT = 0 
GO TO t^0 
C 

15 X2 = X 
Y2 = Y 
J = 1 

GO TO 30 
C 
C 

20 IF(ABS(Y) - COVTSn 21,21,22 

21 J = 2 

GO TO ^0 
C 

22 KCOUNT = KCOUNT + 1 
IFIKCDUNT- 200)18,18,23 

23 WRITE(5,24) X,Y,C0VTST 
J = 2 

GO TO 40 
C 

c 

18 IF(ABStYl) -ABS(Y2)) 16,16, 17 

16 Y2 = Y 
X2 = X 
GO TO 30 

C 

17 Y1 = Y 
XI = X 

C 

C 

30 X = X2 - Y2*(X2-X1)/ (Y2-Y1) 

GO TO 40 
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c 

c 

24 F0RMAT(1H ,50H IN THE FALSE POSITION ROUTINE, KCOUNT IS EXCEEDED/ 
1 4H X =E12.5,5X,4H Y =E12.5,5X,9H CUVTST =E12.5) 

C 

40 RETURN 
END 


tIBFTC TFROS 

FUNCTIOM TFRDMS(S,F, ETABR) 

C 

C TO SOLVE FOR TA BY THE FALSE POSITION METHOD 

C TO GET A FIRST APPROXIMATION FOR TA 

C 

SF = F*ETABR/ ( 1 .+F ) 

FTATST =.0001*A8S ( S) 

TA = P0LY(S,8) 

PHIMS = P0LY(TA,2) + SF* POL Y ( TA, 5 ) 

PHIA = P0LY(TA,2)*(S/PHIMS)*»1.2 
TA1= POLYIPHIA.y) 

C 

C bEGiN FALSE POSITION METHOD 

C 

TA2 = TAl + .02*TA1 
JFALSE = -1 
7 IF(JFALSE) 11,12,9 

11 TA = TAl 
GO TO 9 

12 TA = TA2 

9 FTA = P0LY(TA,2) + S F * POLY ( T A, 5 > - S 
CALL FALSE! TA, FTA, FT ATST,JFALSE,K) 

IF (JFALSE-1)7,9,30 
30 TFROMS = TA 
RETURN 
END 


SIBFTC TFROH 

FUNCTION TFROMH(HM,F , ETABR » 

C 

C TO SOLVE FOR TA BY THE FALSE POSITION METHOD 

C TO GET A FIRST APPROXIMATION FOR TA 

C 

SF - F*ETABR/ ( l.+F ) 

FTAIST =.0001*ABS(HM ) 

HA= HM 

TA= P0LY(HA,7) 

PSIH = POLYI TA,4) 

HMSTR = HA + SF*PSIH 
HA = HA*(HM/HMSTR)**1.2 
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TAl = PDLY«HA,7) 

BEGIN FALSE POSITION METHOD 

TA2= TA1+.02*TA1 
JFALSE = -1 
7 IFUFALSE) 11,12,9 

11 TA = TAl 
GO TO 9 

12 TA = TA2 

9 FTA = POLY(TA,l) + S F * POLY ( T A , 4 ) - HM 
CALL FALSE (TA, FTA, FT ATST,JFALSE,K) 
IF(JFALSE- 1)7,9,30 
30 TFRGMH = TA 
RETURN 
END 


$IBFTC POLYA 

FUNCTION P0LY(X,M) 

DIMENSION ICC(B) , IORD( 3) ,C( 300) 

DATA( I CC( I ), 1= I , 8 )/ 133, 146, 161, 174, 187,200,209,226/ 
DATA(I3RD(I),I=1,8)/4,5,A,4,4,2,6,6/ 

DATA(C( I ), 1 = 133,242) /lOO., 107. 26 16,2. 4975E-1,- 2.26 3 8E-5, 
11.96075E-8,-3.675E-12,2000., 116.484, . 2096 1 0 , 2 . 354 7 1 3E - 5 , - 3 . 338388E 
2-9, 1.8 343 E- 13, 6000., 100. ,1.14506, 1 . 4941 3E- 3 , - 1 . 7 98 3 1 E -6 , I . 394 76E -9 
3,-5.8514E-13,1.0156E-16, 1400., 1 . 429 38 , 4 . 43 7 1 22E-4 , - 1 . 489 1 8E-7 , 
43.4627F-1 1 ,-4.513E-l 5. 2.47>^-19, 6000. , 100.,.252 32,-5.44152E-5, 
57.0682E-8 ,-2.017 lE-1 1 , -5 . I E- 16 , 1400 . , . 1 86 1 , 8 . 0 148E-5 , -2 . 3278E-8 , 
63. 4 1635 E- 12,-1 .989E- 16,6000. , 100. , 95 7 . 028 , -1 . 397247E-2 , 2 . 72 8083E -4 
7,-7. 87 4997E-8, 1 . 1 31 1 8 4 E- 1 1 , 2 000 . , 989 . 0299 , 2 . 42 1 0 5 8 E-2 , 1 . 7 1 9862E-4 , 

8- 2.13151E-8,1.027056E-12,6000., 100 . , - , 224336 , 4 . 84768E-4, -7 . 184546E 

9- 8,-4.315008E-12,3.134744E-15,2000.,-4.089761,5.504657E-3,-2.44669 
A9E-6,4.89935lh-10,-3.617383E-14,6000.,100.,3.986078E-2,3.562965E-4 
6,-6.51282 IE- 8, 2000.,. 206 39 79,1 . 7 64343E-4 , - 1 . 68 3 1 3 76-8 , 6000 . , 100. , 
C-62. 685 51 6,-6.922 1196F-1 , 1.9220060E-2,- 1.1 3456llE-5,-9. 8360870E-8, 
D2. 24 58 1 96 E-10, -1.-^50 96 78 £-13, 614., 1946. 1462,-8.5195086,2.6225174E- 
E2.-2.5543759E-5. 1 . 1 47 9667E-8 ,-6 . 4275539E-1 3 , -6 . 082 041E- 1 6 , 2000. , 
FI. 4, 41 5. 37 102, 58.8657 53, -312 .50006, 16. 52 8755 ,- 29 . 769562 , 5 . 6792908 , 
G53. 206 712, 1.9338 4, 55 3. 84976, 12 3. 43567, - 1 74 . 32 1 29 , 22 . 1 84326 , -1 04. 70 
H82 3, -43. 3 2 23 74, 8 2. 93 8 049 ,2.3/ 

J=1 ABS ( M) 

K=I ecu ) 

10 L=I0RD(J)+2 

11 IF(X-C(K) ) 12,13,13 

12 K=K-L 

GG TO 11 

13 K=K+L 

IF(X-C(K)) 15,15,13 

15 IF(M) 20,16,16 

16 L=L-1 
P0LY=C1 K-1 ) 
no 17 N=2,L 
KN = <-N 
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noon 


17 P0LY=P3LY*X+C(KN) 

GO TO 22 

20 L=L-2 

P0LY=FL0AT(L)*C(K-1) 

DO 21 M=2,L 
KN=K~N 
LN=L-N4 + 1 

21 P0LY=P3LY*X+FL0AT(LN)*C(KN) 

22 RETURN 
END 


$IBFTC FM3S 
C 

FUNCTION FMM(F3,G3,S) 


TO SOLVE AN FQUA. IN FM3 BY NEWTON'S 

J = 0 

S = 2.0*(G3-1. )/(G3+l. ) 

SI = l./S 

FI = ! F3/SQRT(G*G3)) **S 
F2 = 3. 5*(G3-1. )*F1 
X = JF1+.25*F2}**S1 
C 

10 XS = X**S 

F = FI +F2*X*X -XS 
OF = 2.*F2*X - S»XS/X 

XSV = X > 

11 X = X -F/DF 
C 

IF( .331*A3S(XS) - ABS ( F ) ) 13, 13, 12 

12 IF( .03305*ABS( X + XSV) - ABS(X-XSV)) 13,13,15 
C 

13 J = J+1 
IF(J-130) 10,lu,14 

14 WRITE(6,17) X,XSV,F 
C 

15 FMM = X 
C 

16 RETURN 

17 FORMATdH ,6H FM3 =E12.5,5X,8H FM35V =E12.5,5X,4H 
END 


$IBFTC ENTR3 

FUNCTION ENTRUP(T,F,EIABR) 
PHIA = P0LY(T,2) 

SF = F*ETABR/ ( l. + F ) 


METHOD 


F =tl2. 



PSI = P0LY(T,5) 

PHIM = PHIA + SF*PSI 
ENTROP = PHIM 
7 RETURN 
END 


$IBFTC ENTHA 

FUNCTION ENTHAL(T,ETABR,F) 
C 
C 

HA = PDLY(T,1) 

PSIH = POLY(T,4) 

4 SF = F*ETABR/( l.+F) 

HF = SF*PSIH 

HM = HA+riF 

ENTriAL = HM 
6 RETURN 
END 


$IBFTC FULAR 
C 

FUNCTION FULAIR(T1,T2,F1,ETARR) 

C 

KCT = 0 
F2T = .02 
3 KCT= KCT+1 
F3 = FI + F2T 
H2 = ENTHALI T2,F3,ETABR ) 

HI = ENTHAL(T1,F3,ETARR ) 

F2 = (1. +F1)*(H2 -HI )/( 18630. *ETABR +H1-H2) 

6 IF(ABS(F2 - F2T) -l.OE-6) 9,9,7 

7 F2T = F2 

8 IF(KCT -100) 3,3,11 

11 WRITE(6,12) F2,F2T,H1,H2 

12 FORMATdH ,36H IN SUB. FULAIR, F2 CAN NOT BE FOUND/IH , 

1 4HF2 =E12.5,5X,AHF2T=E12.5,5X,4HH1 =E 1 2 . 5 , 5X , 4HH2 =E12.5) 

9 FULAIR = F2 
10 RETURN 

END 


$IBFTC ADIAB 


SUBROUTINE ADIA8T(T1,P1,FTAT,F.T2,P2,DH,ETABR) 
C 
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1 SI = EMTRnP(Tl,F,ETABR) 

R = SASCON(F) 

S2 = SI + ( R/778. 161 )*AL0G{ P2/P1 ) 
T2J = TFRnMS{S2tF,ETABR ) 

HI = FNTHAL(T1,F,ETABR) 

2 H2J = ENTHAK r2J,F,ETABR) 
IF(P2-P1) 4,5,3 

3 OH = ( H2J-H1 ) /ETAT 
GO in 6 

5 OH =0.3 
GO TO 6 

4 OH = ETAT* ( H2J-H1 ) 

6 H2 = Hl+OH 

T2 - TFRnMH(H2,F,FTABR) 

10 RETUR'l 
END 


tIBFTC IDENT 

SUBBOJTINE IDENT 

CnM BONI/XCPIDX/NBLONO, I RUNO , C P I D ( 8 ) 

DAT'\ R JNMES,BLKRES/6HRUN N0,6HbLK NO/ 
CALL PL0T(O.,0.,-3) 

FBLOND=NBLnNO 

CALL SyMB0L(0.0,1.0,.2,CPl0(l),90.,24) 
CALL SYMBOL (0.5, 1.0, . 2 , R UNME S , 90 . , 6 ) 
CALL symbol (.5, 2. 5, .2, I RUNG, 90., 6) 

CALL SYM80L( .5,4., .2,8LKMES,90.,6) 

CALI. NUMBER(.5,6.,.2, FRLC;NO, 90. ,-l ) 
CALL PLOT ( 5. ,0. ,-3 ) 

RETURN 

END 


tIBFTC GA^MAF 

FUNCTION GAMMA(T,F,ETABR ) 

CP = SPHTCP ( T, F,ETABR ) 

GAMMA = CP/ (CP-GASCON( F I/778.16) 
10 RETURN 
END 


UBFTC SPHTC 
C 

EUNCTION SPHTCP(T,F, ETABR) 
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c 


CPA = P0LY(T,3) 

SF = F*ETABR/ ( l.+F ) 
PSICP = P0LY(T,6) 

CPM = CPA + SF*PSICP 
SPHTCP = CPM 
10 RETURN 
END 


tIBFTC GASCN 

FUNCTION GASCGN(F) 

R = (53.349 + 55.088*F»/(1.+F» 
GASCON = R 
15 RETURN 
END 


SIBFTC VISCOS 

FUNCTION VISC(T) 

DIMENSION A(50) ,V(50 ) 

DATA ( V( J 1, J= 1,50 / 73. 8, 136. 0, 1 85. 2, 227. 2, 2 64. 7, 2 99. 2, 331. 3t 
A361 .4, 389.8,417. 1 , 443 . 5 , 469 . 5, 495. 1 , 5 19 . 7 , 543 . 6 , 56 7 . 0 , 589 . 8 , 6 1 2 . 1 
8633.9,555.3,676.3,697.0, 717.3,737.3,757.0,776.5,795.6,814.5,833.2 
C851.6, 869.8,887.8, 90 5. 6, 92 3.2,940.6,957.9, 974.9,991.8,1008.6, 

0102 5.2,1041.6, 1058.0, 107 4. 1, 1090.2, 1106.1, 1121.9,1137.5,1153.1, 
61168.5,1183.8/ 

TKELVN= T*5./9. 

A( 1 ) = 100. 

DO /O N=2,50 
A(N)= A(N-1)+100. 

70 CONTINUE 

IF ( TKELVN.LE.AI 1 ) ) GO TO 71 
IF { KELVN.GF .A( 50 I ) GO TO 72 
K= 2 

75 IF (TKELVN.LT.A(K) ) GO TO 74 
K= K + 1 

GO TO 75 

74 TOP= AIK)-TKELVN 
DELTT= TOP/100. 

DV= V(K)-V(K-1) 

DIFF= DELTT-*OV 
VISM= V(K)-DIFF 
VISC=VISM/14882000. 

GO TO 80 

71 VISC=\/( 1 ) /14882000. 

WRI TE(6,76) 

76 FQRMAT(32H TKELVN IS LESS THAN 100 DEGREES) 

GO TO SO 

72 VISC=\/( 50) /14882000. 

WRITE! 5,77) 

77 F0RMAT(36H TKELVN IS GREATER THAN 5000 DEGREES) 

80 RETURN 

END 



noon 


SIBFTC CGMP 

listing gf alter deck 


♦ ALTER 143»14<h 
♦ENDAL 


C 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE COMP( I OUT, I WGHT , FMS , ALPH5 , DTT , FMR 5 , TBH , W1 , FN, VL ) 
COMMON/DATl/ ADAS, OMEGAS, 

1 DTLOSS, ETABR,ETAF,ETAN,ETAT,FKV, 

2 FM11,FM2,FM2B,FM3S, FNB , FNS , FPR , H2P , HHV , PI , 

3 P1413, PAMB, POPS, PP2B,R,RHRT, RING, TAMB, 

4 T1,T4,UFT, W10BR,WFWA1 ,WLRAR,CT, 

5 CM, FKS8, FKSN ,X I ALS,CCM, AEXP, 

6 THARC.COQS, AR.F14 

COMMON/DAT2/ ITURB, IETAT,IWRITE,PIE,FK,TRM,XIAL6, 

1 FLLC,W1S2.P20P1,W1S1,FN3 

LIFT FAN DESIGN POINT PROGRAM JULY 1968 HALLER-AUER 


IF lETAT = 1, DO NOT INTERATE ON ETAT 
IF lETAT = 0, INTERATE ON ETAT 

IF ITURB =1, USE SINGLE STAGE TURBl ROUTINE WITH OGV 
IF ITURB = 2, USE 

SINGLE STAGE TRUBINE ROUTINE - STATORLESS 
IF ITURB = 3, USE TWO STAGE TRUBINE WITH OGV 
IF ITURB = 3, THE ADDITIONAL INPUT FK IS NEEDED 


IF IWIRTE = 2, THEN COMPUTE WITH THREE VALUES OF W1 
THE THIRD VALUE OF W1 (W1S3) IS COMPUTED. THE FIRST AND 
SECOND VALUES ARE READ IN 

THE THIRD VALUE OF FN (FN3) IS READ IN. THE FIRST AND 
SECOND VALUES ARE COMPUTED. 

IF IGUT = 0, DO NOT WRITE OUTPUT 
IF lOUT = 1, WRITE OUTPUT 


IF(IOUT) 81,81,73 
73 WRITE(6,106) FMS , P20P I , W 1 , FN3 
81 PIE4= 4. /PIE 
G = 32.16 
GAM = 1.4 
XIBET1=AEXP 


TO CONVERT ALPH7 TO RADIANS 


GAMCl = GAM-1. 

PDPS = 1. - .1*(1. - 1. /ADAS/ADAS)*! 1. 

1 GAMC1*FMS*FMS)**(GAM/GAMC1) ) 

GAMC2 = GAMCl/GAM 
GAMC3 = 2. /GAMCl 


l./( 1. + 0. 5^ 
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ZLEAK =0. 

W2 = W1 
T3 = 

TIO = TAMB 
Til = TIO 
P2 = P1*P20P1 
T2 = T1 ■- DTLOSS 
R2 = GASC0N(WFWA1 ) 

CP2 = SPHTCP (T2»WFWA1,BTABR) 

GAM2 = CP2/ (CP2-R2/778.16) 

GAM2C1 = GAM2 + 1. 

GAM2C2 = GAM2 -1. 

GAM2C3 = 0.5*GAM2C1/GAM2C2 
GAM2C^ = 0.5*GAM2C2 
FF2C1 = FM2*SQRT(GAM2*G) 

FF2 = FF2C1/(1. +.5*GAM2C2*FM2*FM2)**GAM2C3 
A2 = W2*S0RT(T2*R2)/P2/FF2 
C DEBUG P2,T2,R2.CP2.GAM2C3»FF2C1,FF2,A2 

A2B1 = (1.+GAM2C4*FM2B*FM2B)/( 1.0 + GAM2CA*FM2*FM2) 

A2B = A2*FM2/FM2B*A2B1**GAM2C3/FNB 
02B = S0RT(PIE4*A2B) 

WFWA3 = FULAIR(T2.T3.WFWA1,ETABR) 

WFWA3P = WFwAl + WFWA3 
W3 = W2*( 1.+WFWA3) 

W4 = W3 

R3 = GASC0N(WFWA3P) 

CP3 = SPHTCP(T3,WFWA3P,ETARR) 

GAM3 = CP3/ (CP3-R3/778. lb) 

GAM3C1 = GAM3+1. 

GAM3C2 = GAM3-1. 

GAM3C3 = 0.5*GAM3C1/GAM3C2 
P3P2B = 1.0 - PP2B 
P3 = P2*P3P2B 
C 

FF3 = W3*SQRT(T3*R3 )/P3/A2B 
C 

FM3= F>1M( FF3,GAM3,G) 

C 

A3S1 = (1. + 0.5*GAM3C2*FM3S*FM3S)/( 1.+GAM3C2*FM3*FM3) 

A3S = A2B *FM3/FM3S»A3S1**GAM3C3/FNS 

D3S = SQRT{PIE4*A3S) 

C 

03S = P3*< l.-l./ { 1.+ ( GAM3-1. )*FM3S*FM3S/GAM3) **(GAM/ { GAM-1. )) ) 
P4 = P3 - 0MEGAS*Q3S 
C 

GAMS = 1.35 
GAMSl = GAMS +1. 

GAMS2 = GAMS-1. 

GAMS3 = 0.5*GAMS1/GAMS2 
GAMS4 = GAMS/GAMSl 
GAMS5 = GAMS2/GAMS 
GAMS6 = GAMS/GAMS2 
FFSl = 1. +0.5*GAMS2 *FMS*FMS 
FFS = FMS*SQRT(GAMS»G)/FFS1**GAMS3 
PSD = PAMB 
FFD = FFS/PDPS/ADAS 
C 

FMD= FMMt FFD, GAMS. G) 
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C 


23 PD = PAMB*(1. + 0.5*GAMS2 • FMD*FMD ) **GAM S6 
PS = PD*PDPS 


BEGIN ETAT CALCULATION 


KETA =D 

^3 CALL ADIABTlT4.P4,ETAT,WFWA3P,TS,PS,DHACT,ETAfaB) 
DHACT = ABS(DHACT) 

RS = GASCON ( WFWA3P ) 

GAMS = GAMMA( TS,WFWA3P,ETABR) 

GAMSl = GAMS +1. 

GAMS2 = GAMS-1. 

GAMS3 = 0. 5*GAMS 1/GAM52 

GAMS4 = GAMS/GAMSl 

GAMS5 = GAMS2/GAMS 

GAMS6 = GAMS/GAMS2 

FFSl = 1. +0.5*GAMS2 *FMS*FMS 

FFS = FMS*S0RT(GAMS*G)/FFS1**GAMS3 

PSD = PAM8 

FFD = FFS/PDPS/ADAS 

FMD = FMM( FFD, GAMS, G ) 

33 PD' = PAMB*(1. + 0.5*GAMS2 • FMD*FMD ) ** GAM S6 
PS = PD*PDPS 
PSS = PS/FFS1**GAMS6 

CALL A0IABT(T4,P4,ETAT,WFWA3P,TS,PS,DHACT,ETABR) 
DHACT = ABS(DHACT) 


BEGIN ZLEAK INTERATION 

30 WT = W4*DHACT 

PLl = 1. + 0.5*GAMS2*FMO*FMO 

PL = PD*I1. - WLBAR*( 1. - 1 ./PL1**GAMS6> ) 

FML = S(3RT(2./GAMS2* ( (PL/PAMB)**GAMS5-1. ) ) 

TL = TS 

TSL = TL/(1. + .5*GAMS2*FML*FML) 

VL = FML*SQKT( GAMS*G*RS*TSL ) 

FNT =(W4*( l.-ZLEAK)*VL/G)*CT 
TD = TS 

AS = W4*( l.-ZLEAK)*SQRT(TS»RS)/PS/FFS 
AD = W4*( l.-ZLEAK)*SQRT( TO*RS)/PD/FFD 
PIO = PAMB 

PllO = l.-WiObR*Il. -l./Il. + .2*FM11*FM11)**3.5) 
Pll = Plin • PIU 
P13 = FPR*P11 

CALL ADIABT(T11,P11, ETAF,0.,T13,P13,0HFACT,0. ) 

WF = WT 

DHFACT = ABS(UHFACT) 

WIO = WF/DHFACT 

PSll = Pll/(1. +0.2*FM11*FM11)**3.B 
TSll = Tll/(1. + .2*FM11*FM11) 

Vll = 49. 1*SQRT( TSll )*FM11 
P14 = P1413*P13 

FM14 = SQRT(GAMC3*(P14/PAMB)**GAMC2 - GAMC3) 
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P151^ = 1. -WL6AR*{1. +.2*FM14*FM14)**3.5) 

PS14 = PAMB 
PS15 = PAMB 

P15S15 = P1514*P14/PS14 

FM15 = S0RT(GAMC3*(P15S15)»*GAMC2 -GAMC3) 

FM13 = FMll + FKV*(FM15-FM11 ) 

PS13C = 1. +.2*FM13»FM13 
PS13 = P13/PS13C**3. 5 
TS13 = T13/PS13C 
V13 = A9. 1*FM13*S0RT { TS13) 

T14 = T13 
T15 = T13 

TS14 = T14/( l.+.2*FM14*FM14) 

V14 = ^9. 1*FM14*5CRT ( TS14) 

TS15 = T15/( l. + .2*FM15*FMI‘i) 

V15 = 49. 1*FM15*SQRT ( TSl 5) 

FNF =(W10 *V15/G)*CT 
FN = FNT4-FNF 


FAN AREAS 

All = <J10*R*TS11/PS11/V11 

A13 = rJ10*R*TS13/PS13/Vl 3 

A14 = W10*R*TS14/PS14/V14 

A15 = W10*R*TS15/PS15/V13 

DFT = 1 . 128*SQRT ( All/ ( 1. - RHRT*RHRT>) 

OFR = DFT *KHRT 

OTR = DFT + RING + RING 

XIBET2=0TR 


NEW CARDS ADDED 3-7-69 

127 IF(KETA-l) 120,121»121 

120 AS53 = AS 
GO TO 124 

121 IFHTURB -2) 122,122,123 

122 AS55 = AZ5 
GO TO 124 

123 AS53 = AS 

124 IFITHARC - (360. - .01)) 126,129,129 

TO COMPUTE DTT FOR PARTIAL ADMISSION 

126 DTT = SQRT(360.*PIF4*AS55/THARC + DTR*DTR) 

GO TO 128 


129 DTT =1. 1284*SQRT(DTR*DTR/PIE4 +AS55) 

128 TRR = DTR/OTT 
C 

TBH = 3.5* (DTT -OTR ) 

DT = .5*(DTR + DTT) 

U = UFT*DT/DFT 
ETATSV = ETAT 
IF(IETAT) 133,133,40 
133 GO TO (41,132,52.998) .ITURB 
C 

132 CALL TURB2(T4,W4,P4,PSS, DFT.OT, FMS.ETAN, 

2 F3, FMll, Pll, PAMB, 


113 



c 

c 


c 


3 

4 

5 

6 

7 

8 
9 
1 
2 
3 
1 

GO TO 53 


TS5 ,V5,FM5,ALPH5, VZ5,VU5,P5,U, BETAS, 

R5, FMR5,R6,ETAB,TS6, VZ6,FMZ6,FMR6,BETA6, 

RJ 6, OMEGA, ALPH6, P6,PS,T6,FMSCAL, WLEAK, 

PHI ,PS L0,SPR,ZLEAK,ETATL1, ETAT, 

FLAM,FiNUBR,ETABR, 

CM,FKSB,FKSN,XIALS, 

AZ5 ,FMfe,FMZ5,PSFAN,RE,RU5, TS, TSS, VO , VU6, X I M , 
XIN.XIP.XIS, ZETAB,ZETAM,ZETAP, 

ZETAS,TSP, DHTTI, CP6 , DH , AEXP , DTR , 
TBH.THARC, C00S,AR,F14) 


52 CALL 
1 
1 
2 

3 

4 

5 

6 

7 

8 
9 
1 
2 

3 

4 

5 

6 

7 

8 


TURB3 (T4, W4,F4,PSS, OFT, FMS.ETAN, 

WFWA3P, 

FMlljPll, PAMB,FKS,SWID,FNSML, 

ETABR,FK, 

ALPH5, ALPH6, ALPH7, ALPH8, AZ5, AZ7, AZ8,AZS, 
BETAS, BET A6, BET A7, BETAS, CP6, CPS, DELBTl, 
0ELBT2,DHN1,DHN2,DLHB1,DLHB2,DT,ETAB1, 
ETAB2, ETANI,£TAN2,ETAT,ETATL1, 
F<M,FKM1,FKM2,FM5,FM6,FM7,FM8,FMR5, 
FMR6,FMR7, FMR8,FMZ5,FMZ5SV,FMZ6,FMZ7, 
FMZ7SV,FMZ8,GAM5,GAM6,GAM8,P5,P6, 

P7, PH, PH I, PS, PS5,PS6,PS7,PS8,PS10,R5,R6,R7, 
R8 ,RH0S,RHn6,RH07,RH08,RU5,RU6,RU7, 
RU8,SHP,SPR, 

T6, T8, TS,TSS, TSP, TSS , T S6 , TS7 , TS8 , U , 

VS, V 7, VS, VU5, VU6, VU7, VU8, VZ8, VZ5, 

VZ6 ,VZ7, VZS, WIDEAL.WLEAK, XIBETl, 

XI BET2 ,ZLEAK, 

OMEGA,CM,FKSB,FKSN,XI ALS,TRM,XI AL6,CCM,FLLC J 


GO TO S3 


C 


41 CALL 
2 

3 

4 

5 

6 

7 

8 
9 


TURBl (T4, V<4, P4,PSS, OFT,DT, FMS,ETAN, 

WF(JA3P,FM11,P11,PAMB, 

TS5,V5,FM5,ALPH5,VZ5, VU5 , PS , U , BE TA5 , 

R5 , FMRS, R6 ,ETAB,TS6, VZ6,FMZ6, FMR6, BETAS, 
RJ6,0MEGA, ALPH6, P6,P SC HK,T6, FMSCAL, WLEAK, 
PHI ,PS11,SPR,ZLEAK,ETATL1, ETAT, 

FLAM, FNUBR,ETABR, 

CM,FKSB,FKSH,XIALS,CCM, 


AZ5,FM6, FMZ5,PSFAN,RE,RU5,TS, TSS, VO, VU6,XIM, 
XI N, XI P,XI S,ZETAAS,ZETAB,ZETAM,ZETAP,Z£TASS, 
ZETAS,TSP,DHTTI , CPS, DH, AEXP, DTR, 

TBH, THARC,C00S,AR,F14) 


GO TO 53 
998 CALL TURB4I 
1 
2 

4 

5 
S 

7 

8 
9 


T4, W4, P4.PSS, OFT,DT, FMS,ETAN, 

F3, FM11,P11,PAMB, 

TSS ,V5,FM5,ALPH5, VZ5, VU5,P5,U,BETA5, 

R5 , FMRS,RS,ETAfl, TSS, VZS,FMZS,FMR6,BETA6, 
RUS , OMEGA, ALPH6, PS, PS, TS, FMSCAL, WLEAK, 
PHI ,PS10,SPR,ZLEAK,ETATL1, ETAT, 

FL AM, FNUBRjETABR, 

CM, FKSB, FKSN, XIALS.CCM, 
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1 AZ5,FM6,FMZ5.PSFAN,RE.RU5»TS,TSS.VQ,VU6,XIM, 

2 XIN, XIP, XIS,ZETAAS,ZETAB,ZETAM,ZETAP,ZETASS, 

3 ZETAS,TSP,DHTTI , CP6, DH» AEXP » DTR, 

1 TBH,THARC,C00S,AR,F14) 

53 IFIABS(ETAT-ETATSV) -.001) AO, 40.42 
42 KETA = KETA +1 

IFCKETA -100)43,43,44 

44 WRITE(6,45) ETAT.ETATSV 

45 FQRMATdH ,52H IN MAIN PROGRAM, ETAT CAN NOT BE FOUND AT STATE. 41 

1 /IH ,7H ETAT =E1 2. 5 , 5X , 7HETASV =E12.5) 

40 WFSTR = 14.696*W11*SQRT( T11/518.7)/P11 
IFIIETAT) 61,61,63 
61 GO TO 163,63,64,63), ITURB 

63 IF(IOUT) 75,75,74 

74 WRITE15,67) A1 1 , A13, A 14, A 1 5, A2, A2B, A3S , AD, AS , D28,D3S, 

1 0FR,DFT,DTR,DT,DTT,ETATSV,FM13,FM14,FM15, 

2 FM3,FMD,FML, FN,FNF,FNT,P11,P13,P14,P2,P3, 

3 P4,PL,PD,PS,PSS,PS11,PS13, T14,T2,TBH, 

4 TRR.TSL 

WRI rE15,68 ) TSP.TSll ,TS13, TS 14, TS 15 , U, VI I , VI 3 , V 14 , VI 5 , 

1 VL.W3, WFWA3,WFWA3P,WT,W10 

V6=SQRT(VU6**2+VZ6**2 ) 

WR ITE( 6,69) ALPH5, ALPH6, AZ 5 , BE TA5 , B ETA6 , E T AB , E T AN , F TA T , 

1 ETATLl, FLAM,FM5,FM6,FMR5,FMR6, 

2 FMSCAL,FMZ5, FMZ6,FNUBR, P5,P6,PHI, 

3 PS11,V6 ,PSFAN,R5,R6,RE,RU5,RU6,SPR,T6, 

4 TSb,TS,TSS,V5, VO, VU5, VU6,VZ5,VZ6,WLEAK,XIM,XIN 

WR1TE15,70) XIP.XIS. ZETAB.ZETAM, ZETAP.ZETAS, 

1 ZLEAK,TSP,DHTTI,TS5,CP6,0H 

75 GO TO 147,47,65,47), I TURB 
47 IF(IOJT) 65,65,89 

89 WRITE(6,48) OMEGA , ZETAAS ,ZET ASS 
GO TO 65 

64 IF(IOUT) 65,65,76 

76 WRITE(6,67) A 1 1 , A 13, A 1 4 , A1 5 , A2 , A2 8 , A 35 , AD , AS , 02B , D 3S , 
10FR,DFT,DTR,DT,DTT,ETATSV.FM13,FM14, FM15, 
2FM3,FMD,FML,FN,FNF,FNT,P11,P13,P14,P2,P3, 
3P4,PL,PD,PS,PSS,PS11,PS13,T14,T2,TBH,TRR,TS1 

C 

GAM7=FK 

WRITE! 6,68 ) TSP , TS 11 , T S 1 3, TS 14, TS 1 5 , U , V 11 , V 1 3 , V14, VI 5 , 

IVL, W3, WFWA3, WFWA3P,WT ,Win 
C 

V6=S0RT { VU6**2+VZ6**2 ) 

V8=SQRT(VU8**2 + VZ8**2 ) 

C 

DEBUG T4,W4,ALPH5,ALPH6,ALPH8,AZ5,AZ7,AZ8,AZS,BETA5, 

1BETA6, BET A7, BETAS, CP6, CPS 
C 

DEBUG OHNl, DHN2, ETAfll,ETAB2, ETANl , ETAN2 , E TA T , ETATL I , 
1FM5,FM6,FM7,FM8,FMR5, FMR6,FMR7,FMR8,FMZ5,FMZ6,FMZ7 
C 

DEBUG FMZ8,GAM5,GAM6 , GAMS, P5,P6,P7,P8,PHI,PS5,PS6, 
lPS7,PS8,PS10,R5,R6,R7,R8,RH05,RHn6,RH07,RH08,RU5, 
2RU6,RJ7,RU8,SHP,SPR 

DEBUG T6,T8,TS,rSS,TSP,TS5, TS6, rS7,TS8, V5, V6, V7, V8, VS, 

1VU5 ,VU6, VU7, VU8,VZ8, VZ5, VZ6, VZ7,WIDEAL,WLEAK, 

2ZLEAK,0MEGA,GAM7, ALPH7 

65 WF3P = WFWA3P*W2 
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SFC = <^F3P/FN 
BPR = kJlO/w'4 
C 

CALL ADIABT(Tl,Pl,1.0tFl ,TS1 .PAMH, OHFNG.ETABR) 

DHEMG = ABS(DHENG) 

VENG = SQRT(2.*G*778. 16*DHENG) 

C the CT which multiplied FNENG was REMOVED 3-7-69 

C 

FNE MG = W1*VENG /G 
AUG = FN/FNENG 
C 

IF(IOJT) 77,77,78 

78 WRITE(6,71) WF STR , W F 3P , SFC , BPR , AUG 
C 

C 

77 IF(IWRITE) 94,94,80 
80 IFIlWRITE - 1) 84,84,83 

83 Wl= W1S2 
FNl = FM 

IF(IOUT) 79,79,86 

86 WRITE(6,66) 

WRITE(5,82) W1S2 

79 IWRITE = 1 
GO TO 31 

84 FN2 = FN 

W1S3 = WlSl + IW1S2 - W1S1)*(FN3 - FND/IFN2 - FNl) 

W1 = W1S3 
IWRITE = 0 
IF(IOUT) 81,81,87 

87 WRITE(6,66) 

WRITE(6,85) W1S3 
GO TO 81 

C 

C 

94 IF( lOJT) 95,95,88 

88 CONTINUE 


67 

FORMAT ( IHO, 7H 

All =E12.5,5X,7H 

A13 =E12.5,5X,7H 

A14 

= E12.5 

1 

5X,7h 

A15 =E12.5,5X,7H 

A2 =E12.5/ 



2 

IH ,7H 

A2B =E12. 5, 5X, 7H 

A3S =E12.5,5X,7H 

AD 

=£12.5, 

3 

5X, 7H 

AS =E12.5 

/ 



4 

IH ,7H 

D26 =E12.5,5X,7H 

DBS =E12.5,5X,7H 

UFR 

=£12.5, 

5 

5X, 7H 

DFT =E12.5,5X,7H 

DTR =E12.5/ 



6 

IH ,7H 

DT =E12.5,5X,7H 

DTT =E12.5,5X,7HETATSV 

=E12.5, 

7 

5X, 7H 

FM13 =E12.5,5X,7H 

FM14 =E12.5/ 



8 

IH ,7H 

FM15 =E12.5,5X,7H 

FM3 =E12.5,5X,7H 

FMD 

=£12.5, 

9 

5X,7H 

FML =E12.5,5X,7H 

FN =E12.5/ 



1 

IH , 7H 

FNF =E12.5,5X,7H 

FNT =E12.5,5X,7H 

Pll 

=E12.5, 

2 

5X , /H 

P13 =E12.5,5X,7H 

P14 =E12.5/ 



3 

IH ,7H 

P2 =E 12. 5, 5X ,7H 

P3 =E12.5,5X,7H 

P4 

=E12. 5, 

4 

5X , 7H 

PL =E12.5,5X,7H 

PD =£12.5/ 



5 

IH , 7H 

PS =£12.5, 5X, 7H 

PSS =E12.5,5X,7H 

PSll 

=E12.5, 

6 

5X, 7H 

PS13 =E12.5 

/ 



7 

IH ,7H 

T14 =E12.5,5X,7H 

T2 =E12.5,5X, 7H 

TBH 

=£12.5, 

8 

5X, 7H 

TRR =E12.5,5X,7H 

TSL =E12.5> 



68 1 

FORMATdH ,7H 

TSP =ei2.5,5X,7H 

TSll =E12.5,5X,7H 

TS13 

=E12.5, 

1 

5X, 7H 

TS14 =E12.5,5X,7H 

TS15 =E12.5/ 



2 

IH ,7H 

U =E12.5, 5X,7H 

Vll =E12.5,5X,7H 

V13 

=£12.5, 
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3 

5X,7H 

V14 =E12.5,5X,7H 

V15 =£12.5/ 


4 

IH ,7H 

VL =E12.5,5X,7H 

W3 =£12.5, 


5 

5X, 7HWFWA3 =E12.5, 5X , 7HWFW A.3P = 6 1 2. 5 , 5X , 


6 

7H WT 

=E12.5/1H , 7H WIO =E12.5) 


69 1 

FORMAT! IHO, 7HALPH5 = E 1 2. 5, 5X , 7HALPH6 =E12.5,5X,7H 

AZ5 =£12.5, 

1 

5X.7HBETA5 = E 1 2 . 5, 5X , 7H8ET A6 =E12.5/ 


2 

IH ,7H 

ETAB =E12.5,5X,7H 

ETAN =E12.5,5X,7H 

ETAT =£12.5, 

3 

5X,7HETATL1=E12.5/ 



4 

IH ,7H 

FLAM =E12.5,5X,7H 

FM5 =E12.5,5X,7H 

FM6 =£12.5, 

5 

5X,7H 

FMR5 =E12.5,5X,7H 

FMR6 =E12.5/ 


6 

IH ,7HFMSCAL=E12.5,5X,7H 

FMZ5 =E12.5,5X,7H 

FMZ6 =£12.5, 

7 

5X,7HFNUBR =£12.5 

/ 


8 

IH ,7H 

P5 =E12.5,5X,7H 

P6 =E12.5,5X,7H 

PHI =E12.5, 

9 

5X, 7H 

PSll =E12.5,5X,7H 

V6 =£12.5/ 


1 

IH ,7HPSFAN =E12.5,5X,7H 

R5 =E 1 2. 5 , 5x , 7H 

R6 =E12.5, 

2 

5X,7H 

RE =E12.5,5X,7M 

RU5 =£12.5/ 


3 

IH ,7H 

RU6 =E12.5,5X,7H 

SPR =E12.5,5X,7H 

T6 =E12.5, 

4 

5X, 7M 

TS6 =E12.5,5X,7H 

TS =£12.5/ 


5 

IH ,7H 

TSS =£12.5,5X,7h 

V5 =E12.5,5X,7H 

VO =E12.5, 

6 

5X, 7H 

VU5 =E12.5,5X,7H 

VU6 =£12.5/ 


7 

IH , 7H 

VZ5 =E12.5,5X,7H 

VZ6 =E12.5, 5X,7HWLEAX =E12.5, 

8 

5X ,7H 

XIM =E12.5,5X,7H 

XIN =E12.5) 
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FORMAT! IH ,7H 

XIP =E12.5,5X,7H 

XtS =E12.5,5X,7HZ£TA8 =E12.5, 

1 

5X,7HZETAM =E12.5 

/ 


2 

IH ,7HZETAP =E12. 5, 5X.7HZETAS =E12.5, 


3 

5X,7HZLEAK =£12.5, 




^ 5X,7H TSP =E12.5/ 

5 IH ,7H0HTTI =F12.5,5X,7H TS5 =E12.5,5X, 

6 7H CP6 =E12.5,‘SX,7H DH =E12.5) 

71 FORMAH IH ,7HWFSTR =E12.5,5X,7H WF3P =E12.5,5X,7H SFC =E12.5, 

1 5X,7H BPR =E12.3,5X,7H AUG =E12.3) 

66 FORMAT( IHK) 

82 FORMATl IHO.yH Wl (2) =E12.3) 

85 FORMAT( lH0,9h Wl (3) =E12.5) 

106 FORMArtlHl, 7H FMS = E 1 2 . 5 , 5 X , 7HP 20P 1 =E12.5,5X,7H Wl =E12.5, 

1 5X,7H FN3 =£12.5) 

A8 FOR'^AKIH f7HUMEGA = £ 1 2 . 5, 4X , 8HZE T A A S = E 1 2 . 5 1 4 X f 8H ZE T A SS =E12.5) 
95 RETUR'J 
END 
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APPENDIX I 


PROGRAM NOMENCLATURE 
Program Inputs 


FORTRAN 

Engineering 

Definition 

symbol 

symbol 


AR 

ARrp 

turbine blade aspect ratio 

AEXP 

b 

exponent in seal leakage equation 

CCM 

CC 

turbine rotor axial velocity multiplier 

CM 

^^sh 

turbine blade Mach number correction coefficient 

CT 

Cp 

thrust coefficient (turbine and fan) 

COOS 

^T 

turbine blade solidity 

F14 

^pa 

partial-admission loss coefficient 

HV 

HV 

lower heating value of fuel, Btu/lbm 

H2P 

H° 

liquid fuel inlet enthalpy, Btu/lbm 

OMEGAS 

Ksc 

scroll-loss coefficient 

FKSB 

(I^se)B 

rotor blade secondary- loss coefficient 

FKSN 

(I^se)N 

stator secondary- loss coefficient 

FKV 

^F 

fan internal loss factor 

TLH 

(L/H)j) 

turbine exit duct length to height ratio 

FLLC 

X 

two-stage turbine work split factor 

FM2 

M2 

inlet duct exit Mach number 

FM2B 

^2b 

burner inlet flow Mach number 

FM3S 

^3sc 

scroll inlet flow Mach number 

FMS 

Mg 

turbine exit stator exit Mach number 

FMll 

^10 

fan inlet axial Mach number 

FNS 

*^A 

number of scroll arms per inlet duct 

FNB 

% 

number of burners 

PI 

Pi 

inlet duct inlet total pressure, psia 
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FORTRAN Engineering Definition 

symbol symbol 


P20P1 

P 2 /P 1 

inlet duct pressure ratio 

PAMB 

Po 

ambient pressure, psia 

FPR 

Pis/Ph 

fan-stage pressure ratio 

P1413 

Pi4/Pl3 

fan outlet duct pressure ratio 

PP2B 

Ap’/p^ 

burner pressure drop fraction 

R 

R 


air gas constant, 53.35 ft-lbf/(lbm)(°R) 

RHRT 

Y 

F 

fan hub-tip ratio at rotor inlet 

RING 


fan-shroud carrier height, in. 

TAMB 

T, 

0 

ambient temperature, °R 

TRM 

TRM 

two- stage turbine mean line radius change factor 

T1 

T 

t 

1 

inlet duct inlet total temperature, °R 

T4 

T 

f 

4 

turbine inlet total temperature, °R 

DTLOSS 

AT^ 

inlet duct temperature loss, °R 

UFT 

U 

F,t 

fan tip speed, ft/sec 

W1 

w 

1 

inlet duct weight flow, Ibm/sec 

WFWAl 

f 

a 

1 

fuel-to-air ratio at inlet duct 

ZIL 


1 

3 

interstage leakage ratio 

ETABR 

^b 

burner efficiency 

ETAF 

'^F 

fan- stage adiabatic efficiency 

ETAN 


initial value of inlet stator efficiency 

ETAT 

^7rp 

initial value of turbine total-to-total efficiency 

WIOBR 

CO. 

10 

fan inlet loss coefficient 

WLBAR 

CO. 

L 

louver loss coefficient 

THARC 

e 


turbine admission arc, deg 

XIAI5 


inlet or interstage stator profile loss 
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Program Outputs - Main Program 


FORTRAN 

Engineering 

Definition 

symbol 

symbol 


All 

All 

2 

fan flow area at rotor inlet, in. 

A13 

Ai3 

2 

fan flow area at stator exit, in. 

A14 

Ai4 

2 

fan flow area at duct exit, in. 

A15 

Ai5 

2 

fan flow area at louver exit, in. 

A2 

^2 

2 

turbine flow area at inlet duct exit, in. 

A2B 

A 2 b 

2 

turbine flow area at burner inlet, in. 

ASS 

Assc 

2 

turbine flow area at scroll inlet, in. 

AD 

Ad 

2 

turbine flow area at exit duct, in. 

AS 

As 

turbine flow area at exit stator exit, in. ‘ 

D2B 

^ 2 b 

burner diameter, in. 

D3S 

^3sc 

scroll inlet diameter, in. 

DFR 

‘*F,h 

fan rotor diameter at hub, in. 

DFT 

*^F,t 

fan rotor diameter at tip, in. 

DTR 

^T,h 

rotor diameter at turbine hub, in. 

DT 

'^T.m 

rotor diameter at midturbine, in. 

DTT 


rotor diameter at turbine tip, in. 

FN 

F 

calculated total thrust, Ibf 

FNF 

^F 

calculated fan thrust, Ibf 

FNT 

Fp 

calculated turbine thrust, Ibf 

FM13 

^13 

Mach number at fan stator exit 

FM14 

Mi4 

Mach number at fan duct exit 

FM15 

^15 

Mach number at fan louver exit 

FM3 

M 3 

Mach number at burner exit 

FMD 

Mj3 

Mach number at turbine duct exit 

FML 


Mach number at turbine louver exit 

Pll 

Pil 

total pressure at fan rotor inlet, psi 


120 



FORTRAN 

Engineering 

Definition 

symbol 

symbol 


P13 

Pl3 

total pressure at fan stator exit, psi 

P14 

PU 

total pressure at fan duct exit, psi 

P2 

P2 

total pressure at inlet duct exit, psi 

P3 

P3 

total pressure at burner exit, psi 

P4 

P4 

total pressure at turbine inlet, psi 

PL 

Pl 

total pressure at turbine louver exit, psi 

PD 

Pb 

total pressure at turbine duct exit, psi 

PS 

Ps 

total pressure at turbine exit stator exit, psi 

PSS 

Ps 

static pressure at turbine exit stator exit, psi 

PSll 

Pll 

static pressure at fan rotor inlet, psi 

PS13 

Pl3 

static pressure at fan stator exit, psi 

T2 


total temperature at inlet duct exit, °R 

T14 

Ti4 

total temperature at fan duct exit, °R 

TBH 

Hb 

turbine blade height, in. 

TRR 

Yt 

turbine hub-tip ratio 

TSL 

Tl 

static temperature at turbine louver exit, °R 

TSll 

Til 

static temperature at fan rotor inlet, °R 

TS13 

Ti3 

static temperature at fan stator exit, °R 

TS14 

Ti4 

static temperature at fan duct exit, °R 

TS15 

Ti5 

static temperature at fan louver exit, °R 

U 

. 

turbine midwheel speed, ft/sec 

Vll 

Vll 

axial velocity at fan rotor inlet, ft/sec 

VI 3 

Vl3 

axial velocity at fan stator exit, ft/sec 

V14 

Vl4 

axial velocity at fan duct exit, ft/sec 

VI 5 

Vl5 

axial velocity at fan louver exit, ft/sec 

VL 

Vl 

axial velocity at turbine louver exit, ft/sec 

W3 

^3 

mass flow at burner exit, Ibm/sec 
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FORTRAN 

Engineering 

Definition 

symbol 


symbol 


WFWA3 

f 

a 

b 

fuel-to-air ratio added in burner 

WFWA3P 

a 

4 

total fuel-to-air ratio at burner exit 

WT 

P 

T 

turbine power output, Btu/sec 

WIO 

w 

10 

fan air flow rate, Ibm/sec 



Program Output - Turbine Subroutine 

AUG 

AUG 

thrust augmentation ratio 

AZ5 

A 

5 

2 

rotor inlet annular flow area, in. 

BPR 

BPR 

bypass ratio 

CP6 

C 

p,6 

specific heat at rotor exit, Btu/(lb)(°R) 

SD 

D 


exit stator diffusion factor 

DH 

<1 

inlet stator actual enthalpy drop, Btu/lbm 

DHTTl 

id 

isentropic total-to-total enthalpy drop, Btu/lbm 

FM5 

M 

5 

absolute rotor inlet Mach number 

FM6 

M 

6 

absolute rotor exit Mach number 

FMR5 

M 

r, 5 

relative rotor inlet Mach number 

FMR6 

M, 

r,6 

relative rotor exit Mach number 

FMSCAL 

M, 

1 

S (calculated) 

stator exit Mach number (calculated) 

FMZ5 


axial rotor inlet Mach number 

FMZ6 

^x,6 

axial rotor exit Mach number 

P5 

P’5 


total pressure at rotor inlet, psi 

P6 

P6 


total pressure at rotor exit, psi 

PSll 

Pll 

static pressure at fan rotor inlet, psi 

PS FAN 

Pex 

static pressure at bellmouth (seal exit), psi 
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FORTRAN 

symbol 

Engineering 

symbol 

Definition 

R5 

W5 

relative velocity at rotor inlet, ft/sec 

R6 

Wg 

relative velocity at rotor exit, ft/sec 

RE 

Re 

Reynolds number 

RU5 

"^u,5 

tangential component of relative velocity at rotor 
inlet, ft/sec 

SFC 

SFC 

specific fuel consumption of tip-turbine fan, 
Ib-fuel/hr 
lb-thrust 

SPR 

SPR 

seal pressure ratio 

T6 


absolute total temperature at rotor exit, °R 

TS6 

^6 

static temperature at rotor exit, °R 

TS 

'^S 

total temperature at exit stator exit, °R 

TSP 

T^,id 

isentropic total temperature at exit stator exit, °R 

TSS 

Ts 

static temperature at exit stator exit, °R 

TS5 

T5 

static temperature at rotor inlet, °R 

V5 

V5 

absolute velocity at rotor inlet, ft/sec 

V6 

Ve 

absolute velocity at rotor exit, ft/sec 

VO 


ideal jet velocity, ft/sec 

VU5 


tangential component of the absolute velocity at 
rotor inlet, ft/sec 

VU6 

Ve 

tangential component of the absolute velocity at 
rotor exit, ft/sec 

VZ5 

^X,5 

axial component of velocity at rotor inlet, ft/sec 

VZ6 

^x,6 

axial component of velocity at rotor exit, ft/sec 

WFSTR 

^eq 

fan corrected mass flow rate, Ibm/sec 

WF3P 

^f,4 

total fuel flow, Ibm/sec 

WLEAK 

^leak 

seal leakage flow rate, Ibm/sec 

ZLEAK 

Z 

ratio of inlet seal leakage to turbine flow rate 
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FORTRAN 

Engineering 

Definition 

symbol 

symbol 


ALPH5 

“5 

rotor absolute inlet angle, rad 

ALPH6 

“6 

rotor absolute exit ar^le, rad 

BETAS 

^5 

rotor relative inlet angle, rad 

BETAS 

^6 

rotor relative exit ar^le, rad 

ETAB 

^B 

turbine blade row efficiency 

ETAN 


inlet stator blade row efficiency 

ETAT 

?7ip 

turbine total-to-total efficiency 

ETATLl 

leak 

turbine total-to-total efficiency (with leakage) 

FLAM 

X 

speed-work parameter 

FNUBR 

V 

blade-jet speed ratio 

OMEGA 


stator total pressure-loss coefficient 

PHI 

e 

seal leakage expansion function 

ZETASS 


stator total-loss coefficient 

XIS 

se^N 

inlet stator secondary-loss coefficient 

ZETAB 

^B 

rotor total-loss coefficient 

ZETAM 

sh^B 

rotor shock-loss coefficient 

ZETAP 

^^p^B 

rotor profile- loss coefficient 

ZETAS 

se^B 

rotor secondary- loss coefficient 

XIM 

sh^N 

inlet stator shock-loss coefficient 

XIN 


inlet stator total-loss coefficient 

XIP 


inlet stator profile- loss coefficient 


Additional Output for Two-Stage Turbine 

AZ7 

A? 

2 

second- stage rotor inlet annular flow area, in. 

AZ8 

-^8 

2 

second- stage rotor exit annular flow area, in. 

AZS 


2 

exit stator exit annular flow area, in. 
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FORTRAN 

symbol 

Engineering 

symbol 

Definition 

CPS 


specific heat at second- stage rotor exit, 
Btu/(lbm)(°R) 

DHNl 

^^N,5 

actual enthalpy drop across first-stage stator, 
Btu/lbm 

DHN2 

7 

actual enthalpy drop across second-stage stator, 
Btu/lbm 

FM7 

M,^ 

absolute second- stage rotor inlet Mach number 

FM8 

Mg 

absolute second-stage rotor exit Mach number 

FMR7 

^r,7 

second-stage relative rotor inlet Mach number 

FMR8 

00 

second-stage relative rotor exit Mach number 

FMZ7 

^x,7 

second- stage rotor axial inlet Mach number 

FMZ8 

^x,8 

second-stage rotor axial exit Mach number 

P7 

P? 

total pressure at second rotor inlet, psi 

P8 

P8 

total pressure at second rotor exit, psi 

PS5 

P 5 

static pressure at first rotor inlet, psi 

PS6 

Pe 

static pressure at first rotor exit, psi 

PS7 

P 7 

static pressure at second rotor inlet, psi 

PS8 

P8 

static pressure at second rotor exit, psi 

R7 

w, 

second-stage rotor inlet relative velocity, ft/sec 

R8 

Wg - 

second-stage rotor exit relative velocity, ft/sec 

RU7 

^^u,7 

second-stage rotor inlet tangential component of 
relative velocity, ft/sec 

RU8 

00 

second-stage rotor exit tangential component of 
relative velocity, ft/sec 

T8 


absolute total temperature at second rotor exit, °R 

TS7 

^7 

static temperature at second rotor inlet, °R 

TS8 

Tg 

static temperature at second rotor exit, °R 

V7 

V7 

absolute velocity at second rotor inlet, ft/sec 
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FORTRAN 

symbol 

Engineering 

symbol 

Definition 

V8 

Vg 

absolute velocity at second rotor exit, ft/sec 

VS 

Vs 

stator exit velocity, ft/sec 

VU7 

^u,7 

second-stage rotor inlet tangential component of 
absolute velocity, ft/sec 

VU8 


second-stage rotor exit tangential component of 
absolute velocity, ft/sec 

VZ7 

Vx,7 

second-stage rotor inlet axial velocity, ft/sec 

VZ8 

00 

second-stage rotor exit axial velocity, ft/sec 

WIDEAL 

^T,id 

ideal stage power, Btu/sec 

SHP 

P, 

actual power output, Btu/sec 

ALPH7 

arj 

second-stage rotor absolute inlet angle, rad 

ALPH8 

“8 

second-stage rotor absolute exit angle, rad 

BETA7 

Pij 

second-stage rotor relative inlet angle, rad 

BETA 8 


second- stage rotor relative exit angle, rad 

ETABl 

^B,6 

first-stage rotor blade row efficiency 

ETAB2 

^B, 8 

second- stage rotor blade row efficiency 

ETANl 

5 

first-stage stator efficiency 

ETAN2 

%, 7 

second-stage stator efficiency 

GAMS 

‘J's 

ratio of specific heats, first rotor inlet 

GAM6 


ratio of specific heats, first rotor exit 

GAM7 


ratio of specific heats, second rotor inlet 

GAM8 

'>'8 

ratio of specific heats, second rotor exit 

RH05 

P5 

O 

first-stage rotor inlet density, Ibm/ft 

RH06 

Pb 

3 

first-stage rotor exit density, Ibm/ft 

RH07 

Prj 

o 

second-stage rotor inlet density, Ibm/ft 

RH08 

Pb 

g 

second-stage rotor exit density, Ibm/ft 
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Figure 1. -Single-stage tip-turbine fan model used in analysis showing major compo- 
nents and station designations. 




Scroll -7 

/ 


Support frame- 








START 



Final output for 
desired thrust 


Figure 4. - Main program iteration flow chart. 
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STflTCIR EXIT MRCH NO.. MS 

Figure 5. - Typical turbine parametric data plot. 
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O Single stage with impulse exit stators 

Single stage without exit stator 
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Exit stator exit axiai fAach number, Nlj 


Figure 6. - Variation of inlet duct flow rate with exit stator 
exit axial Mach number. Total thrust, 10 000 pounds; 
inlet duct inlet pressure, 117.6 psia; turbine inlet tem- 
perature, 1900° R; inlet stator exit angle, ?75°. 
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Figure 7. - Comparison of louver exit velocity for four 
turbine stage configurations. Total thrust, 10 (X)0 
pounds; inlet duct inlet pressure, 117.6 psia; turbine 
inlet temperature, 1900" R; inlet stator exit angle, 
<75°. 



Exit stator exit axial Mach number, Mj 


Figure 8. - Variation of turbine tip diameter with exit stator exit axial 
Mach number. Total thrust, 10 000 pounds; inlet duct inlet pres- 
sure, U7.6psia; turbine inlet temperature. 1900° Rj inlet exit 
angle, <75°. 
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Figure 11. - Fan rotor schematic. 



Figure 12. - Single-stage tip-turbine velocity diagrams. 
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Figure 13. - Two-stage tip-turbine velocity diagrams. 
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Figure 23. - Variation of exponent in equation (F5t for expansion function, single 
element seal. 
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